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Cell  cycle  dynamics  and  the  physiology  of 
saxitoxin  biosynthesis  in 
Alexandrium  fundyense  (Dinophyceae) 

by 

Caspar  Taroncher  Oldenburg 

Submitted  in  partial  fiilfillment  of  the  requirements  for  the  degree  of 
Doctor  of  Philosophy 


ABSTRACT 

The  mechanism  of  saxitoxin  (STX)  biosynthesis  in  marine  dinoflagellates  of  the  genus 
Alexandrium  is  still  unknown.  The  aim  of  this  thesis  was  to  analyze  novel  aspects  of 
toxigenesis  during  the  cell  cycle  in  Alexandrium  and  to  apply  molecular  techniques  to  gain 
new  insights  on  the  genetics  and  regulation  of  STX  biosynthesis. 

Synchronized  cultures  of  A.  fundyense  were  studied  to  determine  the  dynamics  of  toxin 
production  throughout  the  cell  cycle.  Toxin  production  was  discontinuous,  was  induced 
by  light  and  always  occured  during  a  period  of  approximately  eight  to  ten  hours  in  early 
Gi.  Analysis  of  the  cell  cycle  dynamics  suggests  the  existence  of  two  transition  points:  one 
at  the  beginning  of  Gj,  which  is  light-dependent  and  holds  the  cells  in  a  Go-like  period, 
and  a  second  one  at  the  end  of  Gj,  which  is  size-dependent  and  arrests  the  cells  in  Gj.  A 
model  of  the  cell  cycle  of  A.  fundyense  is  proposed  in  which  progression  through  the  cell 
cycle  can  be  arrested  at  two  different  transition  points  located  in  Gj  and  toxin  production  is 
induced  by  light  during  Gj . 
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The  effects  of  temperature  and  phosphate  limitation  on  the  linkage  between  changes  in 
the  duration  of  the  cell  cycle  stages  and  toxicity  were  studied  in  semi-continuous  cultures  of 
A.  fundyense.  A  direct  correlation  between  Gi  duration  and  toxin  content  was  observed, 
along  with  a  clear  uncoupling  of  toxin  accumulation  from  the  S  and  G2  phases  of  the  cell 
cycle.  In  both  experiments,  toxin  production  rates  remained  constant  for  the  respective 
range  of  conditions,  implying  that  the  variations  in  toxin  content  observed  were  a  result  of 
increasing  periods  of  biosynthetic  activity.  Phosphate  limitation  enhanced  toxin  production 
rates  and  affected  interconversions  among  STX  derivatives  in  several  ways:  oxidations  to 
yield  the  hydroxy-series  of  STXs  were  phosphate-dependent  while  sulfatation  reactions 
were  not. 

Differential  Display  (DD)  analysis  was  applied  to  the  identification  of  genes  that  were 
up-  or  downregulated  during  toxigenesis  in  synchronized  cultures  of  A.  fundyense.  Three 
genes  were  isolated:  i'-adenosylhomocysteine  hydrolase,  methionine  aminopeptidase  and  a 
histone-like  protein.  None  could  be  directly  correlated  to  toxigenesis  but  instead  relate  to 
general  cellular  metabolism. 
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INTRODUCTION 


Saxitoxins  in  dinoflagellates: 
Origin,  physiology  and  molecular  biology 


The  increase  in  occurrence  of  harmful  algal  blooms  and  their  rapid  dispersal  into  new 
areas  of  the  world  over  the  past  decades  has  prompted  research  aimed  at  describing  this 
phenomenon  and  ultimately  understanding  the  biological  and  evolutionary  significance  of 
toxin  production  by  phytoplankton.  The  division  Dinoflagellata  (Greek  dinos,  whirling) 
(Fensome  et  al,  1993),  includes  some  of  the  geographically  most  widely  distributed 
species  of  toxic  microalgae.  Some  of  the  common  characteristics  of  these  toxic 
dinoflagellates  are  bi-flagellation,  a  permanently  condensed  nucleus,  a  neritic  or  estuarine 
distribution,  and  the  production  of  toxic  secondary  metabolites  (Taylor,  1987). 

One  of  these  groups  of  compounds  comprises  saxitoxin  (STX),  a  positively  charged 
perhydropurine  alkaloid,  and  its  over  twenty  known  derivatives  (Fig.l)  (Shimizu,  1993). 
STX  is  a  potent  neurotoxin  that  actively  blocks  Na^ -channels  and  consequently  causes 
paralysis  of  the  neuromuscular  system  (Narahashi,  1988;  Lipkind  and  Fozzard,  1994), 
resulting  in  a  symptomatology  called  Paralytic  Shellfish  Poisoning  (PSP)  (Meyer  et  al, 
1928).  PSP  outbreaks  occur  in  littoral  waters  throughout  the  world  and  have  a  major 
impact  on  public  health  and  industry  (Anderson,  1989).  Their  effects  range  from  marine 
mammal,  fish  and  seabird  mortality  to  human  intoxication  and  fatalities  as  a  result  of 
ingesting  contaminated  shellfish  or  fish  (Geraci  et  al,  1989;  Shumway,  1990).  PSP  toxins 
also  induce  alterations  of  the  marine  trophic  balance  through  adverse  effects  on  larvae  and 
other  life  history  stages  of  fisheries  species  (White  et  al,  1989;  Montoya  et  al,  1996). 


ORIGIN  OF  SAXITOXIN 

PSP  was  first  related  to  a  planktonic  dinoflagellate  in  the  genus  Gonyaulax  in  1937 
(Sommer  et  al,  1937).  It  was  proposed  that  the  toxic  phytoplanktonic  cells  were  ingested 
by  filter-feeders  that  were  not  affected  by  the  active  compounds,  resulting  in  a  net 
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accumulation  of  SIX  and  amplification  of  the  toxicity  toward  higher  elements  of  the  food 
chain.  Since  then,  other  species  of  dinoflagellates  in  the  genera  Alexandrium 
(=Gonyaulax,  Protogonyaulax)  (Schantz  et  al,  1975;  Cembella  et  al,  1987;  Steidinger  and 
Moestrup,  1990),  Pyrodinium  (Maclean,  1973;  Harada  et  al.,  1982)  and  Gymnodinium 
(Estrada  et  al.,  1984;  Oshima  et  al.,  1987)  have  been  shown  to  also  produce  saxitoxins.  In 
addition  to  these  dinoflagellates,  several  other  microorganisms  such  as  the  cyanobacteria 
Aphanizomenonflos-aquae  and  Anabaena  circinalis  and  the  bacteria  Moraxella  sp..  Bacillus 
sp.,  Acinetobacter  sp.  and  Pseudomonas  sp.  are  also  able  to  synthesize  saxitoxins  (Alam  et 
al.,  1978;  Kodama  et  al,  1988;  Ogata  et  al,  1990;  Humpage  et  al,  1994;  Franca  et  al, 
1996;  Levasseur  et  al,  1996).  However,  toxin  production  by  bacteria  remains  a 
controversial  issue.  A  variety  of  methods,  including  HPLC,  mouse  neuroblastoma  assay 
and  ELISA,  have  been  used  to  determine  the  presence  of  SIX  in  bacteria,  but  unequivocal 
proof  by  mass  spectrometric  analysis  is  still  needed  to  confirm  the  chemical  structure  of  the 
compounds  involved  (Gallacher  et  al,  1996). 

The  fact  that  some  of  the  above  mentioned  bacteria  can  be  found  in  close  association 
with  dinoflagellates,  either  as  endosymbionts  or  extracellularly,  has  raised  a  controversy  on 
their  possible  role  in  the  production  of  toxins  generally  attributed  to  the  dinoflagellates 
(Rausch  de  Traubenberg  and  Lassus,  1991).  This  so  called  'Bacterial  Hypothesis'  has 
been  substantiated  by  a  suite  of  observations.  The  presence  of  intracellular  bacteria  in 
Alexandrium  and  their  role  in  STX  production  were  first  documented  in  1982  (Silva,  1982; 
Singh  et  al ,  1982).  Deteetion  of  STX  in  the  water  column  in  the  absence  of  dinoflagellates 
has  been  correlated  to  the  presence  of  bacteria  in  the  genera  Moraxella,  Pseudomonas, 
Aeromonas  and  Acinetobacter  (Kodama  et  al,  1990;  Sakamoto  et  al,  1992;  Levasseur  et 
al,  1996).  Laboratory  studies  with  two  bacteria  isolated  from  Alexandrium  cultures, 
Moraxella  sp.  and  Pseudomonas  sp.,  have  shown  that  saxitoxin  accumulation  patterns 
under  nutrient  stress  in  bacteria  follow  similar  trends  to  those  observed  in  the  dinoflagellate 
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Fig.  1  Structures  of  the  saxitoxin  (SIX)  and  neosaxitoxin  (NEO)  series  of  toxins. 

(NEO  is  the  l-AT-hydroxy  derivative  of  STX  (OH  underlined)).  The  basic 
structure,  13-deoxy-decarbamoyl-STX,  is  determined  by  atoms  1  to  13  (compare 
Fig.  3;  STX  is  the  6-carbamoyl  derivative  of  13-deoxy-decarbamoyl-STX  (see 
text)).  Ri  represents  the  ll-P(a)0-sulfatation  site,  R2  represents  the  13- 
hydroxy,  13-0-carbamoyl  or  A^-sulfo-13-O-carbamoyl  substitution  site. 
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under  similar  limiting  conditions  (Kodama,  1990;  Doucette  and  Trick,  1995;  Franca  et  al, 
1996).  However,  in  none  of  the  documented  cases  of  putative  toxin  production  by  bacteria 
did  the  toxin  levels  measured  reach  values  that  could  account  for  toxin  concentrations 
typically  observed  in  bloom  situations  or  in  culture  (Gallacher  et  al,  1996).  Moreover, 
toxin  production  in  axenic  dinoflagellate  cultures  has  been  widely  demonstrated  (Singh  et 
al.,  1982;  Boczar  et  al,  1988;  Kim  et  al,  1993),  while  addition  of  Moraxella  to 
Alexandrium  tamarense  cultures  has  also  been  observed  to  result  in  an  actual  decrease  in 
toxicity  of  the  dinoflagellate  (Ogata  et  al,  1996).  One  additional  argument  against  the 
bacterial  hypothesis  comes  from  dinoflagellate  cyst  mating  experiments  which  suggest  that 
inheritance  of  both  the  ability  to  synthesize  toxin  and  the  toxin  profiles  themselves  follows 
a  1:1  segregation  pattern  for  the  FI  generation.  This  rules  out  a  direct  involvement  of 
bacteria  in  the  process,  as  this  would  result  in  a  random  distribution  of  toxic  phenotypes  in 
the  FI  generation  (Oshima  et  al,  1993;  Oshima,  1995).  All  this  evidence  suggests  that 
both  bacteria  and  dinoflagellates,  as  well  as  cyanobacteria,  are  capable  of  producing 
saxitoxin  but  the  scope  of  their  interactions  and  the  evolutionary  process  underlying  the 
acquisition  of  this  capability  remain  to  be  clarified. 

BIOCHEMISTRY  OF  SAXITOXIN 

Over  4,000  alkaloids  have  been  described  in  organisms  throughtout  the  different 
kingdoms  and  extensive  work  has  been  done  on  elucidating  their  possible  precursors  and 
putative  anabolic  pathways  (Leete,  1969).  Saxitoxin  is  a  positively  charged, 
perhydropurine-like  alkaloid.  Saxitoxin  derivatives  represent  different  combinations  of  1- 
A-hydroxy,  ll-a-P-O-sulfate,  13-hydroxy,  13-0-carbamoyl  and  A-sulfo-13-O-carbamoyl 
substitutions  of  the  basic  skeleton,  13-deoxydecarbamoyl-STX,  that  cover  a  broad  range  of 


toxicities  depending  on  their  affinity  to  Na^-channels  (Fig.  1)  (Hall,  1982;  Shimizu, 
1993).  A  putative  biosynthetic  pathway  for  STX  has  been  proposed  following  a  series  of 
feeding  experiments  with  labeled  compounds  (Shimizu,  1993).  Given  the  purine-like 
structure  of  STX,  a  first  approach  consisted  in  providing  the  dinoflagellate  Gonyaulax 
catenella  (=Alexandrium  catenella)  with  known  precursors  of  the  purine  biosynthetic 
pathway,  but  no  labeling  of  STX  was  observed  (Proctor  et  al,  1975).  Further  attempts  to 
use  dinoflagellates  for  these  studies  were  hampered  by  the  cell's  autotrophic  nature,  which 
makes  them  incapable  of  directly  utilizing  exogenous  organic  molecules,  thus  making  it 
impossible  to  trace  the  fate  of  labeled  compounds  supplied  externally  (Shimizu  et  al, 
1984).  Consequently,  another  organism,  the  cyanobacterium  Aphanizomenon  flos-aquae 
was  used  for  feeding  experiments  and  other  precursors  were  tried  until  it  was  eventually 
determined  that  the  perhydropurine  skeleton  of  STX  incorporates  one  molecule  each  of  the 
basic  aminoacid  arginine  (Arg)  and  acetate,  one  methyl  group  from  methionine,  via  S- 
adenosyl  methionine,  and  three  guanido/ureido  groups  derived  from  Arg  (Shimizu  et  al, 
1984).  This  incorporation  pattern  precluded  a  purine  origin  of  the  perhydropurine-like 
skeleton  of  STX.  Based  on  the  above  scheme,  a  new  biogenetic  pathway  was  proposed, 
consisting  of  a  Claisen-type  condensation  of  Arg  and  an  acetate  molecule,  followed  by 
transamidination  of  a  guanido  group  from  a  second  molecule  of  Arg  with  subsequent 
cyclation  and  methylation  at  the  C-6  position  (Fig  2)  (Shimizu,  1993). 

The  sequence  of  transformations  among  derivatives  can  only  be  inferred  from  indirect 
evidence.  Predominance  of  11-P-epimers  (GTX  III,  GTX  IV,  C2  and  C4)  over  11-a- 
epimers  (GTX  I,  GTX  II,  Cl  and  C3)  suggests  that  only  11-P-epimers  are  synthesized  in 
Alexandrium  and  Gymnodinium  and  that  the  presence  of  11-a-epimers  is  a  result  of  post¬ 
sampling,  in  vitro  epimerization  (Hall  et  al,  1990;  Oshima  et  al,  1993).  Characterization 
of  several  V-sulfotransferases  in  these  two  species  that  use  exclusively  l-V-dehydroxy 
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Fig.  2  Biogenetic  pathway  proposed  for  the  biosynthesis  of  SIX  (reviewed  by 
Shimizu,  1993).  Labeling  studies  were  performed  with  the  dinoflagellate 
Gonyaulax  catenella  (=Alexandrium  catenella)  and  the  cyanobacteria 
Aphanizomenon  flos-aquae.  Arginine,  acetate  and  S-adenosyl-methionine 
(SAM)  have  been  shown  to  be  the  direct  precursors  of  STX. 
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derivatives  as  substrates  (GTX  II,  GTX  III  and  STX)  (Oshima,  1995;  Yoshida  et  ah, 
1996),  and  the  observation  of  a  specific  1-A^-oxidase  in  Alexandrium  that  converts  GTX  II 
and  GTX  HI  to  GTX  I  and  GTX  IV  respectively  (Oshima,  1995),  further  suggest  that  l-N- 
oxidation  is  a  post-sulfatation  process.  A  combination  of  the  above  information,  toxin 
compositions  of  different  organisms  and  the  suite  of  derivatives  that  have  been  isolated 
from  STX  producing  organisms,  allows  us  to  start  assembling  a  picture  of  how  the 
multiple  STX  derivatives  could  be  biochemically  interconnected  (Fig.  3).  New  insights  on 
the  enzymes  responsible  for  the  interconversions  will  be  gained  from  more  detailed  labeling 
studies  and  analysis  of  changes  in  toxin  composition  over  time,  as  toxin  profiles  are  a 
reflection  of  enzyme  activities  involved  in  those  pathways  (Oshima,  1995). 


PHYSIOLOGY  OF  SAXITOXIN  PRODUCTION 

Analyzing  the  variability  in  toxicity  that  can  occur  with  different  growth  conditions  or 
between  different  isolates  of  the  same  species  is  of  fundamental  importance  to  an 
understanding  of  saxitoxin  biosynthesis.  Toxicity  can  be  studied  at  two  different  levels: 
Changes  in  'toxin  content',  i.e.,  in  the  overall  toxicity  of  a  cell  resulting  from  the  integrated 
potency  of  all  toxins  present,  or  variations  in  'toxin  composition',  i.e.,  in  the  number  and 
relative  proportions  of  toxins  present  in  the  cell.  Moreover,  the  organisms  can  be  studied 
in  batch  culture,  therefore  showing  changes  in  toxicity  as  the  cells  go  from  lag  over 
exponential  to  plateau  phase,  or  in  continuous  or  semi-continuous  culture,  where  toxicity 
can  be  determined  under  different  evironmental  conditions  during  balanced  and  adapted 
growth  of  the  dinoflagellates  (Anderson,  1990). 

A  significant  amount  of  research  has  dealt  with  changes  in  toxicity  and  hence  toxin 
production  rates  during  the  growth  cycle  of  dinoflagellates  in  batch  culture.  In  general, 
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toxin  content  per  cell  as  well  as  toxin  production  rates  are  highest  during  early  to  mid 
exponential  growth  of  the  population.  Coinciding  with  the  depletion  of  a  limiting  nutrient 
or  an  increase  in  environmental  stress  that  will  lead  the  population  into  the  non-growth 
stationary  phase,  toxin  content  and  synthesis  rate  decrease  in  parallel  to  a  general  decline  in 
cell  metabolism.  This  is  the  case  for  limiting  inorganic  nutrients  such  as  nitrogen  or  CO2 
and  environmental  factors  such  as  light,  salinity  and  temperature  (Proctor  et  al,  1975; 
White,  1978;  Hall,  1982;  Boyer  et  al,  1987;  Ogata  et  al,  1987;  Anderson  et  al,  1990b; 
Kim  et  al,  1993;  Flynn  et  al,  1994;  Usup  et  al,  1994;  Flynn  et  al,  1996;  Matsuda  et  al, 
1996).  A  widely  accepted  explanation  to  this  phenomenon  invokes  a  reallocation  of 
nitrogen  and  energy  sources  in  the  cell,  assuming  a  possible  role  of  STX  as  a  storage 
secondary  metabolite  for  nitrogen  and  a  reduced  availability  of  precursors  of  STX 
biosynthesis  such  as  arginine  and  ornithine,  which  would  result  in  lower  toxin  content 
(Boyer  et  al,  1987;  Anderson  et  al,  1990b).  One  exception  to  these  dynamics  is  growth 
under  phosphate-limited  conditions.  In  this  case,  toxin  content  and  toxin  production  rates 
increase  significantly  as  cells  enter  stationary  phase,  a  time  when  there  is  also  an  increase  in 
the  ratio  of  toxin  production  to  protein  or  free  amino  acid  synthesis  (Boyer  et  al,  1987; 
Anderson  et  al,  1990b;  Flynn  et  al,  1994;  Flynn  et  al,  1996).  This  decoupling  of  toxin 
production  from  other  cellular  processes  most  likely  mirrors  an  inhibitory  effect  of 
phosphate  deprivation  on  amino  acid,  and  most  importantly,  protein  metabolism  (Flynn  et 
al,  1994).  Such  a  process  could  induce  the  cell  to  allocate  some  of  the  nitrogen  not 
incorporated  into  or  catabolized  from  proteins  toward  highly  nitrogenous  compounds  such 
as  STX.  This  would  serve  as  a  means  of  nitrogen  storage  or  detoxification  from  potentially 
dangerous  levels  of  ammonia  due  to  non  incorporation  into  aminoacids  (Anderson,  1990; 
Flynn  et  al,  1996). 

In  regards  to  variations  in  toxin  composition  over  the  growth  cycle  of  dinoflagellates, 
the  general  observation  is  that  the  relative  amounts  of  the  different  STX  derivatives  remain 
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Fig.  3  Scheme  of  the  information  currently  available  on  transformations  among  STX 
derivatives  in  dinoflagellates.  Only  11-P-epimers  are  represented 
(Denominations  of  11-a-epimers  are  given  in  round  brackets).  Alternative 
names  of  some  compounds  are  indicated  in  square  brackets.  (Broad  arrows 
denote  confirmed  pathways;  dotted  arrow  represents  a  putative  transformation; 
the  compound  enclosed  in  square  brackets  has  not  been  isolated  from 
dinoflagellates). 
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constant  (Boyer  et  al,  1987;  Anderson  et  al,  1990b;  Flynn  et  al,  1994;  Sako  et  al,  1995; 
Matsuda  et  al.,  1996).  Changes  have  only  been  observed  under  extreme  growth-limited 
conditions  in  late  stationary  phase,  probably  due  to  a  general  loss  of  regulatory  capacity  or 
a  change  in  the  chemical  environment  of  the  cell  (Boczar  et  al,  1988),  and  under 
progressive  nitrogen  or  phosphorus  limitation  conditions  sustained  long  enough  for  the 
cells  to  adapt  their  metabolic  machinery  to  the  new  demands  (Anderson  et  al,  1990a; 
MacIntyre  etal,  1997). 

A  fundamental  difference  when  studying  toxin  production  dynamics  in  continuous  or 
semi  continuous  cultures  is  that  the  cells  are  analyzed  once  a  state  of  balanced  exponential 
growth  has  been  reached.  Therefore,  the  influence  of  factors  such  as  growth  rate  and 
residual  internal  pools  of  nutrients  are  minimized.  The  homeostasis  adopted  by  the  cell 
reflects  the  interaction  between  general  metabolism  and  toxin  biosynthesis.  The  effects  of 
nitrogen,  phosphate  and  temperature  limitation  on  STX  metabolism  have  been  studied 
(Anderson  et  al,  1990b;  Matsuda  et  al,  1996).  In  the  case  of  nutrient  limitation,  the 
growth  rates  of  the  organisms  are  correlated  to  the  availability  of  the  limiting  factor  in  the 
medium.  Nitrogen  limitation  induces  a  decrease  in  toxicity  (Anderson  et  al,  1990a; 
Matsuda  et  al,  1996).  Phosphate  limitation  has  the  opposite  effect  in  semi-continuous 
culture  (Anderson  et  al,  1990a),  which  is  also  observed  in  plateau  phase  in  batch  cultures 
(Boyer  et  al,  1987;  Anderson,  1990).  A  direct  correlation  has  also  been  observed  between 
growth  rates  and  toxin  production  rates  or  toxin  levels  (Anderson  et  al,  1990b;  Matsuda  et 
al,  1996). 

In  contrast  to  what  has  been  observed  in  batch  culture  and  during  progression  through 
the  growth  cycle,  toxin  composition  can  change  in  some  dinoflagellates  (Alexandrium 
fundyense  and  A.  tamarense)  depending  on  physico-chemical  parameters  (Anderson  et  al, 
1990a;  MacIntyre  et  al,  1997),  but  not  vary  significantly  in  other  species  (A.  catenella) 
under  similar  conditions  (Matsuda  et  al,  1996).  The  most  plausible  explanation  for  this 
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observation  is  that  the  shifts  in  toxin  profile  in  the  former  examples  were  a  product  of  long 
term  acclimation  of  the  cells  to  new  conditions  of  nutrient  deprived,  yet  balanced  growth 
(MacIntyre  et  al,  1997),  while  in  the  latter  experiments  with  A.  catenella  the  nutrient  stress 
in  the  cells  might  not  have  been  as  limiting  due  to  the  shorter  duration  of  the  experiments 
(Matsuda  et  al.,  1996). 

The  general  picture  that  emerges  from  these  studies  is  that  toxin  production  varies  in 
accordance  to  the  physico-chemical  environment  of  the  cell,  but  is  also  dependent  on 
cellular  processes  such  as  cell  cycle  and  protein  production.  Most  factors  affecting  toxin 
production  also  influence  the  general  metabolism  and  growth  of  the  cell,  e.g.,  temperature 
and  nitrogen  limitation.  On  the  other  hand,  conditions  such  as  phosphate  limitation  clearly 
do  not  affect  toxin  biosynthesis  in  the  same  way  they  modulate  other  cellular  processes, 
resulting  in  an  uncoupling  of  these  metabolic  activities  from  each  other  with  consequences 
for  the  toxicity  and  physiology  of  the  organism. 


GENETICS  OF  SAXITOXIN 

Analogous  to  other  processes  in  the  cell,  toxin  biosynthesis  is  ultimately  determined  by 
a  series  of  regulatory  events  that  modulate  the  activity  of  the  genes  encoding  the  toxin 
biosynthetic  enzymes.  The  identification  of  such  genes  remains  a  daunting  task,  as  the  use 
of  conventional  genetic  tools  such  as  mutational  analysis  or  transformation  with 
dinoflagellates  is  constrained  by  biological  and  technical  factors  including  lack  of 
information  on  the  number  of  copies  of  'toxin  genes',  limitations  in  culturing 
dinoflagellates  on  solid  media  and  the  high  content  of  DNA  per  cell  that  would  make  it 
necessary  to  screen  unrealistic  numbers  of  clones  in  order  to  detect  one  mutant. 
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Nonetheless,  classical  genetic  analysis  to  determine  the  location  of  the  toxin  genes  and 
some  first  attempts  to  use  molecular  techniques  are  underway. 

Alexandrium  is  a  heterothallic,  haploid  organism.  Under  the  appropriate  physico¬ 
chemical  conditions  vegetative  cells  will  produce  gametes.  Two  gametes  of  different 
mating  types  then  fuse  to  form  a  zygote  which  after  a  suitable  period  of  dormancy 
undergoes  meiosis,  resulting  in  tetrads  of  vegetative  cells  that  can  be  recovered  individually 
and  established  in  culture  (Pfiester  and  Anderson,  1987).  Genetic  analysis  of  haploid 
organisms  is  facilitated  by  the  direct  expression  of  the  nuclear  genotype  in  the  phenotype 
and  by  the  presence  of  tetrads,  which  permits  the  quantitative  analysis  of  meiotic  products 
and  the  determination  of  recombination  frequencies.  Heterothallism  is  instrumental  when 
trying  to  determine  if  a  genetic  trait  follows  uniparental  inheritance,  implying  mitochondrial 
or  chloroplastidial  localization  of  the  genes,  or  true  independent  Mendelian  segregation, 
which  would  suggest  nuclear  genes.  Finally,  quantitative  tetrad  analysis  allows  the 
presence  or  absence  of  extranuclear  or  cytoplasmic  inheritance  factors,  such  as  plasmids  or 
bacteria  to  be  assessed,  as  these  would  result  in  random  distributions  of  a  trait  in  the 
progeny. 

Crosses  of  toxic  and  non-toxic  strains  of  Gymnodinium  catenatum  result  in  a  1:1  ratio 
of  toxic  to  non-toxic  cells  in  the  FI  generation  and  a  toxin  profile  in  the  offspring  identical 
to  the  parental  (P)  one  (Oshima  et  al,  1993).  The  toxic  phenotype  further  segregates 
independently  from  mating  type  in  this  species.  Both  these  observations,  non-random 
inheritance  of  toxigenicity  and  independence  from  mating  type,  clearly  suggest  nuclear 
localization  of  the  'toxin  genes'.  Other  studies  in  which  toxic  strains  of  Alexandrium  sp.  or 
Gymnodinium  sp.  were  crossed  with  other  toxic  representatives  of  their  same  species  show 
a  1: 1  Mendelian  inheritance  pattern  for  toxin  profiles  that  is  also  independent  of  mating  type 
(Sako  et  al,  1992;  Ishida  et  al,  1993;  Oshima  et  al,  1993;  Kim  et  al,  1995;  Sako  et  al, 
1995).  In  some  cases,  over  a  50%  of  recombination,  relative  to  the  toxin  profiles,  has  been 
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observed  to  occur  in  the  FI  generation  in  Alexandrium  (Ishida  et  al,  1993).  In  test  crosses 
between  recombinants  the  phenotype  reverts  100%  to  the  parental  (P)  type.  Moreover, 
crosses  between  non-recombinant  FI  individuals  also  result  in  a  high  recombination 
frequency  (>50%)  and  F2  recombinant  phenotypes  identical  to  the  FI  recombinant 
phenotypes  (Ishida  et  al,  1993).  Back-crossing  of  a  recombinant  phenotype  with  one  of 
the  parental  phenotypes  results  in  a  non-recombination  cross.  This  strict  recombination 
pattern  suggests  a  close  link  among  the  STX  transforming  genes/enzymes.  Moreover,  the 
identical  recombinant  phenotypes  and  their  reversal  to  the  original  parental  phenotypes  also 
imply  one  recombination  site  and  the  transfer  of  only  one  specific  gene  or  set  of  genes  that 
is  responsible  for  the  phenotypic  variation.  Close  analysis  of  the  actual  compositional 
changes  resulting  from  recombination  are  unfortunately  too  complex  to  shed  any  more  light 
on  the  nature  of  the  transforming  enzymes. 

In  accordance  with  the  mating  studies  suggesting  that  no  extrachromosomal  factors  are 
involved  in  toxin  production,  efforts  to  physically  identify  plasmids  in  toxic  dinoflagellates 
and  cyanobacteria  have  also  failed  (Boczar  et  al,  1991;  Shimizu  et  al,  1996).  The 
combination  of  all  these  observations  unequivocally  suggests  that  the  genes  encoding  for 
the  STX  biosynthetic  and  interconverting  enzymes  in  dinoflagellates  are  localized  on 
nuclear  chromosomes. 


As  we  look  back,  we  see  that  numerous  studies  have  analyzed  the  chemical  structure 
and  activity  of  STX  and  its  derivatives  as  well  as  the  physiology  of  STX  production  by 
different  dinoflagellates  under  varying  physico-chemical  conditions.  Initial  steps  have  also 
been  taken  to  isolate  enzymes  involved  in  the  interconversions  among  STX  derivatives,  but 
the  fundamental  question  regarding  the  metabolic  role  of  STX  within  dinoflagellates 
remains  unanswered.  Only  a  complete  understanding  of  the  actual  mechanism  of  STX 
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biosynthesis  and  its  regulation  will  reveal  the  function  and  significance  of  STX  production 
in  dinoflagellates.  In  order  to  accomplish  this,  our  focus  has  to  be  directed  toward  new 
approaches  to  the  study  of  STX  biosynthesis  in  dinoflagellates. 

The  aim  of  this  thesis  was  to  study  STX  biosynthesis  in  the  dinoflagellate  Alexandrium 
fundyense  by  combining  high  resolution  physiological  studies  in  the  context  of  the  cell 
cycle  with  analysis  of  genes  expressed  during  toxin  production.  This  approach  was 
intended  at  identifying  genes  and  hence  enzymes  putatively  involved  in  STX  biosynthesis. 
In  the  first  chapter  the  close  link  between  toxin  production  and  the  cell  cycle  in 
synchronized  cultures  of  A.  fundyense  is  decribed.  In  the  second  chapter  the  influence  of 
temperature  and  phosphate  limitation,  two  factors  known  to  affect  STX  production,  on  the 
linkage  between  toxin  biosynthesis  and  the  different  phases  of  the  cell  cycle  is  analyzed. 
The  information  gathered  from  the  preceeding  chapters  is  used  in  chapter  three  in 
combination  with  a  novel  molecular  technique  called  differential  display  to  identify  genes 
expressed  while  toxin  biosynthesis  is  active  in  A.  fundyense. 
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CHAPTER  I 


Toxin  variability  during  the  cell  cycle  of  the  dinoflagellate 
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ABSTRACT 


Cultures  of  the  toxic  dinoflagellate  Alexandrium  fundyense  Balech  were  subjected  to 
conditions  which  induced  two  synchronized  divisions  over  a  period  of  48  hours.  Before, 
during  and  after  this  interval,  toxin  content,  toxin  composition,  and  several  other 
physiological  parameters  were  monitored  every  two  hours  for  94  hours.  Toxin  production 
was  discontinuous,  was  induced  by  light  and  always  occured  during  a  defined  time  frame 
within  the  Gi  phase  of  the  cell  cycle.  Specific  toxin  production  rates  were  positive  for  a 
period  of  approximately  eight  to  ten  hours  in  early  Gi  and  dropped  to  zero  for  the 
remainder  of  the  interphase  and  mitosis.  Analysis  of  toxin  composition  showed  that 
cellular  concentrations  of  all  the  saxitoxin  derivatives  followed  a  similar  pattern  of  increase, 
stabilization  and  decrease  throughout  one  generation  time.  A  putative  sequence  of 
interconversions  between  the  derivatives  could  be  established,  with  C2  as  the  first 
compound  to  appear.  Division  of  a  subset  of  the  population  during  the  first  24  hours  of  the 
experiment  and  the  ensuing  total  synchrony  of  the  culture  suggest  the  existence  of  two 
transition  points  in  the  cell  cycle  of  this  dinoflagellate.  The  first  transition  point,  at  the 
beginning  of  Gj,  is  light-dependent  and  holds  the  cells  in  a  Go-like  period.  The  second 
block  point  at  the  end  of  Gj  is  size-dependent  and  arrests  the  cells  in  Gj.  We  propose  a 
model  of  the  cell  cycle  of  A.  fundyense  in  which  progression  through  the  cell  cycle  can  be 
arrested  at  two  different  transition  points  located  in  Gj  and  toxin  production  is  induced  by 
light  during  Gj.  The  restriction  of  toxin  production  to  a  relatively  short  segment  of  the  cell 
cycle  provides  a  tool  for  comparing  cells  that  are  and  are  not  synthesizing  toxin.  This  has 
important  implications  and  applications  in  the  study  of  genetic  and  biochemical  regulation 
of  toxin  production  in  Alexandrium. 
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INTRODUCTION 


Dinoflagellate  species  in  the  genus  Alexandrium  are  responsible  for  paralytic  shellfish 
poisoning  (PSP)  along  the  northeastern  coast  of  the  United  States  and  Canada  (Maranda  et 
al.  1985;  Cembella  et  al.  1987;  Anderson  1989),  as  well  as  other  temperate  littoral  waters 
throughout  the  world  (Hallegraeff  1993).  PSP  is  caused  by  saxitoxin  (STX),  a  potent 
neurotoxin  produced  by  some  Alexandrium  species,  and  over  twenty  known  STX 
derivatives  which  differ  in  structure  and  toxicity  (Shimizu  1993).  An  important 
characteristic  of  the  STX-producing  dinoflagellates  is  that  the  toxicity  of  a  single  isolate  can 
vary  dramatically  under  different  growth  conditions.  This  variability  is  largely  due  to 
changes  in  toxin  content  (the  integrated  potency  or  total  molar  content  of  toxin  per  cell) 
during  various  stages  of  growth  (Hall  1982;  Boyer  et  al.  1987;  Anderson  et  al.  1990b),  but 
it  can  also  reflect  changes  in  toxin  composition  (the  relative  concentrations  of  the  different 
STX  derivatives)  (Boczeu-  et  al.  1988;  Anderson  et  al.  1990a). 

Most  studies  to  date  have  examined  the  toxin  content  of  dinoflagellates  on  time  scales  of 
days  to  weeks  using  batch  cultures.  These  studies  have  shown  that  toxin  content  varies 
with  nutrient  limitation  (Hall  1982;  Boyer  et  al.  1987;  Anderson  1990)  and  with  changes  in 
salinity  (White  1978),  temperature  (Hall  1982;  Ogata  et  al.  1987;  Anderson  et  al.  1990b)  or 
light  intensity  (Ogata  et  al.  1987).  A  typical  pattern,  observed  for  all  these  conditions 
except  phosphorous  limitation,  is  that  toxin  content  is  highest  in  early  exponential  growth 
and  decreases  as  the  cells  enter  stationary  phase  (Boyer  et  al.  1987;  Boczar  et  al.  1988; 
Anderson  et  al.  1990b).  Furthermore,  specific  toxin  production  rates  were  observed  to  be 
directly  proportional  to  specific  growth  rates,  suggesting  that  toxin  synthesis  might  be 
linked  to  cell  cycle  events  (Anderson  et  al.  1990b). 

Dinoflagellates  have  several  characteristics,  such  as  permanently  condensed 
chromosomes  and  the  absence  of  histones,  that  place  them  as  a  unique  group  among  the 
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eukaryotes  (Rizzo  1987;  Rizzo  1991).  However,  recent  evidence  also  suggests  that  their 
cell  cycle  is  controlled  by  universal  eukaryotic  regulatory  mechanisms  such  as  c<ic2-like 
protein  kinases  (Rodriguez  et  al.  1994;  Van  Dolah  et  al.  1995).  The  typical  eukaryotic  cell 
cycle  initiates  with  the  mitotic  division  of  the  cells  (M  phase)  followed  by  interphase, 
consisting  of  three  segments,  Gj,  S  and  G2.  During  the  first  segment,  Gj,  the  cells  are 
metabolically  very  active.  The  S  phase,  when  DNA  synthesis  occurs,  follows  Gj  and 
extends  until  chromosome  replication  occurs  and  the  DNA  content  of  the  cell  has  doubled. 
The  cell  then  enters  G2,  which  lasts  until  mitosis  occurs  again.  When  conditions  are 
unfavorable,  some  cells  can  withdraw  from  the  regular  cell  division  cycle  and  enter  a 
reversible  cell  cycle  arrest  phase  (Gg)  during  Gj  (Prescott  1976).  The  time  taken  to 
proceed  through  one  complete  cell  division  is  the  generation  time  of  the  cell. 

In  order  to  determine  the  onset  of  a  metabolic  process  within  the  cell  cycle,  as  was  our 
object,  all  the  cells  in  a  culture  must  proceed  through  the  different  stages  and  divide  at  the 
same  time.  Dinoflagellate  division  can  be  phased  by  the  photoperiod  and  occur  during  a 
clearly  defined  interval  such  as  the  dark-to-light  transition  (Chisholm  1981).  However, 
since  dinoflagellate  growth  rates  in  culture  are  typically  less  than  one  division  per  day,  only 
a  portion  of  the  population  will  divide  over  one  photoperiod.  Cells  that  are  not  ready  to 
divide  must  proceed  from  their  respective  positions  in  the  cell  cycle  through  DNA  synthesis 
and,  eventually,  mitosis  during  the  next  or  subsequent  light  period.  Measurements  on 
phased  cultures  thus  represent  an  average  of  the  physiological  status  of  cells  in  several 
different  metabolic  states  at  the  time  of  sampling  (Chisholm  1981).  Synchronized  cultures 
are  therefore  required,  yet  this  is  difficult  to  accomplish  with  dinoflagellates.  Several 
methods  have  been  used  in  the  past  to  synchronize  dinoflagellates.  Gonyaulax  polyedra 
was  synchronized  using  a  sieving  technique  that  separated  cell  cycle  stages  on  the  basis  of 
cell  size  (Homma  and  Hastings  1988).  This  approach,  however,  requires  excessively  large 
volumes  of  culture  in  order  to  obtain  sufficient  biomass  for  the  measurements  of  interest  to 
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our  study.  Amphidinium  carterae  was  synchronized  by  manipulating  the  length  of  the 
light-dark  cycle  (Galleron  and  Durrand  1976),  a  method  that  did  not  work  with 
Alexandriumfundyense  cultures  in  our  hands  (unpublished  data). 

Alternatively,  a  block/release  approach  can  be  used  to  synchronize  cells.  This  method 
is  based  on  the  hypothesis  that  cells  arrested  by  an  external  factor  at  a  specific  point  during 
their  cell  cycle  will  proceed  through  the  rest  of  the  cell  cycle  stages  as  a  homogeneous 
population  after  being  released  from  the  constraint  (Pardee  et  al.  1978).  Such  block  or 
"transition"  points  have  been  identified  during  late  Gj  in  mammalian  cells  (Pardee  et  al. 
1978)  and  yeast  (Johnston  et  al.  1980).  Our  efforts  to  locate  such  points  in  dinoflagellate 
cultures  using  metabolic  inhibitors  have  not  been  successful,  possibly  because  the  cell  wall 
of  dinoflagellates  is  impermeable  to  many  organic  molecules  or  because  the  autotrophic 
nature  of  the  organisms  does  not  facilitate  uptake  of  the  inhibitors  (unpublished  data). 

One  factor  that  can  arrest  cells  at  a  particular  point  of  their  cell  cycle  is  light.  Light  is  of 
primary  importance  for  activating  physiological  processes  as  well  as  determining  the  timing 
of  cell  cycle  events  in  phytoplankton  (Chisholm  1981;  Vaulot  1985).  In  phytoplankton, 
light-dependent  block  points  can  be  passed  by  a  cell  only  after  it  has  been  exposed  to 
sufficient  light.  Such  transition  points  have  been  documented  in  Chlamydomonas  (Spudich 
and  Sager  1980)  Amphidinium  carterii  (Olson  and  Chisholm  1986)  and  Hymenomonas 
carterii  and  Thalassiosira  weisflogii  (Vaulot  et  al.  1986).  If  the  cells  are  deprived  of  the 
minimal  light  requirement  and  there  is  only  one  light-dependent  transition  point,  and  if  the 
length  of  dark  arrest  is  equal  to  or  greater  than  the  entire  cell  cycle  (i.e.  one  generation 
time),  all  cells  should  be  blocked  at  the  same  point.  Once  the  normal  light-dark  cycle  is 
restored,  a  synchronized  population  is  obtained. 

In  this  study,  we  employed  a  block/release  approach,  using  prolonged  darkness  to 
arrest  A.fundyense  in  a  light-dependent  section  of  the  cell  cycle.  Cells  were  released  to  a 
normal  14:10  hr  light:dark  cycle  and  closely  spaced  measurements  made  to  determine  the 
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small-scale  dynamics  of  toxin  production.  The  parameters  measured  included  toxin 
concentration,  toxin  composition,  DNA  content,  cell  size  and  total  protein. 


MATERIALS  AND  METHODS 

The  Organism.  The  dxnoilsigQihtQ  Alexandrium  fundyense  (strain  GtCA29)  was  used 
throughout  this  study.  This  clonal  culture  was  established  in  January  1985  from  a  cyst 
isolated  from  Gulf  of  Maine  sediments,  20  miles  east  of  Portsmouth,  New  Hampshire 
(USA).  Prior  to  the  experiments,  cultures  were  maintained  at  15°C  in  f/2  medium 
(Guillard  and  Ryther  1962)  modified  by  addition  of  H2Se03  (lO’^M)  and  by  reducing  the 
concentration  of  CUSO4  to  lO'^M.  Vineyard  Sound  seawater  (0.2  pm  filtered,  31°/oo 
salinity)  was  used  as  the  medium  base.  Throughout  the  experiment,  irradiance  of  ca.  250 
pEm'^s*'  was  provided  by  cool  white  fluorescent  bulbs  on  a  14:10  h  lightrdark  cycle. 

Synchronization  Experiment.  Three  20L  borosilicate  glass  carboys  containing  17L  of 
sterile  filtered  f/2  medium  were  inoculated  with  acclimated  exponential  phase  cells  to  an 
initial  cell  density  of  ca.  200  cells  ml'^  The  cultures  were  incubated  at  20°C  and  aerated  by 
gentle  bubbling  with  activated-charcoal-filtered  air.  After  completion  of  the  light  cycle  on 
the  third  day  after  inoculation  (ca.  400  cells  ml'^),  the  incubator  lights  were  turned  off  for  a 
period  of  82  hours.  One  hour  before  the  normal  14:10  LD  cycle  was  restored  (to  of  the 
experiment),  samples  for  toxin,  DNA,  cell  concentration,  and  protein  were  taken  from  each 
of  the  carboys.  Samples  were  collected  from  all  three  carboys  at  hrs  2,  4,  12  and  18,  and 
every  two  hours  from  hrs  22  through  94.  The  samples  were  taken  from  each  carboy  by 
blocking  the  aeration  exhaust  port  of  the  carboy,  thereby  pressurizing  the  vessels  and 
forcing  a  sample  of  the  culture  through  a  sample  tube  into  the  sampling  containers.  The 
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cells  were  not  mixed  manually  or  disturbed  by  any  other  physical  means.  During  the  dark 
period  and  in  order  to  ensure  that  the  cultures  were  not  exposed  to  light,  sampling  was 
performed  under  a  6  W,  incandescent  red  light,  and  special  precautions  were  taken  to 
minimize  exposure  to  stray  light. 

Toxin  Analysis.  Between  2.5  and  5  x  10^  cells  were  collected  from  the  sub-samples  by 
centrifugation  (3,000  x  g,  23°C).  The  resulting  pellet  was  rinsed  into  a  centrifuge  tube  and 
resuspended  in  0.5  M  acetic  acid  for  extraction  by  sonification.  The  extracts  were  stored  at 
-20”C  prior  to  analysis.  The  extracts  were  analyzed  by  HPLC  (Oshima  et  al.  1989) 
incorporating  the  modifications  previously  described  by  Anderson  et  al.  (1994). 

Soluble  Protein.  At  least  50,000  cells  were  collected  by  centrifugation  (3000  x  g,  23°C) 
for  protein  quantification.  After  aspiration  of  the  supernatant,  the  pellet  was  stored  frozen 
at  -20®C.  For  analysis,  1  ml  of  0.1  N  NaOH  was  added  to  the  pellet  and  the  suspension 
sonified  for  15  seconds  to  disrupt  the  cells.  Samples  were  then  digested  for  20  min  at 
80”C,  centrifuged,  and  the  supernatant  removed  and  neutralized  with  0.1  volumes  of  1.0  N 
HCl  (Binder  and  Anderson  1990).  A  50  pi  fraction  of  the  neutralized  solution  was  mixed 
with  1  ml  of  bicinchoninic  acid  (BCA)  working  solution  (Pierce,  Rockford  IL),  and  the 
absorbance  of  duplicate  samples  determined  at  562  nm. 

Flow  Cytometric  Analysis.  A  15  ml  aliquot  of  culture  was  preserved  with  5%  formalin 
(v/v)  and  stored  at  4”C.  The  cells  were  centrifuged  (1700  x  g,  23°C,  3  mins)  and  the  pellet 
resuspended  in  10  ml  of  ice  cold  methanol  and  stored  at  4“C  to  extract  chlorophyll.  The 
samples  were  then  further  concentrated  by  centrifugation.  A  subsample  of  the  methanol- 
preserved  material  was  transferred  into  a  1.5  ml  microcentrifuge  tube  and  pelleted  at  2000 
rpm.  The  supernatant  was  removed  and  the  cells  washed  with  1  ml  PBS  (40mM 
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Na2HP04,  22mM  KH2PO4,  85mM  NaCl).  After  centrifugation,  the  pellet  was 
resuspended  in  0.5  ml  propidium  iodide  (PI)  staining  solution  (4  pg  ml'^  PI,  1250  units 
RNaseA  ml'*)  and  allowed  to  react  in  the  dark  for  at  least  two  hours  prior  to  analysis  on  an 
Epics  V  flow  cytometer  (Coulter  Electronics,  Hialeah,  FL)  (Olson  et  al.  1986).  The  results 
were  reduced  by  running  an  inside-out  and  an  outside-in  model  with  trapezoidal  S  phase 
fitting  on  ModFit  5.11  (VERITY  Software  House,  Inc.,  Topsham,  ME). 

Cell  Size.  Cell  diameter  was  determined  by  analyzing  unfixed  cells  (3  ml  of  the  culture) 
on  a  Multisizer  apparatus  (apperture  100pm)  (Coulter  Counter  Ltd.,  GB)  within  30  min 
after  sampling.  Mean  sizes  were  calculated  from  the  Gaussian  distributions  obtained. 

Cell  Concentrations.  Cell  concentrations  were  determined  by  counting  Utermohls- 
preserved  samples  (Utermohl  1958).  At  least  400  cells  per  sample  were  counted  in  a 
Sedgwick  Rafter  counting  chamber.  Cell  densities  represent  the  mean  of  the  averaged 
triplicate  counts  for  each  carboy. 

Data  fitting.  Simple,  segmented  and  smoothed  linear  models  were  fitted  to  the  cell 
growth  and  toxin  concentration  data.  The  fitting  program  uses  the  Marquardt-Levenberg 
algorithm  to  minimize  the  sum  of  the  squares  of  the  differences  between  the  equation  values 
and  the  data  (SigmaPlot,  Jandel  Corp.).  This  is  an  iterative  process  that  lasts  until  the 
model  and  the  data  converge.  Convergence  is  achieved  when  the  absolute  value  of  the 
difference  between  the  square  root  of  the  sum  of  squares  of  the  residuals  is  less  than 
0.0001. 

The  fitting  was  done  using  the  equation: 


38 


/n(r)  =  kt  +  K 


(1) 


where  f(t)  is  a  linear  function  of  time,  k  is  the  growth  rate  (or  toxin  production  rate;  units: 
division  hour"*  and  pmole  ml'*  hour"*  respectively)  and  Nn  is  the  cell  concentration  (or 
toxin  concentration;  units:  cells  ml'*  and  pmoles  ml'*  respectively)  at  the  end  of 
compartment  n-1.  If  X  is  a  parameter  for  cell  concentration  (or  toxin  concentration),  n  is 
the  number  of  the  curve  segment,  and  T  is  a  parameter  for  time  at  the  endpoints  of  the 
different  segments,  the  rate  k  between  Tn  and  Tn-i  will  be  given  by 


Xn-X 


n  -  1 


Tn  - Tn -  1 


(2) 


Since  /(Tn)  =  Xn,  Nn  can  be  defined  as 


Nn  = 


Xn 


1 


Tn  -  XnTn  -  1 


Tn  - Tn  -  I 


(3) 


Incorporating  (2)  and  (3)  into  (1)  yields  a  function  /„(?) 


Xn(Tn  ■?)■!■  Xn  +  l(f  “  Tn  -  l) 
Tn  - Tn -  1 


(4) 


where  f(t)  is  the  linear  function  of  time  that  describes  the  changes  in  cell  number  (or  toxin 
concentration)  in  each  of  the  segments  on  Figs.  lA  and  IB.  Initially,  the  different  values 
for  T  were  arbitrarilly  specified  at  the  light-dark  and  dark-light  transitions.  After  several 
iterations  of  the  model,  the  time  points  T  defining  the  beginning  and  ending  of  each 
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segment  were  determined  from  the  best  least-square  fit  to  the  entire  data  set.  The  results 
were  then  plotted  as  a  segmented  curve  together  with  the  actual  data. 

Abbreviations  used  throughout  this  text  are:  SIX  =  saxitoxin;  NEO  =  neosaxitoxin; 
GTX  1,2, 3, 4, 5  =  gonyautoxins  I,  II,  III,  IV  and  V;  C2  =  toxin  C2;  Go,  Gi  and  G2  =  Gap 
0,  Gap  1  and  Gap  2  phases  of  the  cell  cycle;  S  =  synthesis  phase;  M  =  mitosis. 


RESULTS 

All  results  represent  the  mean  values  obtained  by  averaging  data  from  three  independent 
cultures  grown  in  parallel  and  started  with  the  same  acclimated  inoculum. 

Synchrony  and  cell  growth.  When  the  dark-arrested  cultures  were  returned  to  a  10: 14 
hr  light-dark  cycle,  cell,  numbers  remained  constant  for  18  hours,  with  an  average 
concentration  (+  SD)  of  723  (±  61)  cells  ml'*  (Fig.  1  A).  By  hr  26  a  new  stable  population 
size  of  1012  (±51)  cells  ml'^  was  observed,  a  40%  increase  from  initial  levels.  This  was 
followed  16  hrs  later  by  a  major  surge  in  division,  with  concentrations  doubling  to  1999  (± 
135)  cells  mT^  over  an  8  hr  interval.  This  cell  density  was  sustained  until  hr  66  when  a 
new  synchronized  division  event  began  that  brought  the  population  to  3465  (+  142)  cells 
mb  ^  a  1.8  fold  increase.  Cell  proliferation  always  started  3  to  6  hours  into  the  dark  phase 
and  lasted  for  approximately  9  hours. 

Linear  regressions  were  calculated  for  each  of  the  intervals  described  above.  Only  the 
sections  delimited  by  hours  18-26,  42-50  and  66-74  had  slopes  that  were  significantly 
different  from  zero  (p<0.001).  The  horizontal,  'non-growth'  segments  had  slopes  which 
were  not  statistically  different  from  zero  (p>0.93,  C.I.  97%).  These  stepwise  linear 
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increases  in  cell  concentrations  were  modeled  by  fitting  a  simple,  segmented  and  smoothed 
linear  model  to  the  data  (see  materials  and  methods).  The  fitted  curve  matched  the  data 
extremely  well,  as  shown  in  Fig.  1  A. 

Cell  Size  Distributions.  Cell  diameter  was  determined  for  hrs  0  through  56,  and  three 
distinct  size  classes  could  be  distinguished  (Fig.  2).  At  the  end  of  the  dark  arrest  period,  hr 
0,  the  cells  had  an  equivalent  spherical  diameter  (±  SD)  of  23  pm  (±  1.2).  This  peak 
corresponds  to  the  size  of  cells  in  Gq  (Fig.  4).  During  the  first  24  hours  these  cells 
increased  in  size  to  29  pm  (±  1.4),  coinciding  with  a  peak  in  Gi.  Over  the  next  8  hours 
their  diameter  remained  constant,  starting  to  increase  again  at  hr  34  and  reaching  a 
maximum  of  33  pm  (+  1.6),  at  hour  44,  corresponding  to  the  peak  of  the  G2  +  M  cell  cycle 
stage.  Subsequently,  the  cells  rapidly  decreased  in  size  to  29  pm  (±  1.3). 

Protein  Content.  Protein  content  per  cell  was  400  pg  (±  65)  in  the  dark  arrested  cells, 
but  increased  once  the  light  was  restored,  exhibiting  a  typical  bi-modal  pattern  through  time 
that  peaked  at  the  end  of  the  light  periods  (hrs  15,  39,  63  and  87),  and  reached  minima  at 
the  end  of  the  dark  periods  (hrs  25,  49  and  73)  (Fig.  3).  The  minimum  was  constant  at 
523  pg  (±  58).  The  first  maximum,  at  hr  15,  was  at  only  608  pg  (+  79).  All  other  maxima 
were  at  1 106  pg  (±  126). 

Cell  Cycle  Stages.  Two  cohorts  of  cells  were  apparent  in  the  early  stages  of  the 
experiment  (Fig.  4).  The  first,  comprising  approximately  30%  of  the  total  population, 
appeared  1 1  hours  into  the  first  light  period,  evidenced  as  an  increase  in  the  percentage  of 
cells  in  S  and  a  concurrent  decrease  in  Gi  cells.  Shortly  thereafter,  a  small  peak  in  G2  +  M 
occured,  after  which  96%  of  the  population  was  in  Gj  phase  at  hr  22.  This  corresponded 
to  a  30%  increase  in  total  cell  concentration  at  that  time.  The  second  cohort  of  cells  (the 
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Fig.l  Cell  division  pattern,  changes  in  toxin  concentration  of  the  population  and 
specific  toxin  production  rates,  Ptox>  Alexandrium  fundyense  Balech 
(GtCA29)  as  a  function  of  time  and  during  synchronous  growth  on  a  light;  dark 
cycle  of  14:10  hours  (dark  periods  indicated  by  shaded  areas).  (A)  Cell 
densities  and  (B)  total  toxin  concentration  in  the  culture  measured  at  2-hour 
intervals  (Error  bars  denote  +  SD).  Data  were  fitted  to  simple  segmented  linear 
models  (continuous  line;  see  text).  (C)  Specific  toxin  production  rates, 

(h‘^),  were  obtained  from  the  fitted  curve  in  (B)  (see  text). 
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Fig.  2  Variation  in  the  cell  size  distribution  of  Alexandrium  fundyense  Balech 
(GtCA29)  over  the  first  48  hours  of  the  synchronization  experiment  (height  of 
areas  is  equivalent  to  number  of  cells).  Area  with  vertical  lines  corresponds  to 
cells  in  Go  (23  ±1.2  pm),  black  areas  to  Gi  (29  ±1.4  pm),  and  the  horizontally 
lined  area  to  G2  (33  ±1.6  pm). 
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Fig.  3  Protein  production  in  Alexandriumfundyense  Balech  (GtCA29)  as  a  function  of 

time  and  during  synchronous  growth  on  a  light:dark  cycle  of  14:10  hours  (dark 
periods  indicated  by  shaded  areas).  Protein  content  was  measured  at  2-hour 
intervals  (Error  bars  denote  ±  SD). 
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Fig.  4  Relative  distributions  of  Alexandriumfundyense  Balech  (GtCA29)  cells  in  the 
different  cell  cycle  stages  as  a  function  of  time  and  during  synchronous  growth 
on  a  light:dark  cycle  of  14:10  hours  (dark  periods  indicated  by  shaded  areas). 
The  different  cell  cycle  stage  fractions  represent  the  percentage  of  total  cells  in 
each  phase  and  were  determined  by  analysing  the  varying  amounts  of  DNA  per 
cell  (CV<5%).  Error  bars  not  shown  for  clarity. 
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remainder  of  the  cells)  did  not  divide  during  the  first  36  hours  of  the  experiment.  This 
larger  cohort  was  arrested  in  what  we  can  define  as  Go  based  on  DNA  content,  protein 
content,  toxin  content  and  cell  size  (see  discussion).  From  the  beginning  of  the  second  L:D 
cycle  on  (hr  26),  the  two  populations  progressed  together  through  the  rest  of  the  cell 
cycles,  as  corroborated  by  the  cell  cycle  patterns  (Fig.  4),  the  cell  size  distributions  (Fig.  2) 
and  the  regular  increases  in  cell  density  (Fig.  lA). 

The  duration  of  the  cell  cycle  phases  was  determined  from  the  variations  in  DNA 
content  over  the  duration  of  the  experiment.  Assuming  that  the  cells  that  go  first  through 
one  cell  cycle  stage  are  the  ones  also  going  first  through  the  subsequent  stages,  and  after 
determining  the  transition  times  between  the  peaks  corresponding  to  the  different  phases 
(Beck  1978;  Yamaguchi  and  Imai  1994),  Gi  was  calculated  to  be  approximately  16  hrs,  the 
S  phase  2  hrs,  and  G2  +  M,  6  hrs.  Cells  in  G2  and  cells  undergoing  mitosis  (M)  have  to  be 
regarded  as  one  group  because  their  cytological  differences  cannot  be  resolved  with  flow 
cytometry. 

Toxin  Content.  When  the  culture  was  removed  from  the  dark  and  the  normal  light-dark 
cycle  re-established,  the  total  toxin  in  the  culture  (toxin  ml’^  (±  SD))  increased  from  32  (+ 
8)  pmol  ml'^  to  62.4  (±  2.2)  pmol  ml'^  over  a  period  of  eight  hours  and  then  remained 
stable  for  16  hours.  Following  this  a  new  surge  in  toxin  production  was  observed  at  the 
beginning  of  the  light  period  at  hr  25.  This  increase  continued  for  8  hours,  reaching  a  new 
constant  level  of  102.6  (±  7.7)  pmol  ml'^  The  same  pattern  of  toxin  increase  and  plateau 
was  observed  during  the  subsequent  L:D  cycles,  with  stable  concentrations  of  189  (±  17.5) 
and  347.4  (±  40.8)  pmol  ml'^  respectively.  The  mean  increase  in  concentration  during 
each  toxin  production  interval  was  90%.  This  pattern  of  toxin  production  coincides  with 
the  dynamics  of  toxin  accumulation  per  cell  as  exemplified  in  Fig.  5  for  hours  20  through 
54.  Toxin  content  per  cell,  as  calculated  from  toxin  concentrations  in  the  culture  and  the 
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corresponding  cell  densities,  followed  a  regular  pattern  of  increase,  stabilization  and 
decrease  throughout  each  generation  time. 

Linear  regressions  were  calculated  for  each  of  the  segments  described  in  the  toxin  ml'^ 
data.  Only  the  sections  delimited  by  hours  1-12,  26-34,  48-60  and  72-82  had  slopes  that 
were  significantly  different  from  zero  (p<0.001).  The  horizontal,  'non-production' 
segments  had  slopes  which  were  not  statistically  different  from  zero  (0.9>p>0.57,  C.I. 
89%).  These  stepwise  linear  increases  in  toxin  concentrations  were  modeled  by  fitting  a 
simple,  segmented  and  smoothed  linear  model  to  the  data  (see  materials  and  methods).  In 
this  case,  X  was  the  parameter  for  toxin  concentration  (Fig.  IB).  The  fitted  curve  is  plotted 
with  the  data  in  Fig.  IB. 


Specific  Toxin  Production  Rates.  The  specific  toxin  production  rate,  pxox  >  was 
defined  as  the  change  in  toxin  concentration  relative  to  the  existing  concentration  in  the 
synchronized  cultures.  This  definition  is  based  on  the  previously  described  specific  toxin 
production  rate  in  asynchronous  cultures,  pxox  (day'^)  (Anderson  et  al.  1990b).  The  rates 
were  obtained  using  the  equation: 


t\  -  to 


(5) 


where  jUxox  is  the  specific  toxin  production  rate  (h‘^)  and  Ti  and  To  are  the  toxin 
concentrations  at  times  ti  and  to  respectively,  as  obtained  from  the  fitted  curve  in  Fig.  IB. 
A  regular  pattern  of  toxin  production,  in  which  toxin  production  rates  are  positive  for  a 
period  of  approximately  8  hours  and  then  drop  to  zero  for  the  remaining  16  hours  of  the 
cell  cycle,  is  revealed  (Fig.  1C). 
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Toxin  composition.  The  relative  proportions  (mole  percent  of  total  toxin)  of  the  five 
main  saxitoxins  remained  essentially  unchanged  throughout  the  experiment  (data  not 
shown).  The  accumulation  patterns  of  the  seven  main  derivatives  of  saxitoxin,  C2,  NEO, 
STX,  GTX2,3  isomers  and  GTX1,4  isomers,  were  very  similar  (Fig.  5).  Analysis  of  the 
exact  onset  of  accumulation  for  each  of  the  toxins  suggests  that  C2  increases  first,  followed 
by  STX  and  the  isomeric  pair  GTX  2,3.  The  concentrations  of  NEO  and  the  pair  GTX1,4 
begin  to  increase  approximately  1-2  hours  later. 


DISCUSSION 

This  study  documents  for  the  first  time  the  extent  to  which  the  toxicity  of  the 
dinoflagellate  Alexandrium  fundyense  varies  on  short  (hourly)  time  scales.  The 
fluctuations  in  toxin  production  were  evaluated  in  the  context  of  the  cell  cycle,  a  framework 
that  allowed  us  to  divide  one  cell  generation  into  a  series  of  shorter  stages  that  could  easily 
be  delineated  using  DNA  measurements.  Synchronization  of  the  cells  ensured  that  virtually 
the  entire  population  was  in  the  same  physiological  state  at  each  sampling  point  and 
therefore  that  the  results  of  the  analysis  of  the  population  as  a  whole  could  be  extrapolated 
to  the  individual  cell. 

Synchrony  and  growth  patterns 

Immediately  following  release  from  the  dark,  Alexandrium  shows  a  lag-phase 
equivalent  to  one  generation  time.  Such  non-proliferation  periods,  previously  described  for 
induced  synchrony  in  yeast  and  other  eukaryotes  (Mitchison  1969;  Pardee  et  al.  1978),  are 
the  result  of  a  Go-like  segment  in  the  cell  cycle  that  prevents  the  cells  from  progressing 
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through  Gi  and  dividing.  The  cells  in  this  stage  are  of  smaller  size  (Fig.  2).  After  a  period 
of  24  hours  the  cells  regain  the  normal  size  of  about  30  pm  in  diameter  for  this  species 
(Fig.  2)  with  96%  of  the  cells  being  in  Gi  (Fig.  4).  Over  the  following  24  hours  cell  size 
increases  to  reach  a  maximum  during  the  peak  of  G2  (Figs.  2  and  4).  Thereafter  the  cells 
decrease  rapidly  in  size  to  values  equivalent  to  Gi  cells.  The  smaller  size  correlated  to  cells 
in  Go  is  never  seen  again  during  the  experiment.  These  data  suggest  the  presence  of  a  size- 
dependent  transition  point  at  the  end  of  Gi,  similar  to  other  eukaryotes  (Mitchison  et  al. 
1991).  The  degree  of  synchrony  achieved  in  the  two  subsequent  divisions,  as  determined 
from  the  relative  proportions  of  cells  in  the  different  cell  cycle  stages  and  the  increase  in  cell 
numbers  (Figs.  3  and  lA),  was  100%  and  90%  respectively.  Two  further  divisions  (data 
not  shown)  show  increments  of  89%  and  80%.  These  numbers  imply  that  even  though 
synchrony  was  achieved  during  the  first  proliferation  step,  it  was  progressively  lost. 
Similar  patterns  have  been  previously  described  for  other  eukaryotes  such  as  yeast 
(Mitchison  et  al.  1991).  The  most  widely  accepted  explanation  for  this  loss  of  synchrony 
is  the  fact  that  the  cells  are  growing  at  their  maximum  possible  growth  rate,  a  situation  that 
is  not  sustainable  as  the  cells  are  not  in  balanced  growth,  when  uptake  equals  metabolic 
needs  (Mitchison  1969).  A  similar  situation  is  true  for  Alexandrium  fundyense.  This 
species  usually  grows  at  a  rate  of  0.3  to  0.5  divisions  per  day  in  culture,  yet  under  the 
specific  conditions  of  our  experiment,  the  cells  were  forced  into  dividing  once  per  day,  a 
pace  that  is  not  sustainable  over  long  periods  of  time.  Such  "explosive  growth",  events 
during  which  populations  of  organisms  can  grow  very  fast,  are  believed  to  be  a  possible 
mechanism  for  bloom  formation  (Smayda,  pers.  comm.). 

Temporal  Pattern  in  Toxin  Production 

In  analyzing  the  variability  in  toxin  production  during  this  study,  we  chose  to  focus  on 
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Fig.  5  Changes  in  total  toxin  and  toxin  composition  per  cell  of  Alexandriumfundyense 
Balech  (GtCA29)  as  a  function  of  time  and  during  synchronous  growth  on  a 
light:dark  cycle  of  14:10  hours  (dark  periods  indicated  by  shaded  areas).  The 
Y-axis  is  discontinuous  between  35  and  55  fmol/cell.  Error  bars  are  not  plotted 
for  the  single  toxin  derivatives  for  sake  of  clarity. 
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two  parameters.  First,  the  variations  of  toxin  per  unit  volume  of  culture  (toxin  ml'^)  were 
determined,  for  this  parameter  is  an  indicator  of  net  toxin  production  and  should  never 
decrease  unless  toxin  is  catabolized  or  leaked  from  the  cell.  Secondly,  toxin  cell'^  was 
derived  from  the  previous  data  to  follow  the  dynamics  of  toxin  accumulation  on  a  per  cell 
basis  (Note  that  we  use  the  term  "toxin  production"  to  mean  toxin  accumulated  in  the  cell, 
ignoring  losses  due  to  leakage  or  catabolism).  The  most  important  conclusion  from  the 
data  presented  here  is  that  toxin  production  is  induced  by  light  and  occurs  in  a  defined 
interval  within  the  Gi  phase  of  the  cell  cycle.  Specific  toxin  production  rates  were  positive 
for  a  period  of  eight  to  ten  hours  in  early  Gi  and  remained  zero  for  the  remainder  of  the 
interphase  and  mitosis  (Fig.  1C),  six  hours  of  which  were  still  in  light.  This  behaviour 
was  not  only  observed  once  the  population  was  synchronized,  but  also  during  the  first  24 
hours,  when  the  cells  were  not  yet  progressing  through  the  "regular"  cell  cycle. 

If  toxin  were  produced  continuously  and  at  a  constant  rate,  toxin  mf*  should  increase 
exponentially  through  time.  In  Fig.  IB,  this  would  plot  as  a  straight  line  with  constant 
slope.  Instead,  the  data  reveal  a  step-like  pattern  similar  in  form  to  the  steps  seen  in  cell 
number  when  counted  at  closely  spaced  intervals  through  time  (Fig.  lA).  This  is  further 
confirmed  when  toxin  accumulation  per  cell  is  analyzed,  as  exemplified  in  Fig.  5.  Toxin 
accumulates  (is  produced)  over  a  period  of  8-10  hours  before  its  concentration  reaches  a 
constant  value.  This  level  is  maintained  until  the  cells  start  dividing  (Fig.  1  A)  and  the  toxin 
is  most  probably  partitioned  between  the  two  daughter  cells  (Fig.  5),  as  toxin  concentration 
in  the  culture  does  not  change  during  mitosis  (Fig.  IB).  The  average  increase  in  toxin 
concentration  over  the  duration  of  the  experiment  is  exponential  but  is  determined  by  an 
hourly  scale,  step-like  pattern  of  toxin  production.  Such  discontinuous  production  or 
activity  patterns  in  synchronized  cultures  have  been  previously  described  in  other 
organisms  and  are  universal  for  cellular  processes  that  are  activated  by  external  or  internal 
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factors  rather  than  being  continuous  throughout  the  cell  cycle  (Mitchison  1969). 

Given  this  discontinuous  pattern  of  toxin  production  and  increase  in  cell  concentration, 
we  decided  to  fit  a  segmented  linear  model  to  each  of  the  data  sets.  The  underlying 
assumption  was  that  the  production  rate,  similarly  to  the  growth  rate,  follows  a  Gaussian 
distribution  during  the  "accumulation"  or  "doubling"  periods,  resulting  in  a  sigmoidal  curve 
for  the  increase  in  toxin  concentration  or  cell  number  (Mitchison  1969;  Vaulot  1985).  For 
simplicity  purposes,  such  sigmoidal  curves  can  be  approximated  by  a  straight  line.  The 
slope  of  this  straight  line  would  represent  the  mean  production  rate  or  growth  rate  during 
that  period.  Initially,  the  time  points  (parameter  T,  see  materials  and  methods)  delimiting 
each  segment  of  the  cell  or  the  toxin  concentration  curve  on  Figs.  1 A  and  IB  respectively, 
were  arbitrarilly  defined  at  the  times  when  the  lights  went  on  or  off.  After  several 
iterations,  the  time  intervals  defining  the  segments  were  adjusted  through  a  least-square, 
best-fit  to  the  data.  The  duration  and  timing  of  each  segment  thus  determined,  closely 
matched  the  segments  we  had  originally  determined  for  calculating  the  slopes  (Figs.  1 A  and 
IB)  and  in  the  case  of  toxin  production,  the  model  also  coincided  with  the  patterns 
observed  in  toxin  accumulation  per  cell  (Fig.  5).  These  results  confirmed  the  observation 
that  toxin  production  occurred  during  an  eight  to  twelve  hour  interval  in  early  Gi. 

Specific  toxin  production  rates  were  obtained  using  the  values  from  the  fitted  curve  in 
Fig.  IB.  A  regular  pattern  of  toxin  production  is  revealed  in  which  toxin  production  rates 
are  positive  for  a  period  of  approximately  8-10  hours  and  then  drop  to  zero  for  the 
remaining  16  hours  of  the  cell  cycle  (Fig.  1C).  In  the  three  consecutive  cycles  we 
analyzed,  toxin  production  was  induced  by  light  at  the  beginning  of  the  light  period  and 
stopped  approximately  eight  to  ten  hours  into  Gi. 

Protein  production  did  not  mimic  toxin  production  throughout  time  (Fig.  3). 
Accumulation  of  soluble  protein  was  also  induced  by  light,  but  it  continued  during  the  light 
period  and  was  reversed  with  the  beginning  of  the  dark  period.  Identical  protein 
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accumulation  dynamics  have  been  previously  observed  in  other  phytoplankton  (Howell  et 
al.  1977;  Hastings  et  al.  1981;  Puiseux-Dao  1981). 

We  conclude  from  these  data  that,  unlike  protein  production,  which  extends  over  the 
entire  interphase  comprising  Gj,  S  and  G2  and  is  solely  dependent  on  light,  toxin 
production  is  induced  by  light,  is  confined  to  an  early  portion  of  the  Gi  phase  only,  and 
ceases  before  the  light  phase  ends. 

Light  Induction 

It  has  been  postulated  in  a  previous  study  that  toxin  production  can  be  influenced  by 
light.  Alexandrium  {=Protogonyaulax)  tamarensis  was  shown  to  increase  its  toxicity  when 
the  growth  rate  was  reduced  by  a  decrease  in  light  intensity  (Ogata  et  al.  1987).  This  result 
is  analogous  to  those  obtained  when  growth  rates  are  reduced  as  a  result  of  nutrient 
limitation  (Proctor  et  al.  1975;  Hall  1982;  Boyer  et  al.  1985;  Anderson  et  al.  1990b)  or 
temperature  (Hall  1982;  Ogata  et  al.  1987;  Anderson  et  al.  1990b).  Furthermore,  a  direct 
correlation  between  specific  growth  rates  and  specific  toxin  production  rates  during  the 
exponential  phase  in  semi-continuous  cultures  has  been  observed,  suggesting  a  possible 
relationship  between  toxin  production  and  the  cell  cycle  (Anderson  et  al.  1990b).  These 
observations  can  be  explained  taking  into  consideration  our  results  that  show  that  toxin 
production  is  clearly  induced  by  light,  lasts  only  for  8-10  hours  and  is  confined  to  the  early 
Gi  phase  (Figs.  1C  and  4).  Cells  of  the  dinoflagellate  Amp/uWmmm  carteri  growing  at  a 
slow  rate  do  so  by  increasing  the  duration  of  the  Gi  phase  while  keeping  the  length  of  the 
other  phases  relatively  constant  (Olson  and  Chisholm  1986).  If  a  cell  is  arrested  in  Gi  for 
several  light;  dark  cycles  it  would  be  conceivable  for  it  to  undergo  several  rounds  of  toxin 
synthesis  before  dividing,  thereby  increasing  its  toxin  content  and  becoming  more  toxic. 
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Toxin  Composition 


Another  aspect  of  toxin  production  in  A.  fundyense  is  the  relative  composition  of  its 
toxin  pool.  It  has  been  reported  in  the  literature  that  the  toxin  composition  of  an  individual 
isolate  of  Alexandrium  remains  constant  under  a  variety  of  environmental  conditions  as 
well  as  through  the  different  culture  growth  stages  (Hall  1982;  Boyer  et  al.  1987;  Cembella 
et  al.  1987;  Ogata  et  al.  1987).  Recently,  however,  it  has  been  observed  that  significant 
changes  in  toxin  composition  can  occur  in  cells  exposed  to  different  stresses.  The  relative 
abundance  of  STX  in  Alexandrium  tamarense  and  A.  catenella  has  been  shown  to  decrease 
as  the  cells  enter  plateau  phase  and  stop  growing  (Boczar  et  al.  1988).  A.  fundyense 
shows  an  increase  in  STX  content  at  higher  specific  growth  rates  (Anderson  et  al.  1990a), 
which  is  consistent  with  Boczar  et  al.'s  results  (1988)  for  A.  tamarense  and  A.  catenella. 

Our  study  also  documents  changes  in  toxin  composition  in  A.  fundyense,  but  these 
variations  occurred  at  different  time  scales  from  those  in  previous  studies.  While  the 
relative  proportions  (mole  percent  of  total  toxin)  of  the  five  main  saxitoxins  remained 
essentially  unchanged  throughout  the  experiment  (data  not  shown),  the  accumulation 
patterns  (pmol  cell'^)  of  the  seven  main  derivatives  of  saxitoxin,  C2,  NEO,  STX,  GTX2,3 
isomers  and  GTX1,4  isomers,  showed  slight  variations  with  respect  to  the  timing  of 
formation  (Fig.  5).  Analysis  of  the  exact  onset  of  accumulation  for  each  of  the  toxins 
suggests  that  C2  starts  building  up  first,  followed  by  STX  and  the  isomeric  parr  GTX  2,3. 
With  a  delay  of  about  one  to  two  hours,  the  concentrations  of  NEO  and  the  pair  GTX1,4 
begin  to  increase.  A  plausible  interpretation  of  these  results  would  involve  a  biosynthetic 
pathway  consisting  of  one  precursor  centered  around  C2  and  several  subsequent 
transformation  steps  (Fig.  6).  First,  C2  would  be  transformed  into  GTX2,3  by  means  of 
an  A-sulfatase.  In  a  second  step,  GTX2,3  could  be  further  transformed  into  STX  by 
another  sulfatase.  Finally,  a  hydroxylase  adds  an  A-hydroxy  group  to  GTX2,3  and  STX, 
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Fig.  6  Putative  sequence  of  events  in  the  biosynthesis  of  SIX  and  its  derivatives. 

Continuous  lines  represent  most  likely  interconversions,  as  deduced  from  Fig. 
5  (see  discussion).  Dashed  lines  represent  other  possible  paths. 
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turning  them  respectively  into  GTX1,4  and  NEO.  The  rough  kinetics  of  the 
transformations,  as  derived  from  Fig.  5,  would  fit  in  such  a  putative  scheme.  These 
pathways  are  largely  speculative  at  this  stage.  The  synchronization  of  cultures  and  closely 
spaced  toxin  analysis  shows  promise  for  elucidating  the  kinetic  profiles  of  the  individual 
toxin  productions. 

A  similar  scheme  of  the  STX  biosynthetic  pathway  has  been  outlined  before  (Shimizu 
1993),  but  no  experimental  evidence  had  been  gathered  prior  to  our  study.  This 
interpretation  has  to  be  carefully  contrasted  with  observations  of  changes  and  variations  in 
toxin  composition  in  other  dinoflagellate  species  such  as  Pyrodinium  bahamense  and 
Gymnodinium  catenatum  (Oshima  et  al.  1993;  Usup  et  al.  1994).  In  these  organisms, 
STX  has  been  postulated  to  be  a  precursor  or  parent  compound  that  is  rapidly  converted 
into  NEO  or  GTX  V,  as  these  last  two  compounds  are  predominant  and  only  negligible 
amounts  of  STX  are  detected  (Oshima  et  al.  1993;  Usup  et  al.  1994).  In  an  integrated 
sheme  of  the  STX  biosynthetic  pathway,  interconversions  among  the  toxins  would  be 
determined  in  part  by  the  presence  or  absence  of  the  appropriate  enzymes,  but  actual  toxin 
production  would  be  the  result  of  the  accumulation  of  only  one  precursor  or  parent 
compound.  This  precursor  could  be  one  of  the  toxins  or  another  unknown  compound. 
Similar  studies  to  this  one  using  dinoflagellates  that  differ  in  their  overall  toxin  composition 
are  necessary  to  more  conclusively  elucidate  the  STX  biosynthetic  pathway. 

Cell  Cycle 

The  original  purpose  of  this  study  was  to  investigate  short  term  variability  in  toxicity  in 
A.  fundyense,  yet  it  also  revealed  some  interesting  aspects  of  the  cell  cycle  of  this  species 
and  the  manner  in  which  the  cells  respond  to  prolonged  darkness.  DNA  synthesis  in 
dinoflagellates  has  been  reported  to  be  continuous  throughout  the  cell  cycle  in  several 
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species  (Karentz  1983;  Rizzo  1987),  or  confined  to  a  discrete  S  phase  in  others  (Chang  and 
Carpenter  1988;  Yamaguchi  1992).  A.  fundyense  clearly  sinthesyzes  DNA  during  a 
discrete  time  period,  phased  to  a  specific  portion  of  the  cell  cycle  (Fig.  4).  This  S  phase 
takes  place  during  the  later  portion  of  the  light  period  and  extends  into  the  early  dark  period 
(Fig.  7). 

An  interesting  observation  relates  to  the  presence  of  two  cohorts  of  cells  at  the  onset  of 
the  experiment  and  the  subsequent  homogeneity  of  the  population  several  light:  dark  cycles 
later.  This  can  be  explained  by  the  presence  of  two  transition  points  in  Gi.  Cells  arrested 
at  the  first,  light-dependent  transition  point  (in  early  Gi,  immediately  following  mitosis) 
proceed  through  a  Go-like  period  associated  with  lower  metabolic  activity  and  smaller  size 
(Figs.  2  and  5).  In  addition,  a  small  portion  of  the  population  (approx.  20%)  was 
apparently  arrested  at  a  transition  point  located  further  into  Gi  and  in  this  case  size- 
dependent.  These  cells  were  slightly  larger  than  the  mean  of  the  cells  at  the  beginning  of 
the  experiment  (Fig.  2)  and  account  for  the  surge  in  S  phase  and  G2  +  M  cells  and  the 
concurrent  increase  in  cell  concentration  during  the  dark  period  of  the  first  light:dark  cycle 
(Figs.  lA  and  3).  Subsequently  the  entire  population  doubled  over  the  next  light:dark 
cycle.  This  means  that  all  the  cells  had  become  synchronized  and  were  proceeding  through 
the  cell  cycle  as  a  single  unified  cohort.  In  two  subsequent  light:  dark  cycles  cell  densities 
increased  by  a  factor  of  1.8  and  1.7  respectively  (Fig.  lA).  This  might  indicate  a  gradual 
loss  of  synchronicity  in  the  culture,  but,  the  cell  cycle  information  showed  in  Fig  3 
indicates  that  the  level"  of  synchronicity  was  still  at  approximately  96%.  Another 
explanation  for  the  lack  of  complete  doubling  in  cell  numbers  could  be  the  loss  of  an 
undetermined  number  of  cells  from  the  sampled  population.  From  the  second  day  on  we 
could  observe  a  cell  deposit  at  the  bottom  of  all  three  carboys  that  was  difficult  to  resuspend 
without  damaging  the  cells. 
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Fig.  7  Diagrammatic  representation  of  cell  cycle  events  in  Alexandrium  fundyense 
Balech  (GtCA29)  cells  during  synchronous  growth  on  a  light:  dark  cycle  of 
14: 10  hours.  Open  bar:  light  phase;  Gi:  gap  1  (post  mitotic);  S:  DNA  synthesis 
phase;  G2:  gap  2  (pre  mitotic);  M:  mitosis. 
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Single,  light-dependent  transition  points  in  Gi  have  been  identified  in  Chlamydomonas 
(Spudich  and  Sager  1980)  and  in  Amphidinium  carterii  and  Hymenomonas  carterii  (Olson 
and  Chisholm  1986;  Vaulot  et  al.  1986).  Two  light-dependent  transition  points  have  been 
characterized  in  Thalassiosira  weisflogii,  one  in  Gi  and  a  second  one  in  G2  (Vaulot  et  al. 
1986).  A  single  nutrient  dependent  control  point  located  in  Gi  has  also  been  identified  in 
Prorocentrum  minimum  (Antia  et  al.  1990)  and  the  simultaneous  presence  of  nutrient 
dependent  control  points  in  Gi  and  in  G2  has  been  described  in  Oxyrrhis  marina  (Whiteley 
et  al.  1993).  To  our  knowledge,  the  data  presented  here  are  the  first  evidence  in 
phytoplankton  suggesting  that  there  may  be  two  different  transition  points  in  Gi,  one  of 
them  light-,  the  other  one  size-dependent. 

An  alternative  explanation  for  the  existence  of  two  cohorts  of  cells  at  the  beginning  of 
the  experiment  that  does  not  require  two  transition  points  in  Gi  is  suggested  by  the  cell 
cycle  model  of  Klevecz  (Klevecz  1976),  who  noted  that  certain  mammalian  cell  lines  had 
generation  times  that  tended  to  be  integer  multiples  of  4  hrs.  In  the  Klevecz  model,  Gi  is 
replaced  by  a  loop  called  Gq  which  has  a  duration  equal  to  the  quantized  intervals  between 
generation  times.  Gq  can  thus  be  viewed  as  a  subcycle  from  which  the  cells  can  reenter  the 
remainder  of  the  cell  cycle,  but  only  after  completely  traversing  Gq  one  or  more  times. 
Since  exit  from  Gq  is  regulated  by  environmental  factors  such  as  light,  the  net  effect  of 
prolonged  darkness  would  be  for  a  population  of  cells  to  be  separated  into  two  or  more 
groups  depending  upon  the  number  of  Gq  subcycles  they  had  traversed  when  darkness  was 
imposed.  Evidence  for  quantized  generation  times  in  the  diatom  Thalassiosira  fluvialitis 
has  been  described,  but  the  distribution  of  generation  times  for  that  species  was  not 
consistent  with  the  Klevecz  model  (Chisholm  and  Costello  1980).  In  the  case  of  A. 
fundyense,  such  a  model  cannot  be  dismissed,  but  the  data  presented  here  would  imply  a 
Gq  of  24  hours.  Therefore,  a  model  of  the  cell  cycle  that  integrates  a  light-dependent  and  a 
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size-dependent  transition  points  seems  to  be  the  easiest  way  to  explain  our  observations 
(Fig.  7). 

Summary 

The  present  study  documents  for  the  first  time  the  dynamics  of  toxin  production  during 
the  cell  cycle  of  Alexandrium  fundyense  (Fig.  7).  Toxin  production  was  shown  to  occur 
over  a  discrete  period  during  Gi  and  to  be  induced  by  light.  The  cell  cycle  contains  two 
transition  points  responsible  for  regulating  the  passage  of  the  cell  through  Gi  (light- 
dependent),  and  the  transition  from  Gi  to  S  (size-dependent),  respectively.  We  also 
suggest  that  C2  is  the  first  derivative  accumulating  in  the  biosynthetic  pathway  of  the  toxins 
in  our  A.  fundyense  isolate,  being  subsequently  transformed  into  the  other  toxins.  The 
restriction  of  toxin  production  to  a  relatively  short  segment  of  the  cell  cycle  provides  a  tool 
for  comparing  cells  that  are  and  are  not  synthesizing  toxin.  This  has  important  implications 
with  respect  to  the  search  for  genes  or  enzymes  involved  in  toxin  production  and  their 
regulation  in  Alexandrium.  In  this  context,  work  is  in  progress  in  our  laboratory  to  identify 
toxin-specific  genes  by  Differential  Display  analysis  (Liang  and  Pardee  1992)  of 
messenger-RNA  components  of  cells  harvested  at  different  points  during  their  cell  cycle. 
In  general,  this  work  presented  here  shows  how  studying  dinoflagellates  at  a  high  temporal 
resolution  and  under  synchronized  conditions  can  be  useful  in  order  to  gain  a  broader 
understanding  of  their  biochemistry  and  physiology. 
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CHAPTER  II 


Coupling  of  saxitoxin  biosynthesis  to  the  Gi  phase  of  the 
cell  cycle  in  the  dinoflagellate  Alexandrium  fundyense: 
temperature  and  nutrient  effects 
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ABSTRACT 


The  linkage  between  changes  in  the  length  of  cell  cycle  stages  and  the  toxicity  of 
Alexandrium  fundyense  Balech  were  studied  in  semi-continuous  cultures.  Growth  rates 
ranging  from  0.031  d'*  to  0.36  d  ^  were  established  at  different  temperatures  or  levels  of 
phosphate  limitation.  In  all  temperature  and  phosphate  limitation  treatments,  Gi  was 
always  the  phase  with  the  longest  duration.  Decrease  in  growth  rate  was  associated  with 
an  increase  in  duration  of  the  different  cell  cycle  stages,  Gi  and  G2  +  M  showing  the  most 
variation  relative  to  the  respective  increase  in  generation  time  and  S  the  stage  showing  the 
least.  Over  the  range  of  temperatures  studied  (20”C  8"C),  the  duration  of  Gi  doubled, 

that  of  G2  +  M  increased  by  a  factor  of  2.7  and  that  of  S  by  3.8.  Over  the  range  of  growth 
rates  studied  under  phosphate  limitation  (0.36  0.031),  the  duration  of  Gj  increased  by  a 

factor  of  5.3,  that  of  G2  +  M  by  43  and  that  of  S  by  10.  The  temperature  treatments 
demonstrated  a  direct  correlation  between  length  of  the  cell  cycle,  and  more  particularly, 
duration  of  the  Gi  phase,  with  toxin  content  (=  toxin  cell  quota).  Under  phosphate- 
limiting  conditions  the  same  correlation  with  Gi  duration  was  observed,  along  with  a  clear 
uncoupling  of  toxin  accumulation  from  the  S  and  G2  phases  of  the  cell  cycle.  In  both  the 
temperature  treatments  and  the  phosphate  limitation  experiments,  toxin  production  rates 
remained  constant  for  the  respective  range  of  conditions,  implying  that  the  variations  in 
toxin  content  observed  were  a  result  of  increasing  periods  of  biosynthetic  activity.  Toxin 
accumulation  was  directly  correlated  to  protein  biosynthesis  in  all  temperature  treatments. 
In  contrast,  toxin  content  showed  little  correlation  with  protein  content  as  phosphate 
limitation  increased.  Significant  differences  in  toxin  composition  were  observed  between 
the  temperature  and  phosphate  treatments;  Total  concentrations  of  GTX II  and  III  and  C I 
and  II  were  significantly  higher  in  the  phosphate-limited  cultures,  while  the  levels  of  STX, 
NEO  and  gonyautoxins  I  and  IV  remained  virtually  unchanged.  We  conclude  that  toxin 
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biosynthesis  in  A.  fundycnse  is  coupled  to  the  Gj  phase  of  the  cell  cycle,  that  toxin 
synthesis  is  not  down-regulated  by  phosphate  deprivation  and  that  interconversions  among 
saxitoxin  derivatives  are  influenced  by  the  availability  of  phosphate.  These  data  confirm 
earlier  observations  that  toxin  production  is  discontinuous  over  the  cell  cycle  of  A. 
fundyense  and  shed  new  light  on  the  role  of  phosphate  in  the  biosynthetic  pathway. 
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INTRODUCTION 


Dinoflagellates  in  the  genus  Alexandrium  are  responsible  for  outbreaks  of  Paralytic 
Shellfish  Poisoning  (PSP)  in  temperate  coastal  areas  around  the  world  (Anderson,  1989). 
PSP  is  caused  by  saxitoxin  (STX)  and  its  over  twenty  known  derivatives,  a  group  of 
highly  nitrogenous  alkaloids  with  a  marked  affinity  and  blocking  capacity  for  sodium 
channels  (Shimizu,  1993).  STX  is  biosynthesized  by  Alexandrium  sp.  and  subsequently 
accumulates  in  bivalves  filter-feeding  on  the  dinoflagellate;  ingestion  of  these  toxic  shellfish 
causes  PSP  symptoms  at  higher  levels  of  the  trophic  chain,  including  humans  (Sommer  et 
al.,  1937).  Knowledge  of  the  mechanism  of  toxin  biosynthesis  and  its  regulation  in 
dinoflagellates  is  crucial  for  a  complete  understanding  of  this  phenomenon. 

Considerable  work  has  been  done  on  the  influence  of  growth-limiting  physico-chemical 
factors  such  as  nutrient  concentrations,  light  and  temperature  on  the  ecophysiology  of  this 
genus  and  its  capacity  to  produce  toxins  (Proctor  et  al.,  1975,  White,  1978,  Hall,  1982, 
Boyer  eta/.,  1987;  Ogata  eta/.,  1987;  Anderson  eta/.,  1990b;  Flynn  eta/.,  1994;  Matsuda 
et  a/.,  1996).  Under  nitrogen  starvation,  light  limitation  or  non-optimal  salinity  conditions, 
toxin  content  (=  toxin  cell  quota)  is  reduced  concomitantly  with  a  decrease  in  the 
availability  of  the  limiting  factor.  Only  in  the  case  of  phosphate  starvation  and  temperature 
limitation  has  the  contrary  been  observed  -  an  increase  in  toxin  content  with  increasing 
stress.  Phosphorus  is  incorporated  into  four  main  organic  compound  fractions  in  the  cells: 
RNA,  DNA,  and  phosphorus  included  in  lipids  and  in  esters  that  have  a  direct  bearing  on 
cell  metabolism.  Furthermore,  phosphates  play  an  important  metabolic  role  as  substrates 
for  the  fine  regulation  of  biochemical  activity  in  the  cells  via  protein  kinases  (Pardee  et  al, 
1978).  Temperature  is  a  general  modulator  of  growth  and  biosynthetic  rates  in 
phytoplankton  and  can  therefore  influence  the  accumulation  of  specific  metabolites  (Jitts  et 
a/.,  1964). 
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In  addition  to  changes  in  the  toxin  content,  toxin  composition  (i.e.,  the  absolute  or 
relative  amounts  of  the  different  STX  derivatives  present  in  the  organism)  has  also  been 
shown  to  vary  under  different  limiting  conditions  (Boczar  et  al.,  1988;  Anderson  et  al., 
1990a;  MacIntyre  et  al,  1997).  These  variations  are  a  result  of  physiological  changes  that 
regulate  the  sequence  of  transformations  leading  to  the  different  derivatives.  These 
pathways  have  not  yet  been  elucidated.  Only  recently  have  several  enzymes  involved  in 
specific  interconversions  among  STX  derivatives  been  characterized,  but  the  information  is 
still  fragmented  with  respect  to  the  sequence  of  interconversions  (Oshima,  1995;  Yoshida 
et  al,  1996).  Observation  of  accumulation  patterns  for  the  different  STX  derivatives  in 
synchronized  cultures  of  A.  fundyense  suggests  a  sequence  and  directionality  to  these 
interconversions  beginning  with  C  I  and  II  (Taroncher-Oldenburg  et  al,  1997),  but  no 
specific  pathway  has  yet  been  outlined,  as  the  direct  precursor  or  parent  compound  for  the 
series  of  the  STXs  has  not  been  determined.  The  sequence  of  interconversions  cannot  be 
rooted  at  a  specific  origin  or  parent  compound  with  our  present  knowledge  (Shimizu, 
1993). 

In  parallel  to  the  influence  of  the  above-mentioned  environmental  factors  on  toxin 
production,  a  more  immediate  effect  of  limiting  conditions  has  been  observed  on  the 
growth  rate  of  toxic  dinoflagellates  (Proctor  et  al,  1975;  White,  1978;  Hall,  1982;  Watras 
et  al,  1982;  Anderson  et  al,  1990b;  Matsuda  et  al,  1996).  The  general  trend  in  this  case 
is  a  reduction  in  growth  rate  with  an  increase  in  the  limitation  or  stress.  It  has  further  been 
documented  for  several  groups  of  phytoplankton  that  such  changes  in  growth  rate  are 
tightly  correlated  to  changes  in  the  duration  of  the  different  cell  cycle  stages  (Vaulot,  1985; 
Olson  and  Chisholm,  1986;  Cetta  and  Anderson,  1990).  Depending  on  the  organism, 
changes  in  the  doubling  time  can  be  reflected  in  changes  in  the  duration  of  specific  phases 
of  the  cell  cycle  (Guiget  et  al,  1984;  Olson  et  al,  1986;  Yamaguchi  and  Honjo,  1989). 
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The  eukaryotic  cell  cycle  is  divided  into  four  main  stages.  Following  the  mitotic 
division  of  the  mother  cell,  the  two  daughter  cells  undergo  a  first  period  called  Gi,  during 
which  they  synthesize  essential  cell  constituents  and  grow  in  size.  Once  appropriate 
internal  and  external  conditions  are  met  the  cell  commits  to  a  new  division.  This  process 
follows  three  steps:  First,  the  cells  replicate  their  DNA,  thereby  doubling  their  genetic 
information,  during  the  S  phase.  A  second  metabolically  active  period  follows,  called  G2, 
in  preparation  for  the  final  phase,  mitosis,  M,  during  which  cells  undergo  cytokinesis  and 
divide  their  genetic  information  among  their  offspring.  Progression  through  the  cell  cycle 
is  highly  influenced  by  environmental  conditions,  as  these  determine  the  metabolic  and 
energetic  status  of  the  cell  (Yamaguchi  and  Honjo,  1989).  Cells  will  grow  slower  if  their 
metabolism  is  not  optimal  and  a  successful  division  process  is  compromised  (Pardee  et  al, 
1978).  Conversely,  metabolic  processes  regulated  by  cell  cycle  events  will  also  be 
affected. 

In  algae,  biochemical  pathways  such  as  DNA  synthesis  and  a  variety  of  specific 
enzyme  activities  are  modulated  in  a  periodic,  cell  cycle  dependent  fashion  (Kates  and 
Jones,  1967;  Vassef  er  al,  1973;  Forde  and  John,  1974).  These  biosynthetic  pathways  are 
influenced  by  conditions  affecting  normal  progression  through  the  cell  cycle  (Mitchison, 
1969;  Mitchison  et  al,  1991).  The  lengthening  of  a  specific  cell  cycle  stage  will  thus  result 
in  longer  metabolic  activities  and  hence  higher  product  accumulations.  In  the  case  where 
the  pathway  of  interest  is  only  influenced  by  the  cell  cycle  dynamics  and  no  other  factors, 
this  would  translate  into  a  constant  ratio  between  the  concentration  of  the  metabolite  and  the 
length  of  the  cell  cycle  phase  of  interest.  This  is  equivalent  to  the  existence  of  a  constant 
cell  cycle-specific  production  rate  of  the  compound. 

Previous  studies  have  shown  an  inverse  correlation  between  growth  rate  and  toxin 
content  in  A.  tamarense  (Proctor  et  al,  1975;  Ogata  et  al,  1987)  and  a  direct  correlation 
between  specific  growth  rate  and  specific  toxin  production  rate  in  the  closely  related 


80 


organism  A.  fundyense  (Anderson  et  al,  1990b).  We  have  also  shown  that  toxin 
production  is  discontinuous  and  restricted  to  a  portion  of  the  Gi  phase  of  the  cell  cycle  in 
synchronously  growing  A.  fundyense  (Taroncher-Oldenburg  et  al.,  1997).  These  results 
lead  us  to  the  conclusion  that  there  should  be  a  correlation  between  the  rates  of  toxin 
production  or  toxin  accumulation  and  length  of  Gi.  The  aim  of  this  study  was  to  determine 
if  there  is  any  direct  correlation  between  variations  in  toxin  production,  measured  as 
accumulation  of  saxitoxin  in  a  cell,  and  changes  in  duration  of  the  different  stages  of  the 
cell  cycle  under  different  conditions  that  affect  growth  rate  and  thus  cell  cycle  progression. 
Such  a  correlation  would  confirm  that  STX  biosynthesis  is  discontinuous  and  restricted  to 
the  Gi  cell  cycle  stage. 


MATERIALS  AND  METHODS 

The  organism.  The  dinoflagellate  Alexandrium  fundyense  (strain  GtCA29)  was  used 
throughout  this  study.  A  culture  was  established  in  January  1985  from  a  cyst  isolated  from 
Gulf  of  Maine  sediments,  20  miles  east  of  Portsmouth,  New  Hampshire  (USA).  That 
culture  was  used  to  isolate  a  single  cell  to  establish  a  clonal  culture.  Prior  to  the 
experiments,  cultures  were  maintained  at  15 °C  in  f/2  medium  (Guillard  and  Ryther,  1962) 
modified  by  addition  of  H2Se03  (10"®M)  and  by  reducing  the  concentration  of  CUSO4  to 
10‘^M.  Vineyard  Sound  seawater  (0.2  pm  filtered,  31°/oo  salinity)  was  used  as  the 
medium  base.  Throughout  the  experiments,  irradiance  of  ca.  250  pEm'^s'^  was  provided 
by  cool  white  fluorescent  bulbs  on  a  14:10  h  light:dark  cycle. 

Temperature  experiment  with  semi-continuous  cultures.  Effects  of  temperature  on 
toxin  content,  toxin  composition,  protein  content  and  cell  cycle  phase  lengths  were  studied 
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using  a  temperature-gradient  bar  (Watras  et  al,  1982),  in  which  a  range  of  seven  different 
temperatures  (8”C,  10°C,  ITC,  13°C,  15“C,  17°C  and  20*C)  was  established  (Table  1). 
Twenty  eight  25  x  150  mm  culture  tubes  containing  25  ml  of  sterile  filtered  f/2  medium 
were  inoculated  with  exponentially  growing  cells  and  acclimated,  four  per  treatment,  to 
each  temperature  over  two  transfer  cycles.  Growth  was  monitored  throughout  the 
experiment  using  in  vivo  fluorescence  measurements  obtained  on  a  Turner  Designs  model 
10  fluorometer  (Turner  Designs,  Sunnyvale,  Ca).  Semi-continuous  cultures  were 
established  after  specific  growth  rates  were  determined  for  the  different  temperatures. 
Following  the  second  transfer,  these  rates  were  used  as  the  respective  daily  dilution  rates  of 
fresh  medium  for  each  treatment  in  order  to  maintain  a  constant  cell  density  in  the  cultures 
over  a  period  of  twelve  days.  The  cell  densities  ranged  from  6,000  cells  ml'^  at  8°C  to 
about  8,000  cells  mT^  at  20°C.  At  the  end  of  the  experiment,  the  four  tubes  per 
temperature  were  combined,  two  by  two.  The  resulting  two  samples  per  treatment  were 
then  harvested  and  aliquots  prepared  for  toxin,  DNA  and  protein  analysis. 

Low  phosphate  experiment  with  semi  continuous  cultures.  To  study  the  effect  of 
phosphate  limitation  on  toxin  content,  toxin  composition,  protein  content  and  cell  cycle, 
semi-continuous  cultures  of  A.  fundyense  were  established  at  20°C  and  in  f/2  medium  in 
which  the  concentration  of  P04^'  was  lowered  from  36|iM  to  1.8|iM  (f/40  [P04^  ]).  Five 
different  dilution  rates  were  applied  to  the  different  cultures,  four  per  treatment,  in  order  to 
provide  five  different  phosphate  limited  growth  rates  (0.5,  0.25,  0.1,  0.075  and  0.05 
divisions  per  day)  (Table  1).  Growth  was  monitored  throughout  the  experiment  using  in 
vivo  fluorescence  measurements  obtained  on  a  Turner  Designs  model  10  fluorometer. 
After  two  transfers  the  cultures  were  kept  at  a  cell  density  of  approximately  3,000  cells  ml'* 
for  a  twelve  day  period,  following  which,  the  four  tubes  per  dilution  rate  were  combined. 
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two  by  two.  The  resulting  two  samples  per  treatment  were  then  harvested  and  aliquots 
prepared  for  toxin,  DNA  and  protein  analysis. 

Flow  cytometric  analysis  of  DNA.  Ten  ml  aliquots  of  culture  were  preserved  with  5% 
formalin  (v/v)  and  stored  at  4°C.  The  cells  were  centrifuged  (1700  x  g,  23°C,  3  mins)  and 
the  pellet  resuspended  in  10  ml  of  ice  cold  methanol  and  stored  at  4°C  to  extract 
chlorophyll.  The  samples  were  then  further  concentrated  by  centrifugation.  A  subsample 
of  the  methanol-preserved  material  was  transferred  into  a  1.5  ml  microcentrifuge  tube  and 
pelleted  at  2000  rpm.  The  supernatant  was  removed  and  the  cells  washed  with  1  ml  PBS 
(40mM  Na2HP04,  22mM  KH2PO4,  85mM  NaCl).  After  centrifugation,  the  pellet  was 
resuspended  in  0.5  ml  propidium  iodide  (PI)  staining  solution  (4  pg  ml"^  PI,  1250  units 
RNaseA  ml'^)  and  allowed  to  react  in  the  dark  for  at  least  two  hours  prior  to  analysis  on  an 
Epics  V  flow  cytometer  (Coulter  Electronics,  Hialeah,  FL)  (Olson  et  al,  1986).  The 
proportions  of  cells  in  Gi,  S  and  G2  +  M  were  determined  as  described  elsewhere  (Cetta 
and  Anderson,  1990).  The  duration  of  the  cell  cycle  stages  was  calculated  using  the 
equations  of  Slater  et  al  (Slater  et  al. ,  1977); 

r(Gi)  =-(tD/ln2)ln(l-[P(Gi)/2]) 
r(S)  =  (tD/ln2)  In  (1  +  P(S)/[1  +  P(G2)]) 
r(G2)  =  (tD/ln2)  In  (1  +  P(G2)) 

where  r(x)  is  the  duration  of  each  cell  cycle  stage,  to  represents  the  doubling  time  of  the 
population  under  balanced  growth  conditions  and  P(x)  is  the  proportion  of  cells  in  each 
specific  cell  cycle  stage. 
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Toxin  Analysis.  Between  1.5  and  2.5  x  10^  cells  were  collected  by  centrifugation 
(3,000  X  g,  23°C).  The  resulting  pellet  was  rinsed  into  a  centrifuge  tube  and  resuspended 
in  0.5  M  acetic  acid  for  extraction  by  sonification.  The  extracts  were  stored  at  -20°C  prior 
to  analysis.  The  extracts  were  analyzed  by  HPLC  (Oshima  et  al,  1989)  incorporating  the 
modifications  previously  described  by  Anderson  et  al.  (Anderson  et  al,  1994). 

Soluble  Protein.  Approximately  30,000  cells  were  collected  by  centrifugation  (3000  x 
g,  23 °C)  for  protein  quantification.  After  aspiration  of  the  supernatant,  the  pellet  was 
stored  frozen  at  -20‘’C.  For  analysis,  1  ml  of  0.1  N  NaOH  was  added  to  the  pellet  and  the 
suspension  sonified  for  15  seconds  to  disrupt  the  cells.  Samples  were  then  digested  for  20 
min  at  80°C,  centrifuged,  and  the  supernatant  removed  and  neutralized  with  0.1  volumes  of 
1.0  N  HCl  (Binder  and  Anderson,  1990).  A  50  pi  fraction  of  the  neutralized  solution  was 
mixed  with  1  ml  of  bicinchoninic  acid  (BCA)  working  solution  (Pierce,  Rockford  IL),  and 
the  absorbance  of  duplicate  samples  determined  at  562  nm. 

Abbreviations  used  throughout  this  text  are:  STX  =  saxitoxin;  NEO  =  neosaxitoxin; 
GTX  1,2,3,4,5  =  gonyautoxins  I,  II,  III,  IV  and  V;  C2  =  toxin  C2;  Gq,  Gi  and  G2  =  Gap 
0,  Gap  1  and  Gap  2  phases  of  the  cell  cycle;  S  =  synthesis  phase;  M  =  mitosis. 


RESULTS 

Growth  and  Cell  Cycle  Dynamics.  The  growth  rate  of  A.  fundyense  decreased  as 
temperatures  decreased  (Table  1).  The  slope  obtained  from  a  linear  regression  of  grovAh 
rate  against  temperature  (a=-1.65;  =  0.93)  matched  the  slope  of  a  theoretical  curve  (a=- 
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Table  1 .  Growth  rates  of  Alexandriumfundyense  grown 
in  semi-continuous  culture  under  phosphate 
limitation  or  at  different  temperatures. 


Temperature  f 

(°C) 

PO3'  limited 
cultures  $ 

dil.  rate  (day^) 

Growth  rate  § 

p  (day-1) 

20 

0.277 

17 

15 

13 

0.186 

11 

0.18 

10 

0.15 

8 

0.121 

20 

0.35 

0.36 

20 

0.17 

0.17 

20 

0.069 

20 

0.052 

20 

0.035 

f  The  daily  dilution  rate  was  equivalent  to  the  growth  rate  of  exponentially 
growing  cultures  at  the  same  temperature 


:jl  Cultures  were  diluted  once  a  day  with  low  phosphate  medium 
§  In  all  treatments  the  standard  deviation  obtained  from  the  four  replicates  was 
less  than  8%  of  the  mean  value 
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1.67)  representing  a  doubling  of  growth  rate  with  a  10°C  increase  in  temperature  (i.e.,  a 
Qio  of  2;  figure  not  shown).  Decrease  in  growth  rate  was  associated  with  an  increase  in 
duration  of  the  different  cell  cycle  stages,  Gi  and  G2  +  M  showing  the  most  variation 
relative  to  the  respective  increase  in  generation  time  and  S  the  stage  showing  the  least  (Fig. 
lA). 

In  the  phosphate-limited  cultures,  growth  rates  decreased  as  phosphorus  supply  rates 
decreased  due  to  reductions  in  the  dilution  rate  of  the  cultures  (Table  1).  Discrepancies 
between  the  dilution  rates  and  actual  growth  rates  are  a  result  of  variations  associated  with 
pipetting  the  small  dilution  volumes.  Only  the  empirically  determined  growth  rates  were 
subsequently  used  for  calculations.  Similar  to  the  temperature  experiment,  decreases  in 
growth  rate  also  resulted  in  a  lengthening  mostly  of  Gj  and  G2  M  relative  to  the  increase 
in  generation  time  (Fig.  IB). 

In  all  temperature  and  phosphate  limitation  treatments,  Gi  was  always  the  phase  with 
the  longest  duration.  The  length  of  Gi  was  comparatively  longer  in  the  phosphate 
depletion  experiment  than  in  the  temperature  treatments  (Fig.  1 A  &  B).  As  an  example,  at 
irC  and  with  a  specific  growth  rate  of  0.18,  Gi  was  59  hours,  G2  +  M  28  hours  and  S 
4.7  hours  long,  while  at  a  growth  rate  of  0.17  under  phosphate  limitation,  Gi  was  79 
hours,  G2  +  M  15  hours  and  S  1.8  hours  long  (Fig.  1). 

Variations  in  Toxin  Content.  Toxin  content  was  higher  at  lower  growth  rates  for  both 
experiments  and  was  higher  under  phosphate  limitation  compared  to  the  replete  medium 
conditions  (Fig.  2A).  In  the  temperature  experiment,  and  over  the  range  of  temperatures 
studied  (20°C  8°C),  the  toxin  cell  quota  increased  by  a  factor  of  2.1.  Over  the  range  of 

growth  rates  studied  under  phosphate-limiting  conditions  (0.36  -♦  0.031),  the  toxin  cell 
quota  increased  by  a  factor  of  4.8.  Toxin  content  under  phosphate  limitation  was 
approximately  four  times  as  high  as  in  some  of  the  treatments  in  the  temperature  experiment 
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at  the  same  growth  rate.  These  toxin  levels  and  differences  in  toxicity  between  nutrient- 
replete  and  phosphate-limited  treatments  coincide  with  observations  by  other  authors 
(Boyer  et  al,  1987;  Anderson  et  al.,  1990b). 

Variations  in  Protein  Concentrations.  Figure  2B  shows  the  variability  in  protein 
concentrations  for  cells  grown  at  different  temperatures.  Similar  to  what  was  observed  for 
toxin  content  (Fig.  2A),  the  levels  of  protein  increased  as  growth  rate  decreased.  A 
different  pattern  was  revealed  in  the  phosphate-limited  cultures  (Fig.  2B).  In  this  case  the 
protein  levels  decreased  at  lower  growth  rates.  At  rates  between  0.2  and  0.4  day'^  the 
levels  of  protein,  649±87  pg  cell'^  were  similar  between  the  two  experimental  treatments 
and  coincided  with  values  previously  reported  in  the  literature  (Boyer  et  al,  1987; 
Anderson  etal,  1990b). 

Correlation  between  Toxin  and  Protein  Contents.  The  ratios  of  toxin  content  to  protein 
per  cell  were  determined  for  each  treatment.  In  the  temperature  experiment  this  ratio  was 
constant  at  a  value  of  0.07910.006  fmol  of  toxin  per  pg  of  protein  (Fig.3).  The  slope  of 
the  linear  regression  through  the  data,  a=0.013,  was  not  statistically  different  from  zero  (P 
>  0.95,  C.I.  93%).  This  ratio  was  not  constant  in  the  phosphate-limited  cultures  but 
increased  as  growth  rate  decreased  (P  <  0.005)  (Fig.3). 

Correlation  between  Toxin  Contents  and  Cell  Cycle  Phase  Durations.  Cell  cycle  phase 
durations  and  toxin  content  were  normalized  in  both  experiments  against  the  lowest  values 
of  the  respective  data  series.  These  normalized  values  showed  a  clear  overlap  in  the 
relative  changes  of  Gi  duration  and  toxin  content  (Fig.4).  A  linear  regression  through  the 

r\ 

normalized  data  resulted  in  a  slope  of  a  =  0.97  (R  =  0.98).  In  both  experiments  the 
relative  increases  in  duration  of  G2  and  S  are  different  from  the  concurring  relative  changes 
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Fig.  1.  Cell  cycle  phase  durations  in  Alexandrium  fundyense  at  different  growth 
rates  in  semi-continuous  culture,  under  nutrient  replete  or  phosphate-limited 
conditions.  (Gi  T,  G2  +  M  □,  S  •).  (A)  Cells  grown  in  f/2  medium  at 
seven  different  temperatures  ranging  from  8  to  20'C  (see  Table  1  for 
correspondence  of  p  with  temperature  of  treatment).  (B)  Cells  grown  at 
different  levels  of  phosphate  starvation  (see  Table  1  for  correspondence  of  p 
with  phosphate  depletion  regimes). 
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Fig.  2.  Changes  in  (A)  total  toxin  levels  and  (B)  protein  concentrations  in 
Alexandrium  fundyense  grown  in  semi-continuous  culture  under 
temperature-limited  conditions  in  f/2  medium  (•)  or  under  phosphate-limited 
conditions  (O)  (Error  bars  represent  SE). 
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Fig.  3.  Ratio  of  toxin/protein  contents  in  Alexandriumfundyense  at  different  growth 
rates  in  semi-continuous  culture  and  under  temperature-limited  conditions  in 
f/2  medium  (•)  or  under  phosphate-limited  conditions  (O). 
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Fig.  4.  Relative  changes  in  cell  cycle  stage  durations  of  Gi  (T),  G2  +  M  (□)  and  S 
(•)  against  toxin  content,  for  Alexandrium  fundyense  in  semi-continuous 
culture  (dashed  line  represents  1:1  correlation).  Data  from  experiments 
under  temperature-limited  conditions  in  f/2  medium  and  under  phosphate- 
limited  conditions  are  pooled.  Data  were  normalized  against  the  lowest  value 
of  the  corresponding  series  (Fig.2). 


94 


in  toxin  content.  In  this  case  the  slopes  were  a  =  0.46  (R^  =  0.97)  and  a  =  0.1  (R^  -  0.88) 
respectively. 

Toxin  Production  Rates.  Four  different  net  toxin  production  rates  (fmol  cell’^  hour’^) 
were  calculated  for  each  of  the  experimental  treatments  (Note  that  production  rate  is  used  as 
an  equivalent  term  for  toxin  accumulation  since  we  do  not  know  the  extent  of  toxin  leakage 
or  catabolism).  First,  the  overall  net  toxin  production  rates  were  determined  (Anderson  et 
al,  1990b;  Taroncher-Oldenburg  etal,  1997).  Under  steady-state  conditions  the  specific 
rate  of  toxin  production  must  equal  the  specific  growth  rate  of  the  organism,  resulting  in  a 
constant  toxin  concentration  per  cell.  Hence,  the  net  toxin  production  rate  per  cell  becomes 
a  function  of  the  doubling  rate  of  the  organism  and  the  steady  state  toxin  concentration. 
This  can  be  expressed  with  the  following  equations 

A^Tox  ~ 

where  ^Utox  (hr  i)  is  the  specific  toxin  production  rate  and  ii  the  specific  growth  rate  (hr  i); 

Rtox  =  k  [Tox]  =  (p/lnl)  [Tox]  =  (AtTox/ln2)  [Tox] 

where  Rtox  (fmol  cell'i  hri)  is  the  net  toxin  production  rate,  k  (hr^)  is  the  doubling  rate 
and  [Tox]  (fmol  cell'i)  represents  the  toxin  cell  quota. 

We  further  determined  the  three  cell  cycle  stage-specific  net  toxin  production  rates,  i.e., 
the  rates  at  which  toxin  would  be  synthesized  if  all  production  occured  exclusively  during 
one  of  the  three  phases  of  the  cell  cycle  (Gi,  S  or  G2  +  M).  These  rates  are  described  by 
the  following  equation 
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Rtox(x)  =  [Tox]/r(x) 


where  Rtox(x)  is  the  net  toxin  production  rate  for  each  cell  cycle  stage  (x  =  Gi,  S  or  G2  + 
M)  and  T(x)  is  the  duration  of  each  cell  cycle  stage. 

The  values  of  Rtox  (fmol  cell'^  hour'^  ±  SD)  remained  constant  at  a  mean  value  of 
0.77+0.12  for  all  of  the  experimental  temperatures  (Fig.  5 A).  The  slope  of  the  linear 
regression  through  the  data,  a=-0.01,  was  not  statistically  different  from  zero  (P  >  0.94, 
C.I.  96%).  A  similar  pattern  was  observed  for  the  five  phosphate  depletion  regimes  that 
resulted  in  a  mean  value  of  3.34±0.20  (Fig.  4B).  A  linear  regression  through  the  data  also 
resulted  in  a  slope,  a=0.06,  that  was  not  significantly  different  from  zero  (P  >  0.89,  C.I. 
92%).  The  results  for  Rtox(x)  (fmol  celf^  hour'^  ±  SD)  for  both  experiments  are  plotted  in 
Fig.  5A  &  B.  Rtox(Gi)  was  constant  with  a  mean  value  of  0.94±0.21  in  the  temperature 
experiment  (Fig.  5A)  and  3.56+0.43  in  the  phosphate  depletion  treatments  (Fig.  5B). 
Rtox(S)  and  Rtox(G2  +  M)  showed  variations  within  and  between  experiments.  Linear 
regressions  through  these  data  showed  that  the  slopes  of  Rtox(S)  and  Rtox(G2  +  M)  were 
significantly  different  from  zero  in  both  experiments  (P  <  0.005  for  all  four  curves).  The 
slopes  of  Rtox(Gi  )  were  a=-0.01  for  the  temperature  experiment  and  a=0.05  for  the  range 
of  phosphate  limitations.  Both  slopes  were  statistically  not  different  from  zero  (P  >  0.97, 
C.I  95%  and  P  >  0.9,  C.I.  93%  respectively). 

Variations  in  Toxin  Composition.  The  changes  in  concentrations  (fmoles  per  cell)  and 
relative  proportions  (mole  percent  of  total  toxin)  of  the  eight  main  saxitoxin  derivatives 
(STX,  NEO,  GTX  I,IV,  GTX  II,III,  C  I,II)  were  determined  for  the  different  growth  rates 
in  both  experiments.  Variations  in  concentration  of  STX  derivatives  within  each 
experiment  were  proportional  to  the  changes  in  total  toxin  concentration,  resulting  in 
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Fig.  5.  Overall  net  toxin  production  rate,  Rtox  (dotted  line),  and  cell  cycle  stage 
specific  toxin  production  rates,  Rtox(Gi)  ( — ▼ — ),  Rtox(G2  +  M)  ( — □ — ) 
and  Rtox(S)  ( — • — ),  for  Alexandrium  fundyense  in  semi-continuous 
culture  under  (A)  temperature-limited  conditions  in  f/2  medium  or  (B) 
phosphate-limited  conditions  (Note  scale  changes  on  ordinates). 
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Fig.  6.  Toxin  composition  expressed  in  fmoles  per  cell  of  the  different  toxin 
derivatives  present  in  Alexandrium  fundyense  cells  in  semi-continuous 
culture  and  under  temperature  limited  conditions  in  f/2  medium  (A)  or  under 
phosphate  limited  conditions  (B).  (DC  I, II,  ■  GTX  I, IV,  A  GTX  II,III,  A 
NEO,  O  STX). 
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Fig.  7.  Toxin  composition  expressed  in  mole  per  cent  of  total  toxin  present  in 
Alexandrium  fundyense  cells  in  semi-continuous  culture  and  under 
temperature-limited  conditions  in  f/2  medium  (■)  or  under  phosphate-limited 
conditions  (□).  The  relative  proportions  of  the  different  toxins  fall  within 
three  distinct  groups  by  either  (a)  remaining  constant,  (b)  increasing  or  (c) 
decreasing  under  phosphate  limitation  (Error  bars  represent  SD;  n=7  for  the 
temperature  experiment,  n=5  for  the  phosphate  treatments). 
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constant  relative  proportions  of  the  toxins  for  each  of  the  two  experiments  (Fig.  5). 
Differences  were  only  observed  in  the  absolute  concentrations  and  mole%  of  the  different 
derivatives  between  treatments.  The  concentrations  of  C  I, II  and  GTX  II, III  were  higher 
by  at  least  an  order  of  magnitude  in  cells  grown  under  phosphate  starvation  compared  to 
cells  at  similar  grovidh  rates  in  replete  medium  (Fig.  6).  At  similar  growth  rates,  no 
differences  in  concentration  of  the  other  toxins  were  observed  between  both  experiments 
(Fig.  6).  Analysis  of  the  relative  proportions  provided  a  different  picture.  The  percentages 
of  C  I,II,  GTX  I,IV  and  the  traces  of  other  derivatives  (Others,  Fig.  7)  did  not  vary 
significantly.  The  relative  proportions  of  GTX  II,III  increased  significantly  under 
phosphate  stress  (12%±3  48%±4)  while  NEO  and  STX  decreased  from  27%+2  and 

21%+1  in  the  temperature  experiment  to  6%±3  and  8%±2  respectively  under  phosphate 
limitation  (Fig.  7). 


DISCUSSION 

In  the  past,  toxin  production  in  Alexandrium  fundyense  has  been  correlated  to  several 
environmental  factors  including  nutrient  availability,  salinity  and  temperature  (Proctor  et 
ai,  1975;  White,  1978;  Hall,  1982;  Boyer  et  al,  1987;  Ogata  et  al,  1987;  Anderson  et  al, 
1990b;  Matsuda  et  al,  1996).  Only  recently  has  the  correlation  of  toxin  production  with  a 
specific  physiological  process,  the  cell  cycle,  been  described  (Taroncher-Oldenburg  et  al, 
1997).  Those  studies  showed  that  in  synchronized  batch  cultures  STX  biosynthesis  is 
turned  on  only  during  the  Gi  phase  of  the  cell  cycle.  In  this  manuscript  we  describe  a 
series  of  independent  experiments  showing  that  toxin  content  is  directly  correlated  to  the 
length  of  the  cell  cycle,  and  more  precisely,  to  the  duration  of  the  Gi  phase  of  the  cell 
cycle.  Semi  continuous  cultures  of  Alexandrium  fundyense  were  utilized,  as  they  allowed 
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the  effects  of  growth  limitation  on  toxin  production  dynamics  to  be  analyzed  while  the 
extracellular  conditions  to  which  the  cells  had  adapted  remained  stable.  In  continuous 
cultures,  as  well  as  semi-continuous  cultures,  the  growth  rate  of  an  organism  is  determined 
by  the  rate  of  supply  of  the  limiting  nutrient  (Monod,  1950;  Fogg,  1965). 

In  one  case,  we  examined  the  extent  to  which  temperature-induced  changes  in  growth 
affected  toxin  production  and  its  linkage  to  cell  cycle  events  under  nutrient-replete 
conditions.  In  contrast,  phosphate-limited  treatments  were  established  to  study  the  effects 
of  nutrient-limited  growth  on  cell  proliferation,  cell  cycle  dynamics  and  toxin  production  at 
constant  temperature.  Phosphate  limitation  is  the  only  known  instance,  besides  low 
temperature,  in  which  saxitoxin  biosynthesis  is  enhanced  compared  to  non-limiting 
conditions  (Hall,  1982;  Boyer  era/.,  1987;  Anderson,  1990). 

The  influence  of  phosphate  availability  and  temperature  on  toxin  and  protein  production 
together  with  cell  cycle  dynamics  and  their  relationships  are  discussed  below. 

The  cell  cycle 

Over  the  range  of  growth  rates  studied,  we  observed  differing  effects  of  temperature 
limitation  and  phosphate  starvation  on  the  lengths  of  the  Gi  and  G2  phases  of  the  cell 
cycle.  As  was  observed  for  the  dinoflagellate  Amphidinium  carterii  over  a  range  of 
decreasing  temperatures  (Olson  et  al,  1986),  the  length  of  both  Gi  and  G2  increased  with 
decreasing  growth  rate  in  A.  fundyense  (Fig.  1  A).  A  similar  pattern  was  observed  under 
phosphate  limitation.  In  this  case,  the  length  of  Gi  was  greater  and  the  G2  phase  was 
shorter  at  all  growth  rates  compared  to  temperature-limited  growth  (Fig.  1 A  &  B).  Similar 
patterns  of  proportional  or  general  expansion  of  all  cell  cycle  phases  under  growth  limiting 
conditions  have  been  observed  in  mammalian  systems  and  yeast  (Rivin  and  Fangman, 
1980;  Guiget  et  al,  1984).  The  mechanistic  explanation  is  that  both  energy  and  nutrient 
limitation  act  on  processes  that  control  cell  cycle  progression  (e.g.,  DNA  replication  and 
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protein  biosynthesis).  Thus,  temperature  limitation  causes  a  general  decrease  in  enzyme 
activity,  resulting  in  lower  growth  rates  and,  in  the  case  of  A.  fundyense,  a  non¬ 
proportional  expansion  of  the  different  cell  cycle  stages. 

Under  phosphate  limitation  we  observed  not  only  a  longer  duration  of  the  Gi  phase, 
but  also  a  disproportional  expansion  of  the  different  phases.  The  most  plausible 
explanation  for  these  patterns  is  that  phosphate  is  an  essential  building  block  of  DNA  and 
hence  is  a  limiting  factor  for  cell  cycle  progression  when  present  at  low  concentrations. 
This  suggests  that  the  mechanism  controlling  cell  cycle  length  under  phosphate  limiting 
conditions  is  a  delay  of  the  transition  from  Gi  to  the  S  phase  as  the  cell  needs  additional 
time  to  accumulate  enough  phosphate  to  ensure  a  successful  DNA  replication  process 
(Pardee  et  al,  1978).  At  lower  growth  rates  induced  by  more  severe  phosphate  limitation, 
a  second  mechanism  might  become  operative.  Phosphate  is  not  only  an  essential  element 
for  DNA  synthesis,  but  also  is  a  necessary  component  of  energy  transfer,  via  ATP  and 
related  compounds,  and  is  critical  to  the  regulation  of  protein  activity  by  means  of  enzyme 
phosphorylation.  At  extremely  low  phosphate  concentrations,  protein  metabolism  in 
preparation  for  mitosis  would  be  restricted,  lengthening  G2.  Low  protein  levels  at  slow 
growth  rates  corroborate  this  interpretation  (Fig.  2B).  Temperature  limitation  and 
phosphate  limitation  thus  affect  the  cell  cycle  dynamics  in  different  ways  by  modulating 
separate  regulatory  machineries. 

Toxin  production  and  nitrogen  pools 

Saxitoxin  biosynthesis  is  dependent  on  the  availability  of  the  amino  acid  arginine  (Arg), 
its  main  building  block  (Shimizu,  1993).  Arg  is  one  of  the  principal  constituents  of 
proteins  and  its  metabolism  is  regulated  by  temperature  and  nutrient  levels,  among  other 
factors,  especially  nitrogen  and  phosphate  (Meister,  1965).  Similarly,  STX  biosynthesis  is 
dependent  on  factors  such  as  the  nitrogen  status  of  the  cells  and  phosphate  availability  in 
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the  medium  (Hall,  1982;  Boyer  et  al,  1987;  Anderson  et  ai,  1990b;  Matsuda  et  al,  1996). 
It  thus  became  of  interest  to  look  at  the  correlation  between  toxin  levels  and  protein 
accumulation  under  the  different  experimental  treatments  in  this  study.  If  toxin  and  protein 
production  were  directly  correlated,  the  ratio  of  the  two  components  would  remain  constant 
under  all  physiological  conditions.  We  saw  such  a  pattern  in  the  temperature  experiment, 
but  not  with  phosphate  limitation,  where  toxin  concentrations  became  progressively 
uncoupled  from  protein  concentrations  in  the  cell  as  the  limitation  increased  (Fig.  3).  This 
ratio  distribution  suggests  selective  maintenance  of  enzymes  of  toxin  biosynthesis  at  low 
phosphate  concentrations.  It  can  thus  be  hypothesized,  that  toxin  production  is 
advantageous  to  Alexandrium  as  it  becomes  progressively  phosphate-limited. 

Previous  work  has  shown  that  during  exponential  growth  of  Alexandrium,  pools  of 
Arg  do  not  show  any  significant  differences  between  phosphate-limited  and  nutrient-replete 
cultures  (Anderson  et  al,  1990b).  One  can  surmise  that  protein  and  toxin  production  could 
be  inversely  related  as  a  result  of  protein  catabolism  and  subsequent  channeling  of  Arg 
through  the  general  Arg  pool  into  STX  anabolism.  Such  a  scenario  would  be  supported  by 
the  fact  that  cellular  levels  of  Arg  as  well  as  those  of  other  amino  acids  are  under  tight 
control  in  order  to  preclude  high  accumulations  detrimental  to  the  cell's  physiology 
(Meister,  1965).  However,  analysis  of  the  variations  in  toxin,  protein  and  Arg  content  at 
low  temperatures  reveals  a  direct  correlation  between  these  factors  (Anderson  et  al , 
1990b)(Fig.  3),  suggesting  that  the  regulatory  pathway  of  STX  biosynthesis  does  not 
represent  a  straight  recycling  of  protein-derived  Arg  into  an  alternative  storage  compound, 
the  toxin.  Accumulations  of  protein  and  STX  at  low  temperatures  in  this  study  show  the 
same  direct  proportionality  (Figs.  2A  &  B).  We  conclude  that  toxin  production  is  not  a 
direct  response  to  a  cell's  need  to  balance  its  nitrogen  pools,  but  the  result  of  a  specific 
pathway  integrated  in  the  general  nitrogen  metabolism  of  the  cell. 
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Toxin  production  and  the  cell  cycle 

Our  previous  work  with  synchronized  cultures  in  nutrient-replete  medium  has  shown 
that  SIX  biosynthesis  is  turned  on  during  the  Gi  phase  of  the  cell  cycle  (Taroncher- 
Oldenburg  et  al,  1997).  The  main  hypothesis  we  wanted  to  test  with  the  experiments 
reported  here  was  whether  toxin  production  and  the  levels  of  toxin  accumulation  were 
related  to  one  specific  cell  cycle  stage  or  another  under  growth  limiting  conditions.  We 
first  set  out  to  determine  if  levels  of  STX  accumulation,  duration  of  the  different  cell  cycle 
stages  and  the  deduced  toxin  accumulation  rates  for  each  stage  were  correlated.  We  first 
compared  the  relative  increases  in  toxin  content  and  cell  cycle  phase  duration  for  both 
experiments  (Fig.  4A  &  B).  A  clear  overlap  was  observed  between  the  trends  and  values 
of  relative  durations  of  Gi  and  the  normalized  toxin  concentrations.  These  results  clearly 
indicated  that  toxin  production  was  related  to  the  length  of  the  Gi  cell  cycle  phase.  We 
then  went  on  to  determine  the  cell  cycle  stage-specific,  as  well  as  overall  net  toxin 
production  rates.  A  comparison  of  these  parameters  showed  similar  patterns  and  values  for 
the  net  and  the  Gi -specific  production  rates  (Fig.  5 A  &  B).  The  other  two  cell  cycle- 
specific  rates  (G2  +  M  and  S)  showed  not  only  a  high  degree  of  uncoupling  from  the 
overall  rates,  but  their  trends  also  contradicted  basic  concepts  in  metabolic  regulation. 
Typically,  the  end  product  of  a  biosynthetic  pathway  accumulates  as  a  function  of 
enzymatic  rates  and  the  length  of  the  time  interval  during  which  synthesis  occurs.  The  two 
extreme  cases  result  in  two  different  scenarios.  In  the  first  case,  high  product 
concentrations  are  derived  from  increases  in  production  rates  as  a  result  of  biochemical 
induction  of  the  responsible  enzymes.  Alternatively,  enzymatic  rates  can  remain  constant 
and  accumulation  of  the  end  product  be  the  consequence  of  longer  periods  of  biosynthetic 
activity.  This  appears  to  be  the  case  with  STX  production,  as  shown  by  our  results  (Fig. 
5).  In  both  the  temperature  treatments  and  the  phosphate  limitation  experiments,  toxin 
production  rates  remained  constant  for  a  range  of  conditions,  implying  that  the  variations  in 
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toxin  content  observed  (Fig.  2A)  were  a  result  of  increasing  periods  of  biosynthetic 
activity.  This  increase  in  the  length  of  the  toxin  production  was  directly  proportional  to  the 
changes  in  duration  of  the  Gi  phase  of  the  cell  cycle. 

The  G2  and  S  specific  toxin  production  rates  fall  into  a  third  category  in  which  an 
increase  in  metabolite  concentration  would  be  related  to  an  actual  decrease  in  production 
rates.  Such  dynamics  are  to  be  expected  only  in  the  case  of  a  storage  product  where 
regulation  of  biosynthesis  is  not  coupled  to  metabolite  accumulation.  One  of  the  possible 
roles  postulated  in  the  past  for  saxitoxin  is  as  a  storage  compound  for  nitrogen  (Loeblich, 
1984;  Boyer  et  al,  1987;  Anderson  et  al,  1990b).  Our  data  show  a  coincidence  of 
decreasing  protein  concentrations  with  an  increase  in  toxin  levels  under  phosphate 
limitation  (Fig.  2),  which  would  support  such  a  hypothesis,  considering  that  nitrogen  was 
not  a  limiting  factor  in  those  cultures  and  protein  anabolism  was  probably  compromised  by 
the  lack  of  phosphate.  However,  work  by  other  authors  showing  no  evidence  for 
reallocation  of  toxin-nitrogen  into  protein-nitrogen  (Hall,  1982;  Boyer  et  al,  1985;  Usup, 
1995;  Flynn  et  al,  1996),  together  with  our  observation  of  a  parallel  increase  of  toxin 
content  and  protein  levels  under  the  different  temperature  limitation  regimes  (Fig.  2A  &  B), 
contradicts  such  a  theory. 

Toxin  production  and  phosphate  availability 

The  role  of  phosphate  in  a  cell's  metabolism  is  a  complex  one.  The  physiological 
changes  we  observe  under  different  phosphate  starvation  regimes  are  the  result  of  the 
cumulative  effects  of  phosphate  limitation  on  the  entire  range  of  physiological  processes  in 
the  cell.  We  observed  that  both  the  net  and  the  Gi -specific  toxin  production  rates  increased 
by  a  factor  of  approximately  four  in  all  phosphate-limited  regimes  compared  to  the  nutrient 
replete  ones  over  the  same  range  of  growth  rates  (Fig.  5  A  &  B).  This  confirmed  previous 
toxin  data  comparisons  for  phosphate-limited  cultures  of  Alexandrium  where  toxin  levels 
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were  consistently  higher,  also  by  a  factor  of  four  (Boyer  et  ai,  1987;  Anderson  et  ai, 
1990b).  This  enhancement  of  toxin  production  suggests  that  phosphate  starvation  is  not 
acting  directly  on  the  biosynthetic  pathway  but  rather  increases  the  rate  of  substrate  supply 
via  down-regulation  of  other  cellular  processes  (Anderson  et  al,  1990b).  Moreover,  the 
rate  of  toxin  production  remains  at  a  constant  level  throughout  the  range  of  phosphate 
limitation  regimes  (Fig.  5B).  This  suggests  that  the  rate  of  approx.  3.3  fmol  celF^  houf^ 
is  the  maximum  rate  of  toxin  production  under  phosphate  stress  for  Alexandrium 
fundyense. 

It  is  of  interest  that  phosphate  deprivation  does  not  limit  toxin  production,  i.e.,  no 
reduction  in  the  toxin  production  rate  is  observed  (Fig.  5A  &  B),  which  implies  that 
phosphate  is  not  an  essential  cofactor  or  substrate  in  saxitoxin  biosynthesis.  This  is  further 
corroborated  by  the  observation  that  toxin  production  is  uncoupled  from  protein 
accumulation  under  phosphate  starvation  (Fig.  3).  Protein  metabolism  is  dependent  on 
phosphate  availability.  If  STX  biosynthesis  were  phosphate  dependent  we  should  observe 
a  similar  pattern  of  variation  as  the  one  recorded  for  protein  metabolism.  Instead,  opposite 
trends  were  followed  (Fig.  2B). 

We  conclude  that  STX  biosynthesis  is  not  directly  regulated  by  phosphate,  but  is 
indirectly  modulated  by  the  effects  of  phosphate  starvation  on  other  cellular  processes, 
such  as  protein  and  amino  acid  metabolism. 

Toxin  production  rate  and  growth  rate 

The  results  reported  here  show  that  net  toxin  production  rate  is  independent  of  growth 
rate  under  the  conditions  studied.  Net  toxin  production  rates  had  been  previously 
determined  under  phosphate  limitation  in  A.  fundyense  (Anderson,  1990).  The  average 
toxin  production  rate  was  of  approximately  75  fmol  cell'*  day'^  (=  3.12  fmol  cell'^  hour^). 
This  value  coincides  with  our  value  of  3.34  fmol  cell'^  hour  ^  In  that  instance  the  toxin 
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production  rates  were  suggested  to  be  correlated  to  growth  rate  based  on  their  increase  at 
faster  growth  rates.  The  discrepancy  between  those  data  and  the  results  reported  here  may 
lie  in  the  high  variability  seen  in  the  first  case.  This  is  further  substantiated  by  the  fact  that 
we  observed  constant  toxin  production  rates  in  nutrient-replete  cultures  and  under  varying 
temperatures.  Toxin  production  rates  have  also  been  shown  to  be  constant  over  a  range  of 
nitrogen  limiting  conditions  (Anderson  et  ai,  1990b).  These  data  suggest  that  once  a 
specific  limitation  threshold  is  reached,  toxin  production  rates  are  not  affected  by  further 
limitation.  This  has  important  implications  regarding  the  production  of  STX  as  it  suggests 
possible  regulatory  mechanisms  of  the  biosynthetic  pathway.  A  more  precise  biochemical 
interpretation  of  the  data  is  not  possible  at  this  point  because  none  of  the  enzymes  involved 
in  the  pathway  has  been  isolated  yet.  A  concerted  effort  in  determining  the  exact  sequence 
of  intermediates  and  the  nature  of  the  enzymes  responsible  for  their  processing  is  necessary 
in  order  to  overcome  this  last  hurdle  of  understanding  the  regulatory  mechanisms  involved 
in  the  biosynthesis  of  STX. 

Toxin  composition 

Toxin  composition  in  Alexandrium  had  long  been  considered  a  stable  feature  of  the 
different  species  or  strains  throughout  a  broad  range  of  physico-chemical  conditions  (Hall, 
1982;  Boyer  et  ai,  1987;  Ogata  et  al,  1987)  and  therefore  used  as  a  chemotaxonomic 
character  (Cembella  et  ai,  1987).  Only  more  recently  have  changes  in  toxin  composition 
over  the  life  cycle  or  under  varying  physiological  conditions  been  documented  (Boczar  et 
al,  1988;  Anderson,  1990;  MacIntyre  et  al,  1997).  Our  results  add  one  more  example  in 
which  a  significant  change  in  toxin  composition  was  observed  when  Alexandrium  was 
subjected  to  phosphate-limited  conditions  and  grown  in  semi-continuous  culture  in  nutrient 
replete  medium  (Fig.  6).  The  most  significant  changes  occured  at  the  level  of  absolute 
accumulation  of  the  carbamate  derivatives  GTX  II,III  and  the  sulfocarbamoylated  pair  C 
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I, II  (Fig.  6).  No  significant  differences  were  observed  for  the  total  concentrations  of  the 
two  other  carbamate  derivatives,  GTX  I,IV,  NEO  or  SIX  (Fig.  6).  Transformations 
among  STX  derivatives,  in  particular  ll-oc-P-O-sulfatation  of  STX  and  16-iV-sulfatation  of 
the  resulting  pair  GTX  II,III  are  enhanced  by  phosphate  starvation,  as  previously  shown 
by  Anderson  et  al.  (1990a).  In  a  similar  scenario  to  the  one  discussed  in  the  context  of 
enhancement  of  total  toxin  production  under  phosphate  limiting  conditions,  the  apparent 
increase  in  sulfatation  activity  is  most  probably  the  result  of  an  increase  in  substrate 
availability  for  these  enzymes,  rather  than  an  induction  of  enzyme  activity.  The 
characterization  of  the  two  enzymes  responsible  for  these  reactions  in  Gymnodinium 
catenatum  showed  that  their  activities  were  not  dependent  on  the  presence  of  phosphate 
(Yoshida  et  al,  1996).  We  can  further  explain  the  lack  of  a  parallel  and  significant  change 
in  the  total  concentration  of  the  pair  GTX  I,IV,  via  the  equivalent  transformation  of  the  1- 
W-hydroxy  STX  derivative  series,  by  the  fact  that  an  1 1-a-p-G-sulfatation  activity  of  NEO 
has  not  yet  been  isolated,  and  that  the  two  sulfatases  of  G.  catenatum  do  not  recognize 
NEO  or  the  pair  GTX  I,IV  as  substrates  (Yoshida  et  al,  1996).  A  plausible  explanation 
for  the  observed  patterns  would  therefore  be  the  presence  of  specific  l-A/-hydroxylases  or 
oxidases  of  the  kind  observed  before  in  Alexandrium  tamarense  (Oshima,  1995)  that  would 
in  turn  transform  STX  and  GTX  II,III  into  NEO  and  GTX  I,IV  respectively.  Such 
enzymes  are  always  ATP  dependent  and  hence  inhibited  in  the  absence  of  phosphate. 

The  resulting  scheme  of  transformations  is  outlined  in  Figure  8.  It  is  necessary  to  point 
out  that  the  arrows  follow  opposite  directions  from  those  we  had  postulated  in  previous 
experiments  with  synchronous  cultures  (Taroncher-Oldenburg  et  al,  1997),  but  a 
comparison  of  both  sets  of  observations  after  reinterpretation  of  the  toxin  accumulation 
patterns  in  those  synchronous  cultures  results  in  convergence.  In  that  instance,  we 
observed  first  an  accumulation  of  C  I,II  followed  by  GTX  II,in  and  STX  and  finally  GTX 
I,IV  and  NEO  with  a  slight  delay.  We  interpreted  this  to  be  the  sequence  in  which  these 
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compounds  were  formed  (i.e.,  that  C  I, II  were  the  'parent'  compounds).  These  data  can 
now  be  reinterpreted  in  light  of  our  current  results  as  a  process  during  which  C  I, II 
accumulates  as  a  result  of  the  rapid  transformation  of  its  direct  precursor,  GTX II, III,  until 
its  equilibrium  concentration  is  reached.  This  then  results  in  the  subsequent  accumulation 
of  GTX  II,III  until  its  concentration  reaches  a  maximum  and  STX,  the  precursor  of  this 
latter  transformation  starts  building  up.  The  hydroxylases/oxidases  involved  in  the 
transformation  of  STX  and  GTX  II,III  into  NEO  and  GTX  I,IV  only  start  showing  higher 
activities  once  the  sulfatation  reactions  that  are  competing  for  their  substrates  stabilize, 
suggesting  that  the  Km  of  the  sulfatases  is  lower  than  that  of  the  hydroxylases/oxidases. 

One  last  comment  addresses  the  divergent  results  obtained  when  toxin  composition  is 
analyzed  as  absolute  values  of  toxin  derivative  cell  quota  or  as  derivative  mole  percent  of 
total  toxin  in  intraspecific  comparisons  between  sets  of  toxin  data  where  the  total  toxin 
content  has  varied  significantly  (Figs.  6  &  7).  In  such  cases  the  relative  amounts  of  the 
different  derivatives  have  to  be  contrasted  with  their  actual  concentrations  in  order  not  to 
distort  the  real  changes  in  composition  occurring  in  the  cells.  This  is  illustrated  by  our 
results;  In  the  temperature  experiment,  all  derivatives  are  loosely  grouped  at  similar  levels 
of  concentration,  C  I,II  consistently  showing  the  highest  concentration  (Fig.  6A).  In 
contrast,  under  phosphate-limiting  conditions  GTX  II,III  are  very  clearly  present  at  the 
highest  concentration,  with  C  I,II  coming  second  and  the  rest  of  the  derivatives  all 
clustering  at  similar  concentrations  to  those  observed  in  the  temperature  treatments  (Fig. 
6B).  These  results  indicate  a  pattern  of  enhanced  accumulation  of  GTX  II,III  and  C  1,11 
under  phosphate  limitation.  Analysis  of  these  same  data  by  normalizing  the  concentration 
of  each  derivative  to  total  toxin  provides  a  different  picture  (Fig.  7).  A  clear  change  can 
still  be  seen  for  GTX  II,III  but  the  proportion  of  C  I,II  remains  constant.  In  addition, 
NEO  and  STX  also  show  a  significant  decrease.  These  differences  between  analytical 
approaches  reflect  the  loss  of  information  associated  with  representing  toxin  concentrations 
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Fig.  8.  Suggested  scheme  of  interconversions  among  the  main  STX  derivatives 
synthesized  by  Alexandrium  fundyense.  The  arrows  show  the  most 
plausible  direction  for  the  transformations.  ©  denotes  possible  inhibition  of 
specific  enzyme  activity. 
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as  relative  amounts  of  the  total  toxin.  While  not  invalidating  mole  %  representations  as  a 
means  for  comparing  different  organisms  or  determining  qualitative,  and  to  a  certain  extent, 
quantitative  variations  for  one  species,  our  results  imply  that  changes  in  the  physiology  of 
toxin  interconversions  can  only  be  inferred  from  the  absolute  values  of  the  toxin  cell  quota 
for  each  derivative. 

Summarizing  the  toxin  composition  portion  of  this  study,  the  patterns  of  change  we 
observed,  together  with  the  fact  that  within  treatments  the  variation  in  toxin  composition 
among  the  different  growth  rates  was  minimal,  suggest  that  phosphate  starvation  affects 
toxin  speciation  at  the  level  of  hydroxylase/oxidase  activities.  In  an  integrated  scheme  of 
the  STX  biosynthetic  pathway,  interconversions  among  the  toxins  would  be  determined  in 
part  by  the  presence  or  absence  of  the  appropriate  enzymes,  but  net  total  toxin  production 
would  be  the  result  of  the  accumulation  of  only  one  precursor  or  parent  compound.  This 
precursor  could  be  one  of  the  toxins  or  another  unknown  compound  (Shimizu,  1993). 
Further  efforts  to  identify  and  isolate  intermediate  compounds  are  necessary  in  order  to 
determine  the  exact  sequence  of  the  biosynthetic  pathway. 
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downregulated  genes  during  saxitoxin  biosynthesis  in 
Alexandrium  fundyense  (Dinophyceae) 
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ABSTRACT 


Genes  showing  a  pattern  of  transcriptional  induction  or  suppresion  that  coincided 
with  changes  in  saxitoxin  (STX)  biosynthesis  in  the  toxic  dinoflagellate  Alexandrium 
were  identified  and  characterized.  First,  differential  display  (DD)  of  mRNA  was  applied 
to  the  comparative  analysis  of  different  toxic  and  non-toxic  strains  of  Alexandrium 
tamarense  and  A.  fundyense.  The  genetic  variation  was  large  enough  to  make  it 
impossible  to  identify  differentially  expressed  genes  over  the  general  variation  between 
DD  patterns  due  to  inter-strain  dissimilarities.  Interspecific  variations  between  A. 
tamarense  and  closely  related  strains  of  A.  fundyense  were  equally  uninformative.  A 
second  approach  consisted  in  DD  analysis  of  cells  in  different  cell  cycle  stages  and  while 
toxin  production  was  turned  on  or  off  during  synchronized  growth  of  A.  fundyense. 
Three  genes  were  isolated:  5'-adenosylhomocysteine  hydrolase  (SAHH),  methionine 
aminopeptidase  (MAP)  and  a  histone-like  protein  (HAf).  SAHH  was  downregulated, 
while  MAP  and  HAf  were  upregulated  during  the  interval  when  toxin  biosynthesis  was 
occurring.  SAHH  and  MAP  showed,  respectively,  about  90%  and  70%  similarity  at  the 
aminoacid  level  with  several  eukaryotic  and  prokaryotic  SAHHs  and  MAPs.  The  partial 
MAP  sequence  also  contained  three  cobalt  binding  motifs  characteristic  of  MAPs.  HAf 
showed  a  60%  similarity  at  the  aminoacid  level  with  two  histone-like  proteins  from  the 
dinoflagellate  Crypthecodinium  cohnii  Biecheler.  None  of  the  three  genes  isolated 
seemed  to  specifically  relate  to  STX  biosynthesis.  One  gene,  encoding  the  putative 
histone-like  protein  HAf,  could  be  a  cell  cycle-related  gene  since  it  encodes  for  a  protein 
that  is  surmised  to  have  DNA  structural  and  regulatory  functions.  The  two  other  genes, 
SAHH  and  MAP,  are  both  involved  in  regulatory  pathways.  As  an  active  component  in 
the  modulation  of  one-carbon  metabolism  and  the  biosynthesis  of  S-adenosylmethionine, 
a  precursor  of  STX,  SAHH  is  the  enzyme  most  closely  related  to  STX  biosynthesis.  This 
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study  further  documents  the  potential  of  applying  DD  to  the  identification  of  genes  that 
are  related  to  physiological  processes  or  cell  cycle  events  in  phytoplankton  under 
conditions  where  small  sample  volumes  represent  an  experimental  constraint.  We  also 
show  for  the  first  time  that  pre-translational  or  transcriptional  regulation  does  indeed 
occur  in  dinoflagellates.  Previous  reports  had  suggested  the  possibility  that  translational 
mechanisms  could  be  the  primary  means  of  circadian  regulation  in  these  organisms.  We 
have  shown  not  only  the  induction  or  suppression  of  the  three  genes  described  above,  but 
also  observed  transcriptional  or  pre-translational  regulation  of  about  twenty  other  genes 
with  varying  DD  patterns. 
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INTRODUCTION 


Marine  dinoflagellates  represent  a  substantial  fraction  of  bloom-forming 
phytoplankton  in  coastal  areas  of  the  world's  oceans  (Anderson,  1989).  In  addition  to  the 
formation  of  red  tides,  some  species  of  dinoflagellates  also  produce  a  range  of  toxins  that 
are  poisonous  to  higher  levels  of  the  food  web  (Shimizu,  1993).  One  group  of  toxins 
synthesized  by  these  harmful  algal  bloom  organisms  causes  Paralytic  Shellfish  Poisoning 
or  PSP.  The  toxins  included  in  this  group  are  called  saxitoxins  and  are  all  cyclic 
perhydropurine  compounds  differing  slightly  in  toxicity  (Hall  et  al,  1990;  Shimizu, 
1993).  Besides  dinoflagellates,  other  organisms  also  biosynthesize  PSPs,  as  for  example 
bacteria  (Silva  and  Sousa,  1981;  Kodama  et  al,  1988)  and  cyanobacteria  (Shimizu  et  al, 
1984;  Negri  et  al,  1997). 

Much  work  has  been  done  on  the  ecology  and  physiology  of  STX  production  in 
Alexandrium  fundyense  (Dinophyceae)  as  well  as  on  determining  the  basic  precursors  of 
this  compound.  As  a  result,  we  now  know  that  saxitoxin  is  synthesized  from  one 
molecule  of  arginine  (Arg),  three  guanido/ureido  groups,  derived  from  Arg  as  well,  an 
acetate  moeity,  probably  via  Acetyl-CoA,  and  one  methyl  group  donated  by  S-adenosyl- 
methionine.  Further  studies  of  the  biosynthetic  pathway  have  been  hampered  by 
difficulties  in  manipulating  this  organism  at  the  biochemical  and  genetic  level. 
Conventional  tools  such  as  mutational  analysis  or  transformation  cannot  be  applied  to 
identify  the  actual  'toxin  genes'  because  of  biological  and  technical  constraints  such  as 
lack  of  information  on  the  number  of  copies  of  the  'toxin  genes',  limitations  in  culturing 
dinoflagellates  on  solid  media  and  the  high  content  of  DNA  per  cell  that  would  make  it 
necessary  to  screen  unrealistic  numbers  of  clones  in  order  to  detect  one  mutant. 

Classic  mating  experiments  with  toxic  and  non-toxic  strains  of  Gymnodinium  sp.  and 
Alexandrium  sp.  have  revealed  that  inheritance  of  toxigenic  potential  and  of  toxin 
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composition  in  dinoflagellates  follow  Mendelian  rules  (Sako  et  al.,  1992;  Oshima  et  ai, 
1993;  Sako  etai,  1995).  The  non-random  segregation  patterns  observed  in  the  crossings 
also  imply  chromosomal  localization  of  the  genes  involved  in  toxin  biosynthesis  and 
interconversions  among  derivatives  (Oshima,  1995;  Sako  et  al.,  1995).  A  series  of 
studies  with  synchronized  cultures  of  Alexandrium  fundyense  has  further  shown  that  STX 
is  synthesized  during  a  discrete  time  period  localized  in  the  Gi  phase  of  the  cell  cycle 
(Fig.  1)  (Taroncher-Oldenburg  et  al,  1997).  The  combination  of  these  two  factors, 
nuclear  genes  and  discontinuous  toxigenesis,  gives  rise  to  the  hypothesis  that  STX 
biosynthesis  could  be  regulated  at  the  transcriptional  level  with  the  genes  responsible  for 
toxin  production  being  activated  in  a  cyclic  pattern.  Furthermore,  comparative  analysis 
of  cells  harvested  when  toxin  production  is  on  or  off  would  allow  the  identification  of 
genes  that  are  expressed  while  STX  is  produced  and  are  therefore  potentially  related  to 
toxigenesis. 

Recent  advances  in  the  development  of  molecular  tools  aimed  at  analyzing 
differences  between  complex  genomes  at  the  level  of  gene  expression  show  great  promise 
in  this  regard.  One  such  tool  is  Differential  Display,  a  genetic  approach  that  allows  the 
analysis  of  differentially  expressed  genes  in  eukaryotes  (Liang  and  Pardee,  1992).  This 
technique  has  been  successfully  used  to  identify  developmentally,  environmentally  or 
hormonally  regulated  genes  in  protozoa,  plants,  animals  and  humans  (Utans  et  al.,  1994; 
Abrahamsen  et  ai,  1995;  Sharma  and  Davis,  1995;  van  der  Knaap  and  Kende,  1995).  In 
this  procedure  (Fig.  2),  the  two  or  more  sets  of  differentially  expressed  messenger  RNAs 
(mRNAs)  to  be  compared  are  used  as  templates  to  generate  the  corresponding 
complementary  copy  DNAs  (cDNAs).  Subsequently,  the  cDNA  fragments  are  amplified 
by  means  of  the  polymerase  chain  reaction  (PCR)  using  a  combination  of  an  mRNA 
specific  oligo(dT)  "anchor"  primer  together  with  an  arbitrary  decamer  primer  in  the 
presence  of  radioactively  labeled  dNTPs.  After  separation  in  denaturing  polyacrylamide 
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Figure  1  Saxitoxin  biosynthesis  and  the  cell  cycle  in  synchronously  growing 
Alexandrium  fundyense.  Lines  show  relative  distributions  of  cells  in  the 
different  cell  cycle  stages  (Gi,  S,  and  G2  +  M,  Positive  toxin 
production  rates  (black  bars)  are  circumscribed  to  the  first  eight  to  ten  hours 
of  the  duration  of  the  Gi  cell  cycle  phase.  No  toxin  production  is  detected  for 
the  remainder  of  the  cell  cycle.  Arrows  denote  sampling  points  for  the 
Differential  Display  analysis  protocol:  Two  during  Gi,  (I)  when  toxin  is 
accumulated  and  (II)  after  it  stops,  one  during  the  S  phase  (III)  and  a  last  one 
during  G2  (IV).  Shaded  areas  correspond  to  the  dark  periods  of  the  14:10 
hour  light:  dark  cycle  under  which  A.  fundyense  was  synchronized. 
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Figure  2  Differential  Display.  The  mRNA  components  to  be  compared  undergo 
successive  steps  of  reverse  trancription  (RT)  and  amplification  through  the 
polymerase  chain  reaction  (PCR)  using  on  mRNA  specific  'anchor'  primer 
and  one  'random'  decamer.  This  is  followed  by  visualization  of  the  PCR 
products  on  a  PAGE  gel.  Differences  in  expression  are  identified,  extracted 
and  cloned  for  further  analysis.  (See  materials  and  methods  for  details) 
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gel  electrophoresis,  differentially  amplified  cDNAs  are  identified  by  autoradiography. 
These  cDNAs  are  further  cloned  and  used  as  probes  to  screen  cDNA  libraries  of  the 
organism  of  interest. 

The  goal  of  this  study  was  to  identify  genes  for  which  expression  was  correlated  to 
toxin  production  in  Alexandrium.  Differential  Display  was  used  to  analyze  mRNA 
expression  patterns  among  cells  differing  in  their  toxigenic  status.  A  first  approach 
consisted  in  comparing  phylogenetically  closely  related  toxic  and  non-toxic  strains  of  A. 
tamarense  and  A.  fundyense.  A  second  approach  consisted  in  the  identification  of  genes 
for  which  expression  cycles  were  correlated  to  the  discontinuous  pattern  of  STX 
biosynthesis  in  A.  fundyense.  Synchronously  growing  cultures  of  A.  fundeyense  were 
harvested  several  times  throughout  their  cell  cycle  at  times  of  variable  toxigenic  status. 
Several  genes  were  identified  and  their  possible  role  in  STX  biosynthesis  explored. 


MATERIALS  AND  METHODS 

Organisms.  Batch  cultures  of  the  non-toxic  strains  of  Alexandrium  tamarense 
PGtl83  (Plymouth,  UK),  PelV  (Galicia,  Spain)  and  WKS-1  (Tanabe  Bay,  Japan)  and  the 
toxic  strains  of  A.  tamarense  GtPP03  (Perch  Pond,  MA),  AtSL12  (St.  Lawrence  estuary, 
Canada),  GtLI21  (Moriches  Bay,  NY)  and  A.  fundyense  GtCA28  (Cape  Ann,  MA)  and 
Gt7  (Bay  of  Fundy,  Canada)  were  used  in  this  study.  Cultures  were  maintained  at  20®C 
in  f/2  medium  (Guillard  and  Ryther,  1962)  modified  by  addition  of  H2Se03  (10‘^M)  and 
by  reducing  the  concentration  of  CUSO4  to  10‘®M.  Vineyard  Sound  seawater  (0.2  pm 
filtered,  31°/oo  salinity)  was  used  as  the  medium  base.  Throughout  the  experiment, 
irradiance  of  ca.  250  pE  m‘^s‘*  was  provided  by  cool  white  fluorescent  bulbs  on  a  14:10  h 
light:dark  cycle.  Cultures  were  harvested  in  mid-exponential  phase. 
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Synchronization  of  Alexandrium  fundyense.  Cultures  of  Alexandrium  fundyense 
GtCA28  (see  above  for  details)  were  synchronized  by  a  dark  induced  block/release 
method  (Taroncher-Oldenburg  et  ai,  1997).  Low-density  (approx.  500  cells/ml), 
exponentially  growing  cultures  (L:D  14:10h)  of  Alexandrium  fundyense,  strain  GtCA28, 
were  left  in  continous  darkness  for  82  hours  and  then  re-exposed  to  the  original 
Light:Dark  cycle.  At  the  end  of  the  dark  period,  all  the  cells  had  divided  once  more  and 
were  arrested  in  the  Gi  phase  of  their  cell  cycle.  Cell  size  and  density  were  monitored 
with  a  Multisizer  apparatus  (Coulter  Counter  Ltd.,  GB)  every  two  hours  for  a  period  of 
two  days.  Samples  for  DNA  analysis  (50,000-100,000  cells)  were  taken  with  the  same 
periodicity  and  preserved  in  5%  formaldehyde  at  -20'C.  The  samples  were  treated  with 
the  DNA-specific  stain  propidium  iodide  (PI)  (4  pg  ml'i  PI,  250  units  RNaseA  ml'O  and 
analyzed  on  a  FaxScan  Flowcytometer  (Becton  Dickinson). 

RNA  Purification.  1-3  x  10^  cells  were  harvested  through  centrifugation  at  4,400g  for 
5  minutes  at  4°C.  The  pellet  was  deep-frozen  in  liquid  N2  and  then  resuspended  in  1  mL 
of  RNA  STAT-60  (TEL-TEST  "B",  Inc.,  TX).  The  suspensions  were  placed  in  a  N2  cell 
disruption  bomb  and  the  cells  ruptured  in  three  consecutive  cycles  of  nitrogen 
decompression  at  2,000  psi  (Parr  Instrument  Company,  Moline,  IL).  Homogenates  were 
processed  following  the  manufacturer's  protocol  (RNA  precipitation  for  1  hour)  with  an 
additional  RNA  precipitation  in  the  presence  of  NaOAc  (pH  5.2,  final  cone.  0.3M)  and 
two  volumes  of  ethanol  for  1  hour  at  -70'C.  RNA  was  quantified  and  stored  at  -70‘C. 

Differential  Display.  Differential  Display  (DD)  of  mRNA  was  performed  as 
previously  described  by  Liang  and  Pardee  (Liang  and  Pardee,  1992;  Liang  et  al,  1993). 
The  reverse  transcription  and  PCR  were  modified  as  follows.  Total  RNA  (0.2|ig)  was 
reverse  transcribed  with  200  units  of  M-MLV  Reverse  transcriptase  (Superscript,  RNase 
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H",  GIBCO-BRL),  40)iM  dNTPs  and  l)iM  of  one  of  the  four  different  degenerate 
oligo(dT)  primers  (NBI,  Plymouth,  MN)  T12VN  (V=dA,dC  or  dG).  Control  reactions 
were  performed  with  total  RNA  in  the  absence  of  RTase  and  with  RNAse  treated  RNA  in 
the  presence  of  RTase.  The  resulting  cDNA  was  then  amplified  by  PCR  with  a 
combination  of  the  respective  degenerate  oligo(dT)  primer  and  one  random  decamer 
primer  (OPERON  10-mer  kit  A,  Operon  Technologies)  and  in  the  presence  of  a-[^^P]- 
dATP  on  a  Perkiii-Elmer  2400  thermocycler.  All  components  of  the  PCR  remained  the 
same  as  originally  described  except  for  a  ten-fold  increase  in  the  concentration  of  random 
decamer  primer  (2)iM).  The  samples  were  denatured  at  94'C  for  15  sec.  directly 
preceding  PCR.  Parameters  for  the  40-cycle  PCR  were  as  follows;  denaturation  at  94'C 
for  15  sec.,  annealing  at  40’C  for  1  min.,  and  extension  at  72'C  for  1  min..  At  the  end  an 
additional  extension  at  72'C  for  10  min.  was  performed.  The  radiolabeled  PCR-products 
were  analyzed  on  6%  denaturing  polyacrylamide  gels. 

Cloning  and  Sequencing.  Bands  identified  as  differentially  expressed  on  the  DD  gels 
were  excised  (details  in  RESULTS),  extracted  from  the  acrylamide  by  boiling  in  H2O  and 
precipitated  in  the  presence  of  0.3M  NaOAc,  0.5  mg/mL  glycogen  and  four  volumes  of 
ethanol.  The  DNA  was  then  reamplified  using  the  same  primer  and  PCR  conditions  as 
above,  only  the  dNTP  concentration  was  increased  10-fold.  The  PCR  product  was 
cloned  using  the  pCR-TRAP  cloning  system  (GenHunter  Corp.,  Brookline,  MA).  Cloned 
fragments  were  sequenced  on  a  LiCor  automated  sequencing  system  using  pCR-TRAP 
specific  primers  (Aidseq,  GenHunter  Corp.,  Brookline,  MA). 

cDNA  Library  Construction  and  Screening.  Approximately  10*  cells  of 
exponentially  growing  Alexandrium  fundyense  were  harvested  by  centrifugation  at 
4,400g  for  5  minutes  at  4'C.  The  pellet  was  deep-frozen  in  liquid  N2  and  then 
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resuspended  in  20  mL  of  RNA  STAT-60  (TEL-TEST  "B",  Inc.;  TX).  The  suspension 
was  placed  in  a  N2  cell  disruption  bomb  and  the  cells  ruptured  in  three  consecutive  cycles 
of  nitrogen  decompression  at  2,000  psi  (Parr  Instrument  Company,  Moline,  IL). 
Homogenates  were  processed  following  the  manufacturer's  protocol  (RNA  precipitation 
for  1  hour)  with  an  additional  RNA  precipitation  in  the  presence  of  NaOAc  (pH  5.2,  final 
cone.  0.3M)  and  two  volumes  of  ethanol  for  1  hour  at  -70°C.  Approximately  one  mg  of 
fresh  total  RNA  was  further  purified  for  poly(A)+  mRNA  with  an  oligo(dT)  column 
(5Prime->3Prime,  Inc.,  Boulder,  CO).  Spg  of  poly(A)+  mRNA  were  used  for 
construction  of  a  lambda  phage  cDNA  library  (ZAP-cDNA  cloning  kit,  Stratagene,  La 
Jolla,  CA).  The  library  was  further  amplified  (  2  x  10^  pfu/pl)  and  approximately  1  x  10^ 
plaques  were  screened  on  each  round. 

DD  probes  for  screening  the  library  were  obtained  by  labeling  the  cloned  cDNA  of 
interest  with  Digoxigenin-ll-dUTP  (Genius  Nonradioactive  Nucleic  Acid  Labeling  and 
Detection  System,  Boehringer  Mannheim,  IN)  in  a  30  cycle  PCR  with  the  following 
parameters:  denaturation  at  94°C  for  30  sec.,  annealing  at  52°C  for  40  sec.,  and  extension 
at  72°C  for  1  min..  At  the  end  an  additional  extension  at  72°C  for  5  min.  was  performed. 
The  final  concentration  of  dNTPs  in  the  reaction  mix  was  optimized  to  the  following 
values:  20  pM  dATP,  dCTP  and  dGTP  respectively,  13  pM  dTTP  and  7  pM  DIG-dUTP. 
The  yield  of  Dig-labeled  PCR  product  was  estimated  as  indicated  by  the  manufacturer 
and  hybridizations  were  conducted  overnight  in  a  standard  hybridization  buffer  (5  x  SSC, 
2%  (w/v)  blocking  reagent,  0.1%  N-lauroylsarcosine,  0.2%  SDS),  with  a  probe 
concentration  of  lOng/mL  and  at  a  temperature  of  65°C.  Positive  plaques  were  detected 
using  the  chemiluminescent  substrate  CSPD  (Boehringer  Mannheim,  IN).  Plaques  were 
further  purified  in  a  second  and  if  necessary  tertiary  screening,  until  single  plaques  could 
be  isolated.  The  phagemid  representing  the  cDNA  of  interest  was  then  excised  in  vivo. 
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cloned  and  sequenced  following  the  manufacturer's  protocols  (ZAP-cDNA  cloning  kit, 
Stratagene,  La  Jolla,  CA). 

RNA  Blot  Analysis.  Ten  pg  of  total,  DNAse  treated  RNA  representing  each  of  the 
four  cell  cycle  stages  described  above,  were  loaded  on  positively  charged  nylon 
membranes  (Boehringer  Mannheim,  IN).  The  DD  probes  were  labeled  as  described 
above.  Hybridizations  were  performed  0/N  in  high  SDS  buffer,  as  recommended  by  the 
manufacturer,  and  at  42°C  (Genius  Nonradioactive  Nucleic  Acid  Labeling  and  Detection 
System,  Boehringer  Mannheim,  IN).  Detection  of  positive  dots  was  performed  as 
detailed  above  for  the  cDNA  library  screenings. 

Sequence  Analysis.  The  nucleotide  sequences  obtained  were  compared  to  the  non- 
redundant  GenBank,  EMBL,  DDBJ  and  PDB  sequence  databases  using  the  Basic  Local 
Alignment  Search  Tool  program  (BLAST)  in  its  version  for  nucleotides  (BLASTN)  or 
aminoacids  (BLASTX)  (Altschul  et  al,  1990).  The  corresponding  amino  acid  sequence 
was  predicted  using  the  universal  genetic  code.  The  deduced  protein  sequences  were 
further  characterized  with  the  sequence  analysis  software  program  Mac  Vector  (Oxford 
Molecular  Limited). 


RESULTS 

Differential  Display  of  toxic  and  non-toxic  Alexandrium.  Differential  Display  of 
several  toxic  and  non-toxic  strains  of  Alexandrium  tamarense  and  two  toxic  A.  fundyense 
(see  MATERIALS  AND  METHODS  for  specific  strain  information)  was  done  using 
T12VC  as  the  degenerate  oligo(dT)  primer  and  OP-03  as  the  random  decamer.  High 
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variability  was  observed  among  the  DD  patterns  of  the  different  strains  of  A.  tamarense 
analyzed  (Fig.3).  Interspecific  variations  between  A.  tamarense  and  A.  fundyense  were 
also  high  (Fig.  3). 

Differential  Display  of  Synchronized  Cultures  of  Alexandrium  fundyense.  Differential 
Display  (DD)  of  mRNA  was  applied  to  synchronized  cells  of  Alexandrium  fundyense  in 
different  cell  cycle  stages  and  while  toxin  production  was  on  or  off  (Taroncher- 
Oldenburg  et  al,  1997).  Two  samples  were  drawn  during  Gi,  one  while  toxin  production 
was  active  (I)  and  a  second  one  once  toxigenesis  was  over  (II),  a  third  sample  was 
constituted  by  cells  in  S  phase  (III),  and  the  fourth  one  represented  cells  that  were 
predominantly  in  G2  (IV)  (Fig.  1).  A  total  of  80  different  anchor  plus  random  primer 
combinations  were  used  to  ensure  total  representation  of  the  mRNA  populations  (Liang  et 
al,  1993).  Twenty  four  bands  were  identified  as  differentially  expressed  at  the  four  time 
points  selected.  Several  different  expression  patterns  were  observed.  Most  bands  were 
present  throughout  the  cell  cycle  (i.e.  not  differentially  expressed),  as  exemplified  by 
band  D  on  Figure  4.  Others  showed  cell  cycle  stage  specific  expression  as  was  the  case 
of  band  E  on  the  same  figure,  which  was  induced  throughout  Gi  but  was  not  detectable 
during  S  and  G2.  Only  those  bands  corresponding  to  genes  showing  suppresion,  band  A, 
or  induction,  bands  B  and  C,  during  toxigenesis  at  the  beginning  of  Gi  were  isolated  for 
further  characterization  (Fig.  4).  Three  bands  were  singled  out  by  these  criteria,  excised, 
cloned  and  subsequently  used  as  probes  to  recover  the  corresponding  cDNAs  from  an 
Alexandrium  fundyense  cDNA  library.  The  DD  expression  patterns  were  also  confirmed 
by  RNA  dot  blot  analysis  (Fig.  5).  The  three  genes  were  similar  to  S- 
adenosylhomocysteine  hydrolase  (SAHH),  methionine  aminopeptidase  (MAP)  and  a 
histone-like  protein  (HAf)  respectively. 
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Figure  3.  Comparison  of  inter-  and  intraspecific  DD  analysis  of  dinoflagellates.  In  the 
left  panel,  three  toxic  strains  of  Alexandrium  tamarense,  three  non-toxic 
strains  of  A.  tamarense  and  two  toxic  strains  of  A.  fundyense  are  compared. 
The  inter-  and  intrageneric  variabilities  are  high.  In  contrast,  in  the  right 
panel,  three  different  cell  cycle  stages  in  synchronyzed  cultures  of  A. 
fundyense  (GtCA28)  are  analyzed.  The  patterns  are  identical  and  hence 
allow  the  identification  of  differentially  expressed  bands  (Fig.  4).  (The  left 
image  is  a  composite  of  lanes  from  one  gel  rearranged  for  clarity). 
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Figure  4.  Differential  Display  comparing  cells  of  Alexandrium  fundyense  in  three 
different  cell  cycle  stages  (Gi,  S  and  G2)  and  while  STX  production  is  turned 
ON  (+)  or  OFF  (-)  (Numerals  I  through  IV  refer  to  the  four  sampling  points 
in  Fig.  1;  double  lanes  represent  duplicate  RT-PCRs).  Band  A  is 
downregulated,  while  bands  B  and  C  are  upregulated  relative  to  toxin 
production.  Band  D  is  constitutively  expressed,  while  band  E  is  Gi  -specific 
(expressed  in  early  and  late  Gi). 
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Figure  5.  RNA  dot  blot  analysis  confirming  the  induction  of  DD  bands  A,  B,  and  C  in 
Alexandrium  fundyense.  Total  RNA  (lOjig  per  dot)  corresponding  to  the  four 
sampling  points  (Numerals  I  through  IV,  Fig.  1)  was  hybridized  with  cDNA 
probes  generated  by  PCR  of  cloned  cDNA  fragments  obtained  from  the 
differential  display  gels. 


Toxin  Productton 


ON  OFF 

f— I  i  . . . ’'""1 

I  II  in  IV 


141 


S-adenosylhomocysteine  hydrolase.  The  131  bp  band  A  was  downregulated  during 
toxin  biosynthesis  (Fig.  4).  Its  sequence  showed  that  only  the  random  decamer  had  been 
applied  during  the  PCR  amplification  step  of  DD  (Fig.  6).  It  has  been  shown  previously 
that  a  significant  proportion  of  DD  bands  can  result  from  the  amplification  of  cDNA  with 
only  one  primer,  either  the  anchored,  mRNA-specific  primer  or  the  random  decamer 
(Liang  et  al,  1993;  Guimaraes  et  al,  1997).  These  bands  are  equally  reproducible  and 
represent  potential  transcriptionally  regulated  genes.  Such  single  primer  bands  are  the 
result  of  the  random  decamer  sequence  and  its  inverse  complement  being  present  on  the 
same  cDNA.  Generation  of  these  bands  resembles  arbitrary  primed  PCR  fingerprinting 
of  RNA,  a  similar  technique  to  DD,  that  relies  on  the  use  of  only  one  primer  for  the 
generation  of  cDNAs  representative  of  expressed  sequences  in  eukaryotic  and  prokaryotic 
cells  (Welsh  et  al.,  1992;  Wong  and  McClelland,  1994). 

The  alignment  of  the  partial  protein,  as  deduced  from  the  rescued  cDNA  fragment, 
showed  a  90%  similarity  with  SAHHs  from  eukaryotes  and  prokaryotes  (tobacco 
(Nicotiana  sylvestris)  85%  similarity,  75%  identity;  yeast  {Schizosaccharomyces  pombe) 
75%  similarity,  62%  identity;  a  cyanobacteria  (Synechocystis  sp.)  63%  similarity,  44% 
identity)  (Fig.  7). 

Methionine  Aminopeptidase.  The  115  bp  band  B  was  upregulated  during  toxin 
biosynthesis  in  Gi  (Fig.  4).  The  fragment  rescued  from  the  cDNA  library  corresponded 
to  the  3'  end  of  the  mRNA  (Fig.  8).  The  alignment  of  the  deduced  protein  with  the 
carboxy  ends  of  three  other  MAPs  showed  a  similarity  of  70%  (yeast  {Saccharomyces 
pombe)  61%  similarity,  44%  identity;  a  cyanobacteria  (Synechocystis  sp.)  61%  similarity, 
42%  identity;  Escherichia  coli  63%  similarity,  46%  identity)  (Fig.  9).  In  addition,  three 
of  the  five  highly  conserved  copper  binding  regions  that  have  been  described  for  this 
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enzyme  were  identified  on  the  A.  fundyense  sequence  (Fig.  10;  shaded  areas).  The  other 
two  metal-binding  domains  fall  within  the  missing  section  of  the  sequence. 

Histone-like  Protein.  The  318  bp  band  C  was  upregulated  during  toxigenesis  and 
coinciding  with  the  beginning  of  the  Gi  phase  of  the  cell  cycle  (Fig.  4).  A  577  bp  cDNA 
with  a  255  bp  open  reading  frame  (ORF)  was  recovered  (Fig.  10).  The  database  search 
aligned  this  sequence  to  two  histone-like  proteins  from  Crypthecodinium  cohnii  Biecheler 
(Fig.  11).  The  combined  similarity  was  60%,  with  each  of  the  two  different  histone-like 
proteins  showing  a  similarity  of  52%  and  an  identity  of  35%.  The  ORF  obtained  from  the 
cDNA  library  frames  a  protein  with  85  amino  acids  and  a  molecular  weight  of 
approximately  9.3  kDa  (Fig.  10). 


DISCUSSION 

In  this  study  we  report  the  identification  of  several  transcriptionally  regulated  genes 
that  were  differentially  expressed  throughout  the  cell  cycle  of  Alexandrium  fundyense  and 
either  up-  or  downregulated  during  toxigenesis.  The  ability  to  obtain  cells  at  different 
times  from  a  single  culture  of  A.  tamarense,  while  toxin  production  was  ON  or  OFF,  was 
essential  to  this  study  and  provided  an  ideal  system  for  identifying  molecular  factors 
related  to  toxigenesis.  The  advantages  of  using  Differential  Display  over  other 
techniques,  such  as  subtractive  hybridization,  to  identify  differentially  expressed  genes 
was  that  it  allowed  the  analysis  of  mRNAs  independently  of  their  prevalence,  only 
required  small  amounts  of  RNA  and  resolved  for  up  and  downregulation  simultaneously. 
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Figure  6.  Nucleotide  and  predicted  amino  acid  sequence  of  the  partial  cDNA  for  a 
putative  i'-adenosyl-homocysteine  hydrolase  in  A.  fundyense.  The  underlined 
sequence  represents  the  DD  band  that  was  originally  isolated.  Stars  mark  the 
two  DD  primer  sites  (A,T  and  A  denote  non-matching  nucleotides  on  the 
'random'  decamer;  T  and  C  mark  non-matching  nucleotides  on  the  reverse 
complementary  sequence  of  the  same  'random'  decamer  (see  discussion)). 
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GTG-AAC~GAC-TGC-GTG-ACA-AAG-TCC“AAG-TTC-GAT-AAT~GTG-TAC-GGG-TGC-CGT~ 

VNDCVTKSKFDNVYGCR 


CAC-TCG-TTG-CCG-GAC-ATC-ATG-CGT-GCC-ACA-GAT-GTC-ATG-ATT-GGG-GGC-AAG- 

HSLPDIMRATDVMIGGK 


CGC-GCC-CTG-ATC-TGC-GGC-TAT-GGC-GAT-GTC-GGC-AAG-GGC-TGT-GCT-TTG-GCC- 

RALICGYGDVGKGCALA 

A  **T  *A*  *** 

ATG-CGC-GGA-GCC-GGA-GCC-CGT-GTC-TTG-ATT-ACC-GAG-ATT-GAC-CCA-ATC-TGC- 
MRGAGARVLITEIDP  IC 


GCG-CTT-CAG-GCG-TGC-ATG-GAG-GGC-TTC-CAG-GTT-GTC-ACT-CTG-GAG-TCC-GTA- 

ALQACMEGFQVVTLESV 

*★*  ★ip*  *Q*  * 

GTT-GGG-GAG-GTG-GAT-ATT-TTT-ACG-ACG-ACC-ACG-GGC-AAC-TTC-AAG-ATC-ATC- 
VGEVDIFTTTTGNFKI  I 


ACT-CTG-GAG”CAT-"ATG-AAA-AAG-ATG“AAG“AAC-AAT~GCC”ATC-"GTC“GGC-AAC-ATT- 

TLEHMKKMKNNAIVGNI 


GGC”CAC-TTC-GAC”AAC-GAG-ATT-GAG-ATG"GCT-GAG-CTG-GAA”GGA-ATG-CCG-GGC- 

GHFDNEIEMAELEGMPG 


ATC-AAG-GTG-GAG-AAC-ATC-AAG-GCA-CAG-GTT-GAC-CGC-TTC-ATA-TTC-CCA-GAT- 

IKVENIKAQVDRFIFPD 


GGC-CAC-GGT-ATC-ATC-GTT-CTC-GCC-TCC-GGC-CGG-CTC-TTG-AAC-TTG-GGC-TGT- 

GHGIIVLASGRLLNLGC 


GCC-ACA-GGC-CAT-CCC-TCC-TTC-GTC-ATG-TCT-TGC-TCC-TTC-ACC-AAT-CAG-GTG- 
ATGHPSFVMSCSFT  NQV 


TTG-GAC-AAC-TTG-ACT-TTG-CTT-CAA 

LDNLTLLQ 
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Figure  7.  Partial  alignment  of  the  Alexandrium  fundyense  deduced  5'-adenosyl- 
homocysteine-hydrolase  amino  acid  sequence  with  those  of  tobacco 
(Nicotiana  sylvestris),  yeast  {Schizosaccharomyces  pombe)  and  a 
cyanobacterium  {Synechocystis  sp.)-  Boxes  denote  higher  than  50%  identity 
(Bold  letters)  or  similarity  (Plain  letters)  among  the  four  oganisms. 
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Figure  8.  Nucleotide  and  predicted  amino  acid  sequence  of  the  3'  end  of  a  cDNA 
encoding  for  a  putative  methionine  aminopeptidase  in  A.  fundyense.  The 
nucleotides  in  lower-case  correspond  to  non-encoding  regions  on  the  3'  end 
of  the  cDNA  strand.  The  underlined  sequence  represents  the  DD  band  that 
was  originally  isolated.  Stars  mark  the  two  DD  primer  sites  (C  and  C  denote 
non-matching  nucleotides  on  the  'random'  decamer;  V  stands  for  A,C  or  G  in 
the  anchored  primer).  The  italicized  triplet  (TGA)  represents  the  first  STOP 
codon  encountered. 
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GTC-AGC-GAT-GTC-GCT-GAG-CGC-GAA-GGG-TTT-GGT-GTT-GTC-AAG-ACC-TTG-GTT- 

VSDVAEREGFGVVKTLV 


GGT-CAT-GGC-ATT-GGC-GAG-TTT-TTC-CAT-GGA-GTG-CCT-CAG-GTC-TTC-CAC-TGT- 

GHGIGEFFHGVPQVFHC 


CGA-AAC-AGC-GAC-AAT-CGA-AAG-ATG-CAG-GAA-GGC-ACC-ACC-TTC-ACC-CTG-GAG- 

RNSDNRKMQEGTTFTLE 


CCA-GTA-CTT-ACA-GAG-GGG-TCA-CGG-GAT-TGG-ATA-ACG-TGG-GAT-GAC-GGG-TGG- 

PVLTEGSRDWITWDDGW 


ACA-GTT-GCT-ACA-AGT-GAC-CAT-GGC-CGC-GCT-GCC-CAA-TTT-GAG-CAC-ACG-CTG- 

TVATSDHGRAAQFEHTL 

*  ***  Q-k*  Qk* 

CTG-ATC-ACT-GCC-AAT-GGC-TGT-GAG-GTC-ATG-ACT-rGAtacatctacaacaacatctaa 

LITANGCEVMT 


cccacaatacaatatacaactcctaccataccaaccactataacctactcccccattcaaaaaaaaaa 

★ Y*  *****★★★* 
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Figure  9.  Partial  alignment  of  the  Alexandrium  fundyense  deduced  methionine 
aminopeptidase  amino  acid  sequence  with  those  of  yeast  (Saccharomyces 
pombe),  a  cyanobacteria  (Synechocystis  sp.)  and  Escherichia  coli.  The 
shaded  areas  correspond  to  three  of  the  five  highly  conserved  copper  binding 
regions  that  have  been  described  for  this  enzyme  (see  text).  Boxes  denote 
higher  than  50%  identity  (Bold  letters)  or  similarity  (Plain  letters)  among  the 
four  organisms. 
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Figure  10.  Nucleotide  sequence  of  the  ORF  encoding  for  a  putative  histone-like  protein 
in  A.  fundyense.  The  nucleotides  in  lower-case  correspond  to  non-encoding 
regions  on  the  5'  and  3'  ends  of  the  cDNA  strand.  The  underlined  sequence 
represents  the  DD  band  that  was  originally  isolated.  Stars  mark  the  two  DD 
primer  sites  (T,G  and  C  denote  non-matching  nucleotides  on  the  'random' 
decamer;  V  stands  for  A,C  or  G  in  the  anchored  primer).  The  italicized 
triplet  {TAG)  represents  the  first  STOP  codon  encountered. 


cacgagcccgttcg-ATG-AGC-CCC-TCG-GCG-ACC-AGA-TCC-CGG-ACG-CGA-ACA-GCG- 

MSPSATRSRTRTA 


GGG-GCA-TCT-CGC-AGG-ACA-AGG-CTG-CGG-ACA-AGT-TCT-GCC-AGA-GGT-TCT-GCC- 

GASRRTRLRTSSARGSA 


CGA-GCT-CCG-GTT-GCC-ACT-GCG-GAG-GTG-AAG-AAG-TTG-GGC-GTC-TTC-ACT-GTC- 

RAPVATAEVKKLGVFTV 


CCA-GGC-CCC-TAC-AAT-ATC-AAG-TTG-AGG-GAG-AAG-CCC-GCC-ACC-CAG-GCC-AGC- 

PGPYNIKLREKPATQAS 


AAG-TGG-AAG-GTG-TGC-CTG-AGA-GTG-GCC-TGG-CAT-GTA-TTC-CCT-TTA-AAG-GCT- 

KWKVCLRVAWHVFPLKA 

*  c** 

TGC-AAG-CTC-TGT-y.AGaacatcaaactqaaaaactttacttaaaattacaactaaaacaacaaaa 
C  K  L  C 


qqctcatctacaqqctcttcttcqcqcttccqatcctqqqccctctaaqqaqaaqaaqqqcaaaatqaaaaq 

taqqaaqtqatqccqqtqatattactcaqctccaqtqqcctqctctcaqaqttqatttqctaqqtcaqtcat 

tttaaactctqtaatqtqcattctacacqcaqtatttactttqtaaatctctcaqqcqatcctcttcattac 

attqtqtqqaqqcqqattaaqacqttqtcaaaaaaaaaaaaaaaaaaaaaaaaa 
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Figure  1 1 .  Alignment  of  the  deduced  amino  acid  sequence  of  the  Alexandrium 
fundyense  histone-like  protein  with  two  different  histone-like  protein 
sequences  from  Crypthecodinium  cohnii.  Boxes  denote  higher  than  50% 
identity  (bold  letters)  or  similarity  (plain  letters)  among  the  three  sequences. 
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Synchronized  Growth  and  Differential  Display.  The  application  of  Differential 
Display  (DD)  to  the  identification  of  genes  that  are  expressed  during  toxin  production  in 
Alexandrium  sp.,  and  hence  potentially  related  to  toxin  biosynthesis,  was  contingent  on 
the  design  of  a  system  that  would  allow  the  comparison  of  toxin-producing  and  non¬ 
producing  cells.  Two  experimental  setups  were  possible:  First,  phylogenetically  close 
toxic  and  non-toxic  strains  of  Alexandrium  sp.  could  be  compared;  alternatively,  cells  of  a 
single  isolate  collected  in  different  physiological  stages  relative  to  toxin  production  could 
be  analyzed.  As  discussed  previously,  mutant  formation  was  not  a  realistic  option. 

Extensive  work  has  been  done  on  the  phylogeny  of  the  Alexandrium  group  using 
ribosomal  RNA  gene  sequences  (Scholin  et  al,  1995).  Based  on  the  phylogenetic 
distance  information  we  ran  some  DD  gels  comparing  toxic  and  non-toxic  strains  of 
Alexandrium  tamarense  that  were  closely  related,  but  their  genetic  variation  was  still 
large  enough  as  to  make  it  impossible  to  identify  differentially  expressed  genes  over  the 
general  variation  between  DD  patterns  due  to  inter-strain  dissimilarities  (Fig.  3). 
Variations  in  DD  patterns  arise  from  differences  in  the  nucleotide  sequences  of  the 
mRNAs.  Changes  in  just  one  nucleotide  in  the  region  binding  to  the  DD  primers  can 
result  in  no  annealing  of  a  specific  primer,  and  hence  no  amplification  of  the  mRNA. 
This  same  mRNA  might  be  amplified  by  a  different  set  of  primers.  Therefore,  the 
absence  of  a  specific  band  does  not  mean  that  the  corresponding  gene  is  not  expressed  but 
that  it  is  not  amplified  by  this  particular  set  of  primers.  The  similarity  among  the  DD 
patterns  has  to  be  of  at  least  90-95  per  cent  to  make  the  comparison  meaningful. 
Variabilities  of  the  magnitude  observed  in  Figure  3  result  in  an  impossibility  to  make  any 
prediction  about  the  actual  regulation  of  the  genes  identified  as  differentially  expressed 
because  of  the  uncertainty  about  the  identity  between  similar  bands.  Interspecific 
variations  between  A.  tamarense  and  A.  fundyense  were  equally  uninformative  (Fig.  3). 
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The  sensitivity  of  DD  to  changes  in  sequences  was  such  that  the  cells  to  be  compared  had 
to  be  truly  identical  in  their  genetic  makeup  to  avoid  a  noisy  background. 

The  alternative  was  to  take  advantage  of  the  discontinuity  in  toxin  biosynthesis  over 
the  cell  cycle  in  synchronously  growing  cells  of  A.  fundyense  that  we  had  previously 
reported  (Taroncher-Oldenburg  et  al,  1997).  When  we  performed  DD  on  samples  drawn 
at  different  times  throughout  the  cell  cycle,  the  patterns  of  the  reverse  transcribed  and 
amplified  mRNA  components  analyzed  were  largely  identical,  resulting  in  a  uniform 
background  over  which  specific,  differentially  expressed  genes  could  be  identified  (Figs. 
3«&4). 

S-Adenosylhomocysteine  Hydrolase.  The  enzyme  5-adenosylhomocysteine  hydrolase 
(SAHH)  has  a  key  role  in  the  regulation  of  the  S-adenosylmethionine  (SAM)  pathway, 
one  of  the  main  elements  of  one-carbon  metabolism  in  the  cell.  It  catalyzes  the  reversible 
cleavage  of  ^-adenosylhomocysteine,  the  side  product  of  SAM-mediated  methylation 
reactions,  into  adenosine  and  homocysteine.  The  equilibrium  of  this  reaction  lies 
distinctly  in  favor  of  the  synthesis  of  5-adenosylhomocysteine.  The  homocysteine 
resulting  from  the  cleavage  reaction  is  subsequently  remethylated  to  yield  methionine. 
This  last  reaction  is  catalyzed  by  homocysteine  methyltransferase,  a  methyl- 
tetrahydrofolate  dependent  enzyme.  In  the  closing  step  of  the  SAM  cycle,  methionine  is 
reactivated  by  binding  to  ATP  to  form  SAM.  By  virtue  of  its  reversibility,  the  function  of 
the  SAHH  is  to  control  the  recycling  of  S-adenosylhomocysteine  into  the  SAM  cycle  to 
ensure  a  steady  supply  of  homocysteine  for  its  transmethylation  into  methionine. 
Similarily  to  our  observation  in  Alexandrium  fundyense,  this  enzyme  has  also  been  found 
to  be  transcriptionally  regulated  in  other  organisms.  SAHH  expression  is  induced  by 
pathogenic  attack  in  parsely  (Kawalleck  et  al,  1992)  and  differentially  expressed  over  the 
progression  of  developmental  stages  in  lucerne  (Abrahams  et  al.,  1995). 
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The  connection  SAHH  downregulation  in  A.  fundyense  might  bear  to  STX 
biosynthesis  lies  in  the  postulated  presence  of  one  SAM  dependent  methylation  step  in 
the  early  phases  of  the  synthesis  of  the  STX  backbone  (Shimizu,  1993).  By  regulating 
the  supply  of  SAM,  this  pathway  could  limit  or  enhance  the  activity  of  the  toxigenic 
enzymes.  Such  a  correlation,  though,  would  most  probably  be  non-specific  for  toxin 
biosynthesis,  since  SAM  is  involved  in  a  wide  spectrum  of  anabolic  processes  in  the  cell. 

Methionine  Amino  Peptidase.  Methionine  aminopeptidase  (MAP)  is  responsible  for 
cleaving  the  amino-terminal  methionine  during  processing  of  most  nascent  proteins.  It  is 
cobalt  dependent  and  its  active  site  contains  two  zinc  ion  binding  domains  (Arfm  et  al, 
1995).  The  partial  Alexandrium  sequence  reported  here  contains  three  of  the  five  highly 
conserved  domains  responsible  for  its  metal-binding  character.  The  function  of  this 
enzyme  in  other  organisms  clearly  relates  to  regulation  of  protein  synthesis  and  enzyme 
activities  in  general  and  is  essential  for  the  viability  of  the  organism.  MAPs  are  present 
and  active  throughout  the  entire  life  cycle  of  a  cell.  Depending  on  the  organism,  the 
number  of  different  MAPs  in  a  cell  can  vary  depending  on  the  developmetal  status  or  cell 
cycle  stage  of  the  cell  (Arfin  et  al,  1995).  Regulation  further  happens  at  both  the 
transcriptional  and  post-translational  levels  in  different  organisms  (Arfm  et  al.,  1995). 

Protein  biosynthesis  regulation  at  this  primary  level  might  have  a  direct  influence  on 
toxin  biosynthesis  via  complex  interactions  between  different  components  of  the  general 
cellular  metabolism,  but  no  straight  correlation  to  STX  biosynthesis  can  be  inferred. 
Other  scenarios  cannot  be  completely  excluded  as  recent  reports  have  shown  that  in 
addition  to  its  protein  editing  function,  MAP  also  acts  as  an  inhibitor  of  eukaryotic 
initiation  factor  2a  (elF-2a)  phosphorylation  and  is  involved  in  regulation  of 
angiogenesis  in  higher  eukaryotes  (Griffith  et  al.,  1997),  implying  that  not  all  the 
functions  of  this  gene  are  known  yet. 
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Histone-like  Protein.  Dinoflagellate  nuclei  are  unique  among  eukaryotes  in  that  their 
chromosomes  remain  condensed  throughout  the  entire  cell  cycle,  in  contrast  to  the 
unwinding  and  relaxation  of  the  chromatin  observed  in  all  other  eukaryotic  organisms 
(reviewed  (Rizzo,  1991)).  In  addition  to  this,  dinoflagellates  also  lack  the  typical  histone 
proteins  and  nucleosome  structures  that  define  typical  nuclei  of  other  eukaryotic 
organisms  (Herzog  and  Soyer,  1981;  Rizzo,  1981).  Nonetheless,  dinoflagellate 
chromosomes  have  been  shown  to  contain  an  alternative  suite  of  basic  proteins,  the  so 
called  histone-like  proteins  (Rizzo  and  Morris,  1984;  Rizzo,  1987).  No  clear  function  has 
been  adscribed  to  these  proteins  yet,  but  the  evidence  gathered  so  far  suggests  that  rather 
than  having  an  intrinsic  structural  or  stabilizing  function  they  might  be  involved  at  some 
level  in  transcriptional  regulation  or,  perhaps,  initiation  (Sala-Rovira  et  al,  1991). 

Band  C  of  the  Differential  Display  analysis  rescued  a  full-length  cDNA  with  a  high 
similarity  at  the  amino  acid  level  with  histone-like  proteins  from  the  dinoflagellate 
Crypthecodinium  cohnii.  Differential  or  cyclic  'de  novo'  synthesis  of  histone-like  proteins 
in  C.  cohnii  has  been  suggested  based  on  immunocytochemical  studies  (Geraud  et  al, 
1991).  No  information  is  available  on  histone-like  proteins  in  A.  fundyense. 

No  direct  correlation  can  be  drawn  in  this  case  between  the  induction  of  the  histone¬ 
like  protein  and  STX  biosynthesis.  On  the  other  hand,  the  isolation  of  this  protein 
illustrates  the  fact  that  DD  analysis  of  any  population  of  synchronized  cells  is,  by 
definition,  an  ideal  system  to  identify  genes  related  to  cell  cycle  progression.  As 
mentioned  above,  histone-like  proteins  are  closely  associated  to  DNA  and  its  structure  or 
some  of  its  functions,  which  means  that  its  biosynthesis  is  probably  regulated  in  a  DNA 
or  cell  cycle  dependent  manner.  Cell  cycle-dependent  expression  has,  for  example,  been 
shown  for  histones  in  plants  (Callard  and  Mazzolini,  1997). 
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Regulation  of  STX  biosynthesis.  Biosynthetic  pathways  can  be  regulated  at  many 
different  levels.  Expression  of  a  specific  enzyme  can  be  regulated  at  the  transcriptional 
level,  through  the  induction  or  suppresion  of  gene  transcription  itself,  or  at  the  pre- 
translational  level,  which  includes  mRNA  editing,  regulation  of  translation  and  post- 
translational  mechanisms  such  as  protein  maturation  and  enzyme  inhibition.  Not  much  is 
known  about  regulation  of  STX  biosynthesis  in  dinoflagellates.  Enhancement  and 
inhibition  of  toxigenesis  have  been  shown  under  varying  environmental  conditions  (Hall, 
1982;  Boyer  et  al,  1987;  Anderson  et  al,  1990).  Toxin  content  and  variations  in 
generation  time  or  length  of  the  cell  cycle  have  been  suggested  to  be  correlated 
(Anderson  et  al.,  1990).  More  recently,  STX  biosynthesis  has  also  been  shown  to  be 
coupled  to  the  Gi  phase  of  the  cell  cycle  (Taroncher-Oldenburg  et  al,  1997).  All  this 
evidence  led  us  to  the  hypothesis  that  STX  biosynthesis  might  be  regulated  in  a  cell 
cycle-dependent  manner. 

Many  cell  cycle-dependent  genes  are  known  to  be  transcriptionally  regulated  in 
eukaryotes  (Price  et  al,  1991;  Callard  and  Mazzolini,  1997).  Based  on  this  and  our 
ability  to  synchronize  cultures  of  Alexandrium  fundyense,  we  set  out  to  identify  genes 
regulated  in  patterns  that  correlated  with  episodes  of  toxigenesis  during  the  early  part  of 
the  Gi  phase  of  the  cell  cycle.  The  results  of  this  search  were  expected  to  fall  within  one 
of  four  possible  categories:  (1)  Gene  products  which  could  be  involved  in  one  of  the 
proposed  steps  of  STX  biosynthesis  (Shimizu,  1993);  (2)  genes  with  known,  cell  cycle 
related  functions  (Callard  and  Mazzolini,  1997);  (3)  genes  encoding  regulatory  factors, 
cell  cycle  related  or  not,  that  could  potentially  activate  toxin  production  (Tandeau  de 
Marsac  and  Houmard,  1993;  Van  Dolah  et  al,  1995;  Custer  et  al,  1997);  and  (4)  novel 
genes  with  unknown  functions  that  could  be  related  to  toxigenesis  and/or  cell  cycle 
progression  (Joseph  et  al,  1994;  Utans  et  al,  1994).  None  of  the  three  genes  we  isolated 
seemed  to  specifically  relate  to  STX  biosynthesis.  One  gene,  encoding  the  putative 


160 


histone-like  protein  HAf,  could  be  a  cell  cycle-related  gene  since  it  encodes  for  a  protein 
that  is  surmised  to  have  DNA  structural  and  regulatory  functions.  The  two  other  genes, 
SAHH  and  MAP,  are  both  involved  in  regulatory  pathways.  One,  MAP,  acts  at  the  more 
general  level  of  post-translational  regulation.  The  other,  SAHH,  regulates  the  supply  of 
methylating  units,  SAM,  in  the  cell  and  hence  affects  not  just  the  methylation  step  in  STX 
biosynthesis  but  also  other  methyl-dependent  pathways  in  the  cells.  No  genes  with 
unknown  functions  were  isolated. 

From  the  results  presented  here  we  can  conclude  several  things  related  to  STX 
biosynthesis  and  its  regulatory  mechanism  in  Alexandrium  fundyense.  First,  and  despite 
the  fact  that  we  did  not  identify  specific  enzymes  for  STX  biosynthesis,  we  cannot 
exclude  the  possibility  that  this  process  is  still  transcriptionally  regulated.  DD  is  the  most 
reliable  and  representative  method  available  to  identify  transcriptionally  regulated  genes 
in  eukaryotes,  but  it  is  based  on  PCR  and  thereby  introduces  a  statistical  probability  of 
missing  and  not  amplifying  rare  transcripts  (Liang  et  al,  1993;  Wan  et  al,  1996). 
Secondly,  the  genes  involved  in  STX  biosynthesis  could  be  post-transcriptionally 
regulated.  This  could  happen  at  the  translational  level,  in  which  case  the  stable 
transcripts  could  be  present  throughout  the  cell  cycle  and  only  translated  at  specific  times, 
or  it  could  be  regulated  at  the  post-translational  level. 

Few  data  exist  on  genetic  regulation  in  dinoflagellates.  The  best  studied  system  is  the 
regulation  of  luciferin  binding  protein  (LBP)  in  Gonyaulax  polyedra  (Machabee  et  al, 
1994;  Mittag  et  al,  1994).  The  synthesis  of  LBP  is  under  circadian  control  and  its 
corresponding  mRNA  is  present  throughout  the  cell  cycle.  Translation  is  repressed  in  the 
presence  of  a  protein  that  binds  specifically  to  a  3'-untranslated  region  (UTR)  of  the 
mRNA.  The  regulation  of  two  other  families  of  proteins,  the  peridinin  chlorophyll-a 
binding  proteins  (PCPs)  and  the  chlorophyll-a-c  binding  proteins  (ACPs),  has  been 
studied  at  a  physiological  level  similar  to  STX  production  (Jovine  et  al,  1992;  Triplett  et 
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al,  1993;  Johnsen  et  al,  1994).  The  activation  pattern  of  the  cell  cycle  regulatory  protein 
CDC2  kinase  has  been  characterized  in  Gambierdicus  toxicus  (Van  Dolah  et  al.,  1995). 
CDC2  kinase  is  constitutively  expressed  throughout  the  cell  cycle  and  activated  only  in 
mitotic  cells.  Evidence  for  the  presence  of  this  ubiquitous  regulator  of  cell  cycle 
progression  in  eukaryotes  has  also  been  gathered  in  Crypthecodinium  cohnii  (Rodriguez 
etai,  1994). 

The  common  denominator  to  all  the  previous  systems  is  translational  regulation.  In 
this  study  we  show  for  the  first  time  that  transcriptional  or  pre-translational  regulation 
does  indeed  occur  in  dinoflagellates.  Previous  reports  had  suggested  the  possibility  that 
translational  mechanisms  could  be  the  primary  means  of  circadian  regulation  in  these 
organisms  (Milos  et  al,  1990).  We  have  shown  not  only  the  induction  or  suppression  of 
the  three  genes  described  here,  but  also  observed  transcriptional  or  pre-translational 
regulation  of  about  twenty  other  genes  with  varying  DD  patterns  (e.g.,  band  E,  Fig.  4). 

Differential  Display  and  Dinoflagellates.  This  study  illustrates  the  potential  of 
applying  Differential  Display  to  the  identification  of  transcriptionally  regulated  genes  in 
dinoflagellates  and  other  phytoplankton.  Many  of  these  organisms  cannot  be  cultivated 
or,  when  finally  cultivated,  do  not  produce  enough  material  for  traditional  molecular 
studies.  This  is  not  a  problem  with  DD  due  to  the  amplification  step  inherent  to  the 
method.  DD  also  circumvents  the  need  to  create  mutants  for  the  analysis  of  specific 
genetic  functions.  In  this  study,  we  adapted  DD  to  the  analysis  of  different  cell  cycle 
stages  in  Alexandrium.  Other  applications  include  the  study  of  different  life  cycle  stages 
of  phytoplankton  or  the  isolation  of  specific  genes  induced  by  environmental  stresses. 
Examples  for  these  applications  on  higher  eukaryotes  already  exist  in  the  literature. 
Sexual-stage  specific  mRNAs  have  been  isolated  in  Eimeria  bovis,  a  parasitic  protozoan, 
and  genes  involved  in  the  different  steps  of  hematopoiesis  have  been  identified  in  mouse 
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(Abrahamsen  et  al,  1995;  Guimaraes  et  al,  1997).  Tissue  specific  markers  have  been 
observed  in  the  rhodophyte  Porphyra  perforata  (Hong  et  al,  1995).  Genes  induced  by 
ozone  stress  in  Arabidopsis  thaliana  and  genes  induced  by  phosphate  depletion  in  rat 
kidney  have  been  isolated  (Sharma  and  Davis,  1995;  Custer  et  al,  1997). 

In  conclusion,  transcriptional  or  translational  mechanisms  could  underly  the 
regulation  of  STX  biosynthesis,  but  it  will  be  necessary  to  gather  more  information  on  the 
biosynthetic  pathway  in  order  to  identify  the  enzymes  involved  and  study  their  expression 
patterns  as  well  as  those  of  the  inferred  genes.  This  study  has  also  documented  the 
potential  of  applying  DD  to  the  identification  of  genes  that  are  related  to  certain 
physiological  or  developmental  processes  in  phytoplankton. 
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The  work  described  in  this  thesis  has  provided  new  insights  into  the  cell  cycle 
dynamics  and  the  physiology  of  STX  biosynthesis  in  the  dinoflagellate  Alexandrium 
fundyense.  This  was  accomplished  by  combining  a  new  approach,  namely  the  study  of 
toxigenesis  in  correlation  to  the  cell  cycle  in  synchronized  cells,  and  the  application  of 
Differential  Display  to  analyze  differences  in  expressed  mRNAs.  In  contrast  to  prior 
studies,  this  approach  allowed  the  analysis  of  toxin  production  over  short  time  intervals  and 
relative  to  a  well-defined  cellular  process,  the  cell  cycle.  It  further  allowed  the  identification 
of  genes  which  had  expression  patterns  related  to  toxigenic  periods.  Overall,  this  work  has 
shown  the  potential  of  using  traditional  physiological  techniques  and  molecular  approaches 
to  study  the  dynamics  and  the  regulation  of  biochemical  pathways  in  phytoplankton. 

The  first  chapter  presents  the  results  from  a  study  of  toxigenesis  in  synchronized 
cultures  of  Alexandrium  fundyense.  For  the  first  time  it  was  possible  to  determine  that 
toxin  production  is  discontinuous  and  coupled  to  the  Gi  phase  of  the  cell  cycle.  The 
differential  accumulation  of  STX  derivatives  over  time  was  also  interpreted  in  terms  of  a 
possible  sequence  of  interconversions  among  these  compounds.  I  further  demonstrated  the 
existence  of  a  discrete  S  phase  and  the  presence  of  two  transition  points  in  the  cell  cycle  of 
A.  fundyense.  These  results  implied  cell  cycle  regulation  of  STX  biosynthesis  and  a 
simple  scheme  of  interconversions  among  STX  derivatives. 

Chapter  11  examines  the  effects  of  temperature  and  phosphate  limitation  on  the  coupling 
of  STX  biosynthesis  to  the  Gi  phase  of  the  cell  cycle.  Over  a  range  of  growth  rates  and 
experimental  treatments,  toxin  accumulation  always  correlated  to  increases  in  duration  of 
Gi.  Therefore,  and  contrary  to  what  had  been  previously  suggested,  toxin  production  rates 
were  independent  of  the  growth  rate  in  A.  fundyense.  Differences  in  toxin  composition 
between  phophate-limited  and  nutrient-replete  cultures  were  further  interpreted  to  show  that 
interconversions  among  STX  derivatives  by  sulfatation  reactions  are  not  phosphate- 
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dependent,  while  those  mediated  by  oxidation  reactions  are  inhibited  by  phosphate 
limitation. 

In  chapter  III,  Differential  Display  (DD)  analysis  of  mRNA  was  used  to  identify  genes 
that  were  up-  or  downregulated  during  the  now  well  characterized  periods  of  STX 
biosynthesis  in  synchronized  cultures  of  A.  fundyense.  Three  mRNAs  were  characterized 
and  their  expression  patterns  confirmed.  This  study  further  showed  unequivocally  that 
transcriptional  regulation  occurs  in  dinoflagellates,  contrasting  with  previous  suggestions 
that  translational  mechanisms  might  be  the  primary  means  of  circadian  regulation  in  these 
organisms.  Despite  the  fact  that  none  of  the  genes  that  were  isolated  was  directly  and 
exclusively  related  to  toxigenesis,  this  work  revealed  the  potential  of  using  DD  and  other 
similar  techniques  for  the  identification  of  differentially  expressed  genes  in  dinoflagellates 
and  phytoplankton  in  general. 

Overall,  this  study  clarifies  several  aspects  of  the  biochemistry,  physiology  and 
regulation  of  STX  biosynthesis  and  opens  new  paths  for  future  research  that  I  will  present 
in  detail  now. 

Biochemistry  of  STX  biosynthesis.  The  main  problem  remaining  to  be  solved  in  this 
field  is  the  elucidation  of  the  biosynthetic  pathway  of  STX.  The  primary  precursors, 
arginine,  acetate  and  S-adenosylmethionine  have  been  determined,  but  these,  especially 
arginine,  could  be  incorporated  in  a  variety  of  ways  which  have  yet  to  be  determined.  This 
has  to  happen  at  two  levels:  (1)  determine  the  direct  precursors  or  first  intermediates  leading 
to  the  synthesis  of  the  STX  backbone,  and  (2),  establish  the  sequence  of  interconversions 
among  STX  derivatives.  The  first  is  the  more  difficult  to  address  as  it  requires  possible 
rethinking  of  the  biogenetic  pathway  proposed  by  Shimizu  and  the  development  of 
alternative  pathways  that  include  all  of  the  available  information  on  precursor  incorporation. 
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Fig.  1  Proposed  scheme  of  interconversions  among  SIX  derivatives.  A  putative  origin  is 
located  around  do-dc-STX,  the  simplest  of  all  known  derivatives  (Broad  arrows 
indicate  proven  conversions  (refer  to  the  INTRODUCTION);  thin  arrows  denote 
putative  reactions). 
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The  recent  isolation  of  simple  derivatives  of  SIX  such  as  deoxy-decarbamoyl  STX 
indicates  that  perhaps  the  first  compounds  produced  by  the  biosynthetic  pathway  are 
centered  around  this  kind  of  simple  molecule.  The  breadth  of  information  on  toxin 
compositions  and  the  few  data  available  on  enzymes  that  specifically  catalyze  some  of  the 
conversions  also  indicate  that  the  first  toxins  to  appear  might  be  those  of  the  STX  series  (no 
OH  on  N-1).  This  hypothesis  is  supported  by  several  facts:  (1)  No  dc-NEO  or  dc-GTX 
I, IV  has  been  isolated  from  any  toxin  synthesizing  organism,  (2)  the  two  sulfatases  so  far 
described  in  the  literature  are  specific  for  STX  and  GTX II, III,  (3)  the  existence  of  a  STX- 
specific  oxidase/hydroxylase  activity,  (4)  the  selective  increase  of  C  I,II  and  GTX  II,III 
under  phosphate-limited  conditions  as  described  in  chapter  II  and  (5)  the  isolation  of  do-dc- 
STX  and  do-dc-GTX  II,  III  but  no  do-dc-NEO  or  do-dc-GTX  I,IV.  The  scheme  of 
interconversions  might  be  further  simplified  by  the  fact  that  11-P-epimers  are  possibly  the 
only  derivatives  synthesized.  All  of  the  preceding  information  can  be  summarized  in  a 
scheme  of  interconversions  in  which  the  first  compound  is  located  somewhere  in  the 
'structural'  vicinity  of  do-dc-STX  and  the  interconversions  among  the  derivatives  follow  a 
few  possible  routes  (Fig.  1).  Only  three  classes  of  enzymes  would  be  involved  in  this 
proposed  pathway;  (1)  oxidases/hydroxylases  for  the  oxydation  of  the  do-dc-compounds 
and  to  generate  the  hydroxy-  or  NEO-series  of  toxins;  (2)  a  carbamoyltransferase  for  the 
carbamoylation  of  the  C-13  residues;  (3)  sulfatases  responsible  for  sulfatetransfers  on  two 
different  positions  of  the  derivatives.  The  question  still  remains  open  as  to  how  exactly  that 
first  parent  compound  is  synthesized  and  what  its  nature  is. 

Physiology  of  STX.  One  of  the  challenges  in  studying  the  physiology  and  dynamics  of 
STX  production  is  the  fact  that  we  still  do  not  know  the  specifics  of  the  biosynthetic 
pathway.  This  makes  it  difficult  to  interpret  some  of  the  information  accumulated  over  the 
last  decades.  Very  careful  experiments  describing  the  changes  in  toxicity  in  response  to 
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variations  in  the  physico-chemical  environment  of  the  cells  provide  an  insight  into  the 
effects  of  these  environmental  variables  on  toxigenesis,  but  the  general  nature  of  these 
factors  and  their  broad  influence  on  the  cell's  physiology  make  it  hard  to  determine  their 
mechanisms  of  action.  This  task  is  further  complicated  by  the  fact  that  regulation  of  toxin 
production  might  happen  at  many  different  levels,  making  it  complicated  to  isolate  the  effect 
of  single  factors  on  its  pathway.  It  will  therefore  be  imperative  in  the  future  to  learn  more 
about  the  biosynthetic  pathway  in  order  to  fully  understand  its  regulation.  One  intriguing 
result  presented  in  this  thesis  is  related  to  the  toxin  production  kinetics  curve  shown  in  the 
first  chapter.  The  shape  of  this  curve  suggests  the  possibility  of  a  biphasic  kinetics  pattern 
in  which  an  initial  burst  in  biosynthetic  activity  is  followed  by  a  steady  state  kinetics  at  a 
lower  rate.  Such  dynamic  patterns  will  be  understood  better  once  the  actual  enzymes 
responsible  for  STX  production  are  isolated  and  purified. 

Genetics  and  regulation  of  STX  biosynthesis.  This  is  still  an  unchartered  territory.  In 
chapter  III  of  this  thesis  I  proposed  a  new  approach  to  the  identification  of  genes  related  to 
STX  biosynthesis  based  on  the  hypothesis  that  toxin  production  is  cell  cycle  regulated  and 
hence  possibly  regulated  at  the  transcriptional  level.  This  study  showed  the  strengths  and 
weaknesses  of  such  an  approach.  As  mentioned  earlier,  we  still  lack  sufficient  knowledge 
about  the  biosynthetic  pathway  to  look  for  specific  enzymes  or  genes.  Blind  approaches  to 
the  isolation  of  putative  'toxin'  genes,  such  as  the  one  described  in  this  thesis,  are  hard 
because  we  still  do  not  have  the  ideal  system  for  singling  out  what  makes  toxic  cells 
different  from  non-toxic  ones.  The  use  of  synchronized  cultures  and  the  fact  that  toxin 
production  is  discontinuous  represents  a  great  leap  forward.  Future  studies  should 
probably  also  look  at  other  organisms,  specifically  cyanobacteria,  which  are  easier  to 
manipulate  genetically.  Mating  experiments  between  toxic  and  non-toxic  strains  of 
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dinoflagellates  also  have  the  potential  to  tell  us  more  about  the  genetics  and  inheritance  of 
toxigenesis. 

Differential  Display  and  phytoplankton.  This  thesis  also  illustrates  the  potential  of  using 
comparative  molecular  techniques  to  study  gene  expression  in  organisms  where  standard 
genetic  approaches  are  not  applicable  due  to  biological  and  physiscal  constraints.  The 
reaction  of  a  microorganism  to  environmental  stress  could  be  monitored  and  specific  genes 
isolated.  I  mentioned  earlier  that  the  study  of  regulation  at  the  molecular  level  requires  a 
good  understanding  of  the  biochemical  mechanisms  underlying  a  metabolic  process.  In 
cases  such  as  nitrogen  and  phosphate  depletion,  quite  a  bit  of  information  is  available  on 
the  underlying  mechanisms.  This  would  make  it  simpler  to  design  and  subsequently 
interpret  experiments  where  variations  in  the  nutritional  status  are  the  only  changing  factor. 
For  many  such  processes,  transcriptional  regulation  is  also  the  main  means  of  regulation. 
Genes  identified  by  DD  or  any  of  the  other  emerging  molecular  techniques  that  allow  the 
identification  of  differentially  regulated  genes,  could  serve  as  markers  for  nutrient  status  of 
the  cells  or  the  environment.  Molecular  meirkers  for  pollution  could  also  be  developed  in  a 
similar  way. 

In  conclusion,  this  thesis  provides  new  insights  on  the  cell  cycle  dynamics  and  the 
physiology  of  saxitoxin  biosynthesis  in  the  dinoflagellate  Alexandrium  fundyense.  The 
correlation  between  these  two  processes  has  been  characterized  and  several  genes  that  are 
regulated  in  a  similar  cell  cycle-dependent  way  were  isolated.  One  of  these  genes,  a 
putative  histone-like  protein,  was  analyzed  in  detail  and  constitutes  an  appendix  to  this 
thesis.  A  possible  scheme  of  interconversion  among  STX  derivatives  is  further  proposed 
and  the  regulatory  function  of  temperature  and  phosphate  on  STX  biosynthesis  were 
analyzed.  Finally,  future  directions  in  the  study  of  STX  production  are  proposed  that 
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converge  on  one  focal  point:  the  primary  importance  to  elucidate  the  sequence  of  steps  of 
the  STX  biosynthetic  pathway  in  order  to  understand  the  significance  of  this  metabolic 
process  for  the  producing  organisms  and  its  regulation,  both,  at  the  biochemical  level  and 
by  the  physico-chemical  environment  of  the  cell. 
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APPENDIX 


Cloning  and  characterization  of  a  putative  histone-like 
protein  in  Alexandrium  fundyense  (Dinophyceae) 
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ABSTRACT 


The  biochemical  characteristics  of  a  putative  basic  histone-like  protein  (HAf)  in  the 
toxic  dinoflagellate  Alexandrium  fundyense  Balech  and  its  expression  pattern  during  the 
cell  cycle  were  determined.  HAf  was  first  identified  by  Differential  Display  (DD)  of 
mRNA  applied  to  the  analysis  of  synchronized  cells  of  A.  fundyense  in  different  cell  cycle 
stages.  Several  mRNAs  expressed  only  at  the  beginning  of  the  Gi  phase  of  the  cell  cycle 
were  cloned  and  used  to  screen  an  A.  fundyense  cDNA  library.  The  deduced  amino  acid 
sequence  of  one  of  these  bands  corresponded  to  HAf.  HAf  is  a  basic  protein  with  a  60% 
similarity  with  two  histone-like  proteins  from  the  dinoflagellate  Crypthecodinium  cohnii 
Biecheler,  was  identified  by  screening  an  A.  fundyense  cDNA  library.  Several  of  its 
properties,  such  as  a  molecular  weight  of  9.3  kDa,  a  basicity  of  24%  and  its  amino  acid 
composition,  were  comparable  to  those  of  two  other  dinoflagellates,C.  cohnii  and  Oxyrrhis 
marina  Dujardin.  The  predicted  secondary  structure,  antigenicity  and  hydropathy  profiles 
were  analogous  to  those  of  two  variants  of  histone-like  proteins  in  C.  cohnii.  Moreover, 
induction  of  HAf  was  observed  only  at  the  beginning  of  the  Gi  phase  of  the  cell  cycle. 
This  pattern  of  transcriptional  regulation  agreed  with  direct  measurements  of  accumulation 
patterns  of  histone-like  proteins  in  C.  cohnii  cultures.  No  defininte  conclusions  can  be 
drawn  about  the  function  of  this  putative  protein,  but  its  similarity  to  other  histone-like 
proteins,  its  basicity  and  its  induction  at  the  beginning  of  the  Gi  phase  of  the  cell  cycle, 
suggest  that  it  could  have  a  role  in  DNA  structure  determination  or  transcription  initiation. 
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INTRODUCTION 


The  distinct  nature  of  the  dinoflagellate  nucleus  is  a  well-documented  phenomenon 
(reviewed  Rizzo,  1987).  Dinoflagellates  represent  the  only  group  of  eukaryotes  lacking 
nucleosome  structures  and  the  associated  histones  (Bodansky  et  al,  1979;  Herzog  and 
Soyer,  1981).  Other  features  include  permanently  condensed  chomosomes,  the  presence 
of  a  nuclear  membrane  throughout  the  cell  cycle  and  an  extranuclear  mitotic  spindle  (Kubai 
and  Ris,  1969).  These  morphological  characteristics  had  long  been  considered  diagnostic 
of  a  primitive  origin  of  the  dinoflagellates  that  placed  them  somewhere  intermediate 
between  prokaryotes  and  'true'  eukaryotes  (Dodge,  1965;  Taylor,  1978).  Only  more 
recently  has  the  phylogenetic  position  of  the  dinoflagellates  been  clarified  through  the 
analysis  of  large  subunit  rRNA  and  their  origin  placed  closer  to  the  emergence  of  other 
eukaryotes  such  as  ciliates  and  yeasts  (Lenaers  et  al,  1990). 

In  addition  to  the  distinct  morphological  features  described  above,  dinoflagellate  nuclei 
also  show  several  biochemical  pecularities.  The  DNA  composition  is  characterized  by  the 
presence  of  a  modified  base,  5-hydroxymethyluracil  (HOMeU),  that  is  not  found  in  other 
eukaryotes  (Rae,  1976).  The  presence  of  unusually  high  levels  of  DNA-bound  metal 
elements,  such  as  copper,  iron,  nickel  and  zinc,  has  been  suggested  to  play  a  structural  role 
in  chromosomes  (Sigee  and  Kearns,  1981).  The  existence  of  a  distinct  group  of  low 
molecular  weight,  basic  proteins  associated  with  nuclear  DNA  has  also  been  reported 
(Rizzo,  1981).  These  proteins  were  termed  'histone-like'  proteins  because,  like  histones, 
they  were  localized  in  the  nucleus  and  were  associated  with  DNA,  but  showed  less  affinity 
for  DNA  and  had  compositional  and  biochemical  differences  relative  to  the  five  typical 
eukaryotic  histones  (Rizzo,  1985). 

No  definite  role  has  been  adscribed  to  histone-like  proteins.  Two  non-exclusive 
hypotheses  have  been  proposed;  a  structural  role  in  DNA  compaction  and  a  regulatory 
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function  in  DNA  transcription  (Sala-Rovira  et  al,  1991).  Immunocytochemical  studies  of 
the  Crypthecodinium  cohnii  histone-like  protein,  HCc,  have  shown  a  homogeneous  and 
close  association  of  this  protein  with  DNA  during  interphase  and  mitosis,  suggesting 
organization  or  stabilization  of  chromosome  structure  (Geraud  et  al,  1991).  HCc  proteins 
also  form  intra-  and/or  intermolecular  dimers  via  disulphide  bridges  between  cysteines  that 
could  have  a  further  DNA  stabilizing  role  (Vemet  et  al,  1990).  In  view  of  the  high  degree 
of  chromosomal  complexity  in  dinoflagellates  and  their  lack  of  regular  eukaryotic 
chromosome  organizers  (nucleosomes)  such  a  structural  function  is  an  attractive 
proposition. 

A  closer  look  at  distribution  patterns  of  HCc  in  chromosomes  showed  that  although 
detectable  over  the  entire  structure,  HCc  was  predominantly  bound  to  extrachromosomal 
DNA  loops  (Geraud  et  al,  1991).  Extrachromosomal  DNA  loops  in  dinoflagellates  have 
been  surmised  to  be  the  transcriptionally  active  regions  of  the  chromosomes  based  on 
radioactive  labeling  experiments,  the  presence  of  Z-DNA  and  the  detection  of  specific 
coding  sequences  in  these  regions  (Sigee,  1984;  Soyer-Gobillard  et  al,  1990;  Anderson  et 
al,  1992).  HCc  has  also  been  observed  in  nucleolar  organizing  regions,  the  sites  where 
ribosomal  transcription  is  initiated  (Geraud  et  al,  1991).  All  this  evidence  suggests  a 
possible  role  of  these  small  basic  peptides  in  transcription  initiation  or  regulation. 

Histone-like  proteins  from  several  dinoflagellates  have  been  characterized  with  respect 
to  amino  acid  composition  and  peptide  mapping  (Rizzo  and  Nooden,  1974;  Jing-Yan, 
1984;  Rizzo  et  al,  1988;  Sala-Rovira  et  al,  1991).  Few  data  exist  at  the  molecular  level, 
as  only  two  sequences  of  histone-like  protein  variants  in  Crypthecodinium  cohnii  have  been 
determined  so  far  (Sala-Rovira  et  al,  1991).  Some  information  is  also  available  on  the 
temporal  fluctuations  of  HCc  over  the  cell  cycle  in  C.  cohnii  (Geraud  et  al,  1991).  In  this 
report  we  describe  a  new  putative  histone-like  protein  from  the  thecate  dinoflagellate 
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Alexandrium  fundyense.  Expression  patterns,  sequence  information  and  a  comparative 
compositional  and  structural  analysis  with  other  histone-like  proteins  are  presented. 


MATERIALS  AND  METHODS 

The  organism.  The  toxic  dinoflagellate  Alexandrium  fundyense  (strain  GtCA29)  was 
used  in  this  study.  Cultures  were  maintained  at  20"C  in  f/2  medium  (Guillard  and  Ryther, 
1962),  modified  by  addition  of  H2Se03  (10‘^M)  and  by  reducing  the  concentration  of 
CUSO4  to  10‘^M.  Irradiance  of  ca.  250  pEm'^s"^  was  provided  by  cool  white  fluorescent 
bulbs  on  a  14:10  h  lightrdark  cycle.  Exponentially  growing  cultures  (approx.  400  cells  ml" 
were  synchronized  by  a  dark  induced  block/release  method,  as  previously  described 
(Taroncher-Oldenburg  etal,  1997).  Samples  for  Differential  Display  (DD)  analysis  were 
drawn  at  different  times  of  the  cell  cycle:  once  at  the  beginning  and  once  at  the  end  of  the 
Gi  cell  cycle  phase,  once  at  the  peak  of  the  S  phase  and  a  last  time  during  G2  (Fig.  1).  The 
mRNA  fractions  were  purified  and  analyzed  by  DD  as  reported  elsewhere  (Chapter  III). 

cDNA  Library  Screening.  DD  probes  for  screening  the  library  were  obtained  by 
labeling  the  cloned  cDNA  of  interest  with  Digoxigenin-ll-dUTP  (Genius  Nonradioactive 
Nucleic  Acid  Labeling  and  Detection  System,  Boehringer  Mannheim,  IN)  and 
hybridizations  were  conducted  overnight  as  previously  described  (Chapter  III).  Positive 
plaques  were  then  excised  in  vivo,  cloned  and  sequenced  following  the  manufacturer's 
protocols  (ZAP-cDNA  cloning  kit,  Stratagene). 

RNA  Analysis.  RNA  dot  blot  analysis  were  performed  as  described  elsewhere 
(Chapter  III).  For  the  RNA  Northern  analysis,  ten  pg  of  total,  DNAse  treated  RNA 
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Figure  1  Cell  cycle  dynamics  in  synchronously  growing  Alexandriumfundyense.  Lines 
show  relative  distributions  of  cells  in  the  different  cell  cycle  stages  (Gi,  S, 
and  G2  +  M,  Arrows  denote  sampling  points  for  the  Differential 

Display  analysis  protocol:  Two  during  Gi,  (I)  early  and  (II)  late  segments, 
one  during  the  S  phase  (III)  and  a  last  one  during  G2  (IV).  Shaded  areas 
correspond  to  the  dark  periods  of  the  14:10  hour  light:dark  cycle  under  which 
A.  fundyense  was  synchronized. 
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representing  each  of  the  four  cell  cycle  stages  described  above  were  loaded  on  a 
formaldehyde  gel,  run  at  a  constant  intensity  of  5  V/cm  in  MOPS  running  buffer  and 
transferred  overnight  to  positively  charged  nylon  membranes  (Boehringer  Mannheim,  IN). 
The  histone-like  specific  probe  obtained  from  the  DD  analysis  was  labeled  as  described 
above.  Hybridizations  were  performed  overnight  in  high  SDS  buffer,  as  recommended  by 
the  manufacturer,  and  at  42°C  (Genius  Nonradioactive  Nucleic  Acid  Labeling  and  Detection 
System,  Boehringer  Mannheim,  IN). 

Sequence  Analysis.  The  nucleotide  sequence  was  compared  to  the  non-redundant 
GenBank,  EMBL,  DDBJ  and  PDB  sequence  databases  using  the  Basic  Local  Alignment 
Search  Tool  program  (BLAST)  in  its  version  for  nucleotides  (BLASTN)  or  aminoacids 
(BLASTX)  (Altschul  et  al,  1990).  The  corresponding  amino  acid  sequence  was  predicted 
using  the  universal  genetic  code.  The  deduced  protein  sequence  was  further  characterized 
with  the  sequence  analysis  software  program  MacVector  (Oxford  Molecular  Limited). 
Secondary  structures  were  determined  by  the  Robson-Gamier  method  (Gamier  et  al, 
1978),  and  antigenicity  and  hydropathy  profiles  were  calculated  following  standard 
procedures  (Kyte  and  Doolittle,  1982;  Jameson  and  Wolf,  1988). 


RESULTS 

In  a  previous  study  targeting  the  identification  of  up-  and  downregulated  genes  during 
toxigenesis  in  the  toxic  dinoflagellate  Alexandrium  fundyense,  a  gene  for  a  putative 
histone-like  protein  was  isolated  (Chapter  III).  The  Differential  Display  pattern  of  this 
protein  showed  accumulation  of  the  corresponding  mRNA  at  the  beginning  of  Gi  (Fig. 
2A),  a  pattern  that  was  confirmed  by  total  RNA  dot  blot  analysis  (Fig.  2B)  and  by 
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Northern  blot  analysis  (Fig.  2C).  The  DD  probe  hybridized  to  a  band  corresponding  to  a 
size  of  approximately  600  base  pairs  (Fig.  2C).  The  corresponding  cDNA  rescued  from  an 
A.  fundyense  cDNA  library  contained  577  nucleotides  and  an  open  reading  frame  (ORF)  of 
255  nucleotides  or  the  equivalent  of  85  aa  (Chapter  III).  The  database  search  aligned  this 
sequence  to  two  histone-like  proteins  from  Crypthecodinium  cohnii  Biecheler.  The 
combined  simileirity  was  60%,  with  each  of  the  two  different  histone-like  proteins  showing 
a  similarity  of  52%  and  a  homology  of  35%.  No  significant  similarities  to  other  proteins 
were  found. 

The  inferred  protein  has  a  molecular  weight  of  approximately  9.3  kDa.  Amino  acid 
composition  and  other  characteristics  of  this  protein  were  analyzed  by  comparison  with 
three  histone-like  proteins  from  different  dinoflagellates  (Table  1;  Fig.  3).  Two  indicators 
of  basicity  were  calculated:  (1)  the  fraction  of  basic  amino  acids  (Arg,  Lys)  per  protein, 
which  was  remarkably  high  at  about  24%  and  conserved  among  the  four  proteins,  and  (2), 
the  ratio  of  basic  to  acidic  amino  acids  ((Arg,  Lys)/(Gln,  Glu,  Asp,  Asn)),  which  was  the 
same  for  HAf  and  one  of  the  Crypthecodinium  cohnii  proteins  and  approximately  2.5  times 
lower  for  the  second  C.  cohnii  protein  and  the  Oxyrrhis  marina  Dujardin  protein  (Table  1). 
Even  though  the  compositional  analysis  of  all  four  proteins  revealed  a  high  degree  of 
similarity  among  them,  a  few  differences  are  apparent:  The  Lys  to  Arg  ratio  shows  a  high 
inter-  and  intraspecific  variability  (Fig.  3;  Table  1);  two  cysteine  residues  are  clearly  present 
only  in  HAf  and  in  HCcl,  while  only  trace  levels  have  been  reported  for  HCc  and  HOm 
(Fig.  3);  the  high  percentage  of  methionine  residues  in  HCcl  corresponds  to  eight  residues 
on  the  protein  while  HAf,  HCc  and  HOm  only  contain  one  methionine  or  none.  In  a 
separate  analysis,  it  was  determined  that  in  contrast  to  the  two  deduced  protein  variants 
from  C.  cohnii  and  the  amino-terminal  of  a  third  peptide,  also  an  HCc  protein,  HAf  only 
contains  one  methionine  in  the  NH2-region,  while  the  other  three  peptides  contain  between 
three  and  six  methionine  residues  (Fig.  4). 
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Figure  2.  Expression  of  HAf.  (A)  Differential  Display  (DD)  analysis  of  Alexandrium 
fundyense  in  early  and  late  Gi,  S  phase  and  G2  (arrow  points  at  HAf; 
numerals  I  through  IV  refer  to  sampling  points  in  Fig.  1;  double  lanes 
represent  duplicate  RT-PCRs).  (B)  Confirmation  of  DD  pattern  by  total  RNA 
dot  blot  analysis.  (C)  RNA  blot  analysis  (Northern)  confirming  the  induction 
of  the  histone-like  protein  (HAf)  in  Alexandrium  fundyense  during  Gi  (RNA 
loading  was  evaluated  by  ethidium  bromide  staining  and  identification  of  the 
18S  and  28S  rRNA  bands  (bottom)). 
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Table  1 .  Properties^  of  histone-like  proteins  from  several  dinoflagellates 


HAf^ 

HCcl‘' 

HCc'^ 

HOm' 

M.W. 

9,319  Da 

12,000  Da 

12,532  Da 

— 

Basic/Total  aa^ 

23.5 

25.7 

23.5 

25.1 

Basic/acidic 

(K,RrtD,N,E,Q) 

4.3 

4.1 

1.6 

1.4 

Lys® 

9.4 

17.7 

17.8 

20.7 

Arg 

12.9 

8 

3.6 

2.3 

Lys/Arg 

0.7 

2.2 

5 

9 

a.  Determined  directly  (O.  marina  &  C  cohniiV)  or  from  deduced  amino  acid 
sequence  (A.  fundyense  8c  C  cohnii) 

b.  Alexandrium  fundyense  (this  study) 

c.  Crypthecodinium  cohnii  (Sala-Rovira  et.  at,  1991) 

d.  Crypthecodinium  cohnii  (Rizzo  &  Morris,  1984) 

e.  Oxyrrhis  marina  (Jing-Yan  1984) 

f.  Ratios  obtained  from  Mole%  of  the  different  aa 

g .  Numbers  represent  Mole% 
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Figure  3. 


Amino  acid  compositions  of  histone-like  proteins  from  several  dinoflagellates 
(Asx  and  Glx  refer  to  Asp  and  Asn,  and  Glu  and  Gin,  respectively). 
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MOLE  %  TOTAL  AMINO  ACIDS 


AMINO  ACID 
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Figure  4.  Alignment  of  the  NH2-terminals  of  HAf,  HCc2,  HCc  and  HCcl  (Stars  denote 
gaps;  Met  residues  are  underlined).  The  HAf,  HCcl  and  HCc2  peptide 
sequences  are  inferred  from  the  cloned  cDNA  and  HCc  was  determined 
directly  from  a  peptide. 
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A.  fundyense 
C.  cohniil 
HCc 

C.  cohniil 


MSPSAT 
MKAMKA*  *  *TKKi^T 
PKMKAAMKATMKi^KA*  *  *TKK;^T 
^PI^KAA*  *  *  *MK^KAPi^KGi™T 
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Figure  5.  Hydropathy  profiles  of  the  histone-like  proteins  HAf  (A),  HCc2  (B)  and 
HCcl  (C).  The  proteins  are  positioned  relative  to  each  other  according  to  their 
sequence  alignment.  Positive  values  indicate  hydrophylicity  and  negative 
values  hydrophobicity.  Numerals  I  through  V  indicate  regions  of  structural 
and/or  sequence  similarity  among  proteins  (see  text  for  details).  NH2  and 
COOH  denote,  respectively,  the  amino  and  carboxy  terminals  of  the  proteins. 
The  dotted  lines  frame  the  smallest  of  the  three  proteins,  HAf.  The  length  of 
the  black  bar  represents  20  aa  on  the  abscissa. 
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Predictions  of  the  secondary  structure  as  well  as  antigenic  index  and  hydropathy 
profiles  were  obtained  for  HAf,  HCcl(a)  and  HCc(b).  Five  distinct  regions  could  be 
defined.  These  domains  are  plotted  on  the  hydropathy  profiles  (Fig.  5).  The  hydropathy 
profiles  revealed  an  overlapping  pattern  of  alternating  hydrophylicity  and  hydrophobicity 
among  the  three  proteins.  Region  I  (HAf  aa  14-33)  is  characterized  by  high  sequence 
variability  among  the  three  proteins.  Structurally,  it  is  defined  by  a  succession  of  P-sheets 

and  the  highest  hydrophylicity  within  each  of  the  proteins.  Region  II  (Haf  aa  35-42)  is 
highly  conserved  and  of  varying  hydropathy  among  the  three  proteins.  It  represents  the 
transition  to  region  III  (HAf  aa  43-52),  a  very  hydrophobic  region  with  P-sheet 
configuration  in  the  middle  of  all  three  proteins.  In  all  three  cases  this  region  also  shows 
the  lowest  antigenicity  index  of  the  entire  protein.  Region  IV  (HAf  aa  54-65)  is  again 
highly  conserved  and  very  hydrophylic.  The  carboxy-end,  region  V  (Haf  aa  75-83),  is 
hydrophobic  and  has  a  low  antigenicity  index.  Region  V  is  constituted  of  an  a-helix.  The 
carboxy-terminal  portions  of  HCcl  and  HCc2  not  overlapping  with  HAf  consist  of 
hydrophylic  P-sheets. 


DISCUSSION 

Isolation  of  a  putative  histone-like  protein  in  Alexandrium  fundyense.  Differential 
Display  analysis  of  different  cell  cycle  stages  and  levels  of  toxigenicity  in  A.  fundyense 
resulted  in  the  identification  of  a  putative  histone-like  protein.  Differential  Display  is  a 
genetic  technique  that  allows  the  identification  of  transcriptionally  regulated  genes  in 
eukaryotes  (Liang  and  Pardee,  1992).  It  is  based  on  an  RT-PCR  (reverse  transcription  - 
polymerase  chain  reaction)  approach  and  the  usage  of  primers  that  specifically  amplify 
mRNA.  Differential  Display  is  highly  sensitive  and  reproducible  and  has  been  successfully 
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used  to  identify  developmentally  and  environmentally  regulated  genes  in  a  variety  of 
organisms  (Utans  et  al,  1994;  Abrahamsen  et  al,  1995;  Sharma  and  Davis,  1995;  van  der 
Knaap  and  Kende,  1995;  Wan  et  al,  1996). 

Toxigenesis  in  synchronized  cultures  of  Alexandrium  fundyense  is  circumscribed  to  a 
period  of  8  to  10  hours  at  the  beginning  of  the  Gi  phase  of  the  cell  cycle  (Taroncher- 
Oldenburg  et  al,  1997).  In  an  effort  to  identify  genes  that  were  up-  or  downregulated 
during  those  periods  of  toxin  biosynthesis,  we  analyzed  the  mRNA  populations  of  cells 
that  were:  (1)  in  Gi  and  producing  toxin;  (2)  in  Gi  and  not  producing  toxin;  (3)  in  the  S 
phase;  and  (4)  in  Gz  (Fig.  1)  (Chapter  III).  As  a  result  of  this  experimental  design,  we 
expected  to  isolate  not  only  cDNAs  of  genes  that  were  putatively  related  to  STX 
biosynthesis,  but  also  cDNAs  of  genes  that  might  bear  some  relation  to  the  beginning  of  the 
Gi  cell  cycle  phase  of  A.  fundyense.  One  of  the  bands  we  identified  rescued  a  full-length 
cDNA  with  a  similarity  of  60%  at  the  amino  acid  level  with  histone-like  proteins  from  the 
dinoflagellate  Crypthecodinium  cohnii  (Chapter  III). 

The  identification  and  characterization  of  this  putative  histone-like  protein  is  significant 
because  molecular  information  is  scarce  on  this  family  of  dinoflagellate  proteins. 
Comparative  biochemical  and  structural  analysis  of  the  A.  fundyense  protein  with  several 
histone-like  proteins  from  other  dinofiagellates  also  sheds  new  light  onto  our  understanding 
of  their  possible  structure-function  relationships. 

Characterization  ofHAf  Histone-like  proteins  are  defined  by  several  characteristics: 
Low  molecular  weight,  high  basicity,  low  protein  to  DNA  ratio  and  their  affinity  for  DNA 
(Rizzo,  1985).  They  share  some  of  these  properties  with  the  'true'  histones  but  still  form  a 
class  by  themselves  as  their  molar  ratio  with  DNA  is  typically  about  20%  or  less  (while  true 
histones  tend  to  be  in  equimolar  concentrations  with  DNA)  and  have  an  affinity  for  DNA 
that  is  significantly  less  than  that  of  histones  (Rizzo,  1985).  The  five  different  types  of 


199 


histones  characterized  in  eukaryotes  are  highly  conserved,  both  in  sequence  and  in 
structure,  while  histone-like  proteins  are  much  more  variable.  'True'  histones  have  very 
well  defined  histone-histone  binding  sites  that  are  not  present  in  histone-like  proteins 
(Rizzo,  1985).  Inter-  and  intramolecular  binding  sites  in  histone-like  proteins  have  been 
suggested  to  be  related  to  the  presence  of  cysteine  residues  in  the  molecule  (Vemet  et  al, 

1990) . 

Histone-like  proteins  from  several  dinoflagellates  have  been  characterized  using 
electrophoresis,  peptide  mapping  and  amino  acid  analysis  (Jing-Yan,  1984;  Rizzo  and 
Morris,  1984)  (reviewed  Rizzo,  1987).  Only  two  complete  histone-like  protein  sequences 
and  one  amino-terminal  peptide  have  been  reported  in  the  literature  (Sala-Rovira  et  al., 

1991) .  A  remarkably  conserved  characteristic  of  these  proteins  is  their  high  basicity 
(approximately  24%)  which  is  responsible  for  their  DNA  affinity  (Table  1).  Other 
characteristics  such  as  the  ratio  of  basicity  to  acidity  show  a  wide  range  and  tell  us  more 
about  the  kind  of  variations  to  expect  among  this  class  of  proteins.  For  instance,  the  two 
C.  cohnii  proteins,  HCcl  and  HCc,  show  a  high  disparity  in  this  latter  parameter  and 
exhibit  respectively  more  similarity  to  HAf  or  HOm  (Table  1).  The  amino  acid 
compositions  are  a  further  example  for  inter-  and  intraspecific  variations  superimposed  over 
a  pattern  of  analogous  functionalities  (Fig.  3).  It  is  of  interest  to  notice  the  antiparallel 
distributions  of  the  main  basic  amino  acids,  Arg  and  Lys  (Fig.  3;  Table  1).  These  patterns, 
resulting  in  the  highly  conserved  basicity  of  the  protein,  the  main  characteristic  of  these 
molecules,  suggest  a  low  evolutionary  pressure  on  the  actual  composition  of  the  proteins 
but  a  very  high  pressure  on  their  functional  basicity.  Comparative  studies  of  their  peptide 
maps  and  electrophoretic  mobilities  have  also  shown  that  these  proteins  are  less  conserved 
than  the  equivalent  set  of  histones  in  other  eukaryotes  (Rizzo  and  Morris,  1984). 

The  molecular  weights  of  histone-like  proteins  range  from  approx.  10  kDa  to  about  17 
kDa  with  two  exceptions  in  C.  cohnii,  where  in  addition  to  six  low  molecular  weight 


proteins,  two  basic  proteins  of  high  molecular  weight  (45  and  50  kDa)  have  been 
characterized  (Vemet  et  al,  1990).  The  deduced  amino  acid  sequence  we  report  in  this 
paper  gives  a  molecular  weight  of  approximately  9.3  kDa  to  HAf,  thereby  becoming  the 
smallest  such  protein  described  in  the  literature.  Insufficient  sequence  information  is 
available  to  allow  for  a  comprehensive  analysis  of  size,  structure  and  functionality  of  these 
proteins,  but  one  interesting  issue  was  raised  by  the  analysis  of  the  amino  terminals  of  the 
two  variants  HCcl  and  HCc2  (Sala-Rovira  et  al,  1991).  The  high  content  in  methionine 
residues  observed  in  the  fost  ten  to  twenty  aminoacids  suggested  the  possibility  of  multiple 
or  alternative  transcription  initiation  sites  that  could  cause  significant  variations  in  the  size 
of  these  proteins  (Fig.  4).  In  contrast,  the  amino-terminal  of  HAf  only  contains  one 
methionine  residue  that  aligns  with  aa  17  of  HCcl  and  aa  6  of  HCc2  (Fig.  4).  Synthesis  of 
peptides  with  different  sizes  from  the  same  mRNA  template  has  been  recently  described  for 
the  luciferin-binding  protein  from  Gonyaulax  polyedra  (Mittag  et  al,  1997).  In  this  case  it 
was  shown  that  the  synthesis  of  two  LBP  variants  of  different  length  was  due  to 
translational  initiation  at  two  different  AUG  codons  (Met  residues)  located  within  the  5'- 
end  of  the  same  ORF.  Such  a  mechanism  could  be  involved  in  the  synthesis  of  the 
different  HCc  variants  that  have  been  biochemically  characterized  in  Crypthecodinium 
cohnii  (Vemet  etal,  1990;  Sala-Rovira  etal,  1991). 

Despite  their  size  differences,  the  three  proteins  exhibit  remarkably  similar  hydropathy 
and  antigenicity  profiles  (Fig.  5).  The  core  stmcture  consists  of  one  small  and  highly 
hydrophobic  region  flanked  by  two  larger  and  veiy  hydrophylic  domains  This  overlap  in 
molecular  characteristics  might  relate  to  the  function  of  these  proteins.  More  sequenee 
information  of  other  histone-like  proteins  is  needed  to  determine  if  indeed  their  DNA 
binding  properties  are  related  to  that  specific  configuration.  Predictions  of  the  secondary 
structures,  though  consistent,  should  still  only  be  considered  suggestive  of  the  trae 
stmctures  because  of  a  lack  of  supporting  stmctural  data  (Nishikawa,  1983). 
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Expression  of  HAf.  Immunocytochemical  studies  to  elucidate  the  distribution  of  the 
DNA-binding  protein  HCc  during  the  cell  cycle  of  Crypthecodinium  cohnii  revealed  a  close 
association  between  the  protein  and  nucleofilaments  as  well  as  extrachromosomal  loops  of 
DNA  (Geraud  et  al,  1991).  Constant  labelling  throughout  interphase  as  well  as  the 
different  stages  of  mitosis  suggests  new  supply  of  protein  at  either  the  very  end  of  the 
mitotic  process  or  the  beginning  of  Gi  (Geraud  et  al,  1991).  Our  results  with  DD  and 
their  confirmation  by  Northern  blotting  (Fig.  2C)  show  clear  induction  of  HAf  expression 
at  the  start  of  the  Gi  phase  of  the  cell  cycle,  corroborating  the  prior  observations  on  C. 
cohnii.  Other  cell  cycle-dependent  genes  have  been  shown  to  be  transcriptionally  regulated 
in  eukaryotes  (Price  et  al,  1991).  Some  histones  are  also  induced  in  proliferating  cells 
(Callard  and  Mazzolini,  1997).  We  conclude  that  histone-like  proteins  are  closely 
associated  with  DNA  and  their  expression  is  regulated  in  a  cell  cycle  dependent  way. 
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